Questions tagged [sequencing]
DNA sequencing is the process of determining the precise order of nucleotides within a DNA molecule.
152
questions
2
votes
0
answers
55
views
Adapter trimming
I am trying to do adapter trimming, alignment and sorting for a range of large scale paired end fastq files. The code I am using is given below:
...
2
votes
0
answers
26
views
What data set was used Gene Expression Data of Postmortem Tissues?
I want to run the experiments "Sample Application to Genotype-Tissue Expression (GTEx) Project Gene Expression Data of Postmortem Tissues" mentioned in "Enter the Matrix: Factorization ...
- 477
1
vote
1
answer
36
views
Lower case vs. upper case nucleotids in sequence vs. dots at the end
What is the difference between lower case and upper case nucleotides in a sequence? My other question is what are the dots at the end of the sequence? Some examples are shown below:
GGgG,GGgG,GGgG,...
3
votes
0
answers
29
views
Can multiplexing in Sequel II SMRTcells reduce the coverage?
I would like to have high depth of coverage of the transcriptome, but I need to analyze several tissues. I was suggested to pool all 5 tissue samples in same SMRT 8M cell. Will this result in a ...
- 247
1
vote
1
answer
31
views
How can Iso-Seq reverse transcriptase artifacts be avoided?
My end goal is annotate a de-novo assembled genome. When trying to select the best method for transcriptome assembly I read Iso-seq was the preferred method. However other people suggested Nanopore as ...
- 247
1
vote
2
answers
98
views
What is the most accurate approach for de Novo sequencing?
I'm trying to decide between PacBio HiFi or Illumina sequencing platforms for sequencing the genome of aChrysina scarab. We want to identify the color pattern loci and also perform a phylogenetic ...
- 247
2
votes
3
answers
75
views
Can we do NGS library preparation using UMI with large amount of DNA input?
We know that in next-generation sequencing (NGS), the unique molecular identifier (UMI) can reduce or eliminate sequencing or PCR errors and result in very high accurate data. Therefore, UMI is widely ...
- 21
1
vote
1
answer
50
views
How can I generate ASV file from nanopore sequencing data?
I am converting the fast5 file to fastq by using guppy basecaller after that by using kraken2 classified the sequence.
Now I am trying to generate ASV file.
Is this possible to generate ASV file from ...
- 21
1
vote
1
answer
60
views
Extracting base sequences from ABI/AB1 sanger sequencing chromatogram
I am trying to understand the sanger sequencing ABI/AB1 file format better, and extract base calls from given signal intensities over time.
As I understand, reading in a raw AB1/ABI file into python, ...
- 11
1
vote
1
answer
68
views
Does the number of RNA reads per cell obtained from the 10X scRNA experiment depend on amount of mRNA in given cell?
As we know, the amount of RNA reads per cell obtained from 10X scRNA experiment vary between cells. I wonder if this is effect of technical issues or does the number of RNA reads per cell obtained ...
- 11
0
votes
0
answers
62
views
Interpreting 'samtools mpileup' output for multiple inputs
I would like to calculate sequencing coverage for a WGS project. Both long and short reads. I've used samtools as following:
...
- 627
0
votes
1
answer
51
views
Why must a maximal non-branching path be a contig?
The following is from Bioinformatics Algorithms:
Fortunately, we can derive contigs from the de Bruijn graph. A path in a graph is called non-branching if in(v) = out(v) = 1 for each intermediate ...
- 127
0
votes
0
answers
41
views
Locate the saved template of minknow
I am using minknow on ubuntu. After saving the template (settings), where is that file locally stored in the system?
Also what is its format/extension?
- 1
1
vote
1
answer
40
views
What are Dummy fastqs?
I had to do a file transfer for in our system to run from fastqs and I came across the term "dummy fastq", not sure what exactly it means and what is the purpose of it in the workflow. Can ...
- 19
1
vote
2
answers
121
views
How to retrieve fasta sequence after local blast?
I have created a Blast database using a reference genome. Then, I have performed a local blast search in command line using a gene of interest. I have obtained some hits with the usual Blasting ...
- 75
0
votes
0
answers
27
views
Drug-DNA complex simulation using AMBER
I want to carry out the simulation of drug-dna complex by placing the drug molecule at particular site in my DNA sequence (intercalating site). I have PDB files of both drug and DNA sequence. I will ...
2
votes
1
answer
117
views
When using cutadapt, can I specify an R2 adapter as optional when I have a required R1 adapter?
I'm using cutadapt 3.5 to trim adapters and perform some filtering on paired-end data.
Both R1 and R2 sequences have 3' adapters that might be found depending on the sequence length, but R1 also has ...
- 671
0
votes
1
answer
67
views
Shotgun analysis: Where and how to start?
I am a complete rookie in shotgun analysis, having only performed marker gene analysis before. I have a pipeline that I want to try but I can't wrap my head around how it all works.
Please correct me ...
- 53
0
votes
1
answer
95
views
How to Read SCF file in Python?
Is there any way that I can read SCF file in python like in R using sangerseqR, I have tried with Biopython, it seems it does not support this format.
- 151
1
vote
1
answer
54
views
Where to find publicly available Sanger chromatography data?
I am looking for any public available databases where I can download Sanger chromatography data ideally in Ab1 or SCF file format. I need extremely large amounts of this, since I want to use it for ...
- 151
2
votes
1
answer
79
views
Low pass sequencing has been reported to detect common variants. How low can one go and get reliable data? Is 2X pass sequencing analysis possible?
I would like to use low pass sequencing to replace a genotyping chip to be able to detect variants up to 0.1 % allele frequency in available population data. What is the minimum depth I can opt for to ...
- 21
1
vote
0
answers
202
views
Illumina vs PacBio vs Nanopore price
Can someone help me with comparing the cost per run of Illumina, PacBio and Nanopore sequencing technologies?
- 31
0
votes
1
answer
42
views
Which format is the most raw format
I wanted information regarding which format is considered the most raw from an illumina sequencer ?
fasta
fastq
bcl
bam
As per my research it should be bcl but I am not sure.
- 11
0
votes
1
answer
49
views
Which sequencing technologies are considered short read
I couldn't find any information on this online. Which of the following sequencing technologies are considered short read technologies? Illumina, Oxford Nanopore, Ion Torrent, Pacific Biosciences
- 19
2
votes
2
answers
239
views
Construct the Overlap Graph of a Collection of k-mers
I am doing a course on graph algorithms in genome sequencing and the current assignment involves building an overlap graph from a bunch of k-mers. Formally the problem is as follows:
Input: A ...
- 21
2
votes
1
answer
339
views
What is the right way of calculating a Phred score by hand?
i am trying to calculate mean Phred scores for my sequencing data, but i feel not very comfortable about it. There are actually two ways of calculating. (I just use an existing sample)
giving: 3 reads ...
- 45
2
votes
2
answers
147
views
Get frequency of index sequences of fastq file
I recently did some sequencing on a MiniSeq and unexpectedly the Undetermined_S0_L001_R1_001.fastq.gz file contained quite a lot of data (about 40% of total). I ...
- 281
0
votes
1
answer
37
views
Modeling number of reads mapped to a gene
I am looking for a probability distribution of a number of reads mapped to a particular gene in metagenomic sequencing (NGS, shotgun, likely illumina).
Naively one could model it via a binomial (or ...
- 357
-1
votes
2
answers
111
views
Filter rows of VCF file for Match=EXACT?
How do I remove rows in a VCF file on 1kGenome column, where Match=EXACT using bash ?
- 11
1
vote
1
answer
114
views
What is a PacBio "movie file"?
I came across references to "movie files" from PacBio sequencing in this paper:
https://www.jimmunol.org/content/204/12/3434
Specifically:
Movie files used to generate results presented in ...
- 671
0
votes
1
answer
52
views
How to determine a correct promoter sequence? (Bisulfite sequencing)
I'd like to ask you for help with an university assignment. First let me post the entire question:
Design primers for bisulfite sequencing of the promoter of the ACE2
gene in microbat Myotis ...
1
vote
1
answer
60
views
Quast duplication ratio and mismatches percent
When analyzing Quast results it seems that it doesn't calculate mismatches and indels in a useful way if the "Duplication ratio" is over 1.
For example, that's what I get for an assembly ...
- 349
6
votes
1
answer
83
views
How to search for high coverage SRA entries
The Question
I want to find high coverage SRA entries, e.g., above 100x.
I guess the best way is to use https://www.ncbi.nlm.nih.gov/sra with an appropriate search term. I don't mind if the search ...
- 261
1
vote
1
answer
80
views
Metagenome simulation from a concatenated FASTA file
I am trying to simulate metagenome sequencing. So I will start with a file with a lot of concatenated genomes. From there, I would like to randomly extract 10,000 sequences of length 200bps.
I got ...
- 41
1
vote
1
answer
118
views
0
votes
0
answers
26
views
DNA sequence error annotation
Suppose I generate a DNA sequence of the following pattern.
AAGTC
And after being passed through a channel with insertion, deletion and substitution errors obtain the following sequence
AAAGGC
Where ...
2
votes
2
answers
178
views
How can I subset WGS data to the level of WES variants?
I would like to compare mutational signatures1 in patients from different studies, however some studies are based on exome seq (i.e. ~20,000 coding variants) and some are from whole genome seq (i.e. ~...
- 203
2
votes
2
answers
182
views
Are sequencing reads written differently when inserts and indexes are sequenced seperately?
I was going over different sequencing library protocols and noticed that the adapter sequence can vary between protocols. In some adapter sequences the index sits between the primer binding site and ...
- 21
1
vote
0
answers
28
views
Fastq dataset after clustering
I'm doing a project for DNA archival storage systems
Is there an open-source dataset used as ground truth for clustering?
E.g
A file containing the ground truth strand sequences, which have been ...
- 111
3
votes
2
answers
90
views
Virus genome: protein inside another protein
The genome sequence for the GenBank virus JX869059.2 (a human betacoronavirus) contains, among others, two proteins:
N protein (range 28566..29807) and
Orf8b (range 28762..29100).
The interval for ...
- 215
1
vote
2
answers
99
views
Bacterial DNA at the tail of transcriptome reads. What does that mean?
I am assembling a transcriptome obtained from the Internet. The transcriptome was extracted from a human cancer tissue that had been previously grafted into a mouse. I have detected that many ...
- 349
0
votes
1
answer
117
views
How to anonymise a bam file
I have a bam file of, for instance, RNA-Seq, which contains patient-identifiable data in the form of Single Nucleotide Polymorphisms (SNPs) throughout.
I would like a technique to take this aligned ...
- 121
1
vote
0
answers
47
views
How many markers or loci are used in forensic DNA profiling and the effect of the DNA sequencing?
I am reading the Wikipedia page on DNA profiling, specifically STR analysis. It says that DNA profiling relies on matching 20 or so loci. It also gives an account of finding false identification. But ...
- 189
2
votes
4
answers
210
views
Separate multiple sequence into different file, one sequence per file
I have a file contain multiple sequence, and I want to separate them by "gene:" into different file.
example:
example.fa
...
- 33
1
vote
2
answers
128
views
How to combine two Genome-wide Association Study (GWAS)? [closed]
I did a GWAS analysis in the past for antibiotic resistance of E. Coli and the results were interesting (matching the literature). I did a new GWAS analysis for some new samples, but the results are ...
- 107
1
vote
1
answer
115
views
merge the matrix from computematrix of deeptools
I have several matrix from computematrix of deeptools. I need to merge two of them using "computeMatrixOperations cbind -m input1.mat.gz input2.mat.gz -o output.mat.gz"
but I am running to error "/...
- 95
2
votes
1
answer
61
views
Assembly reads having a copy of their beginning in their tail
I am analyzing the reads for the SARS-CoV-2 assembly with id SRR11140748. Apparently these reads were obtained with parallel sequencing by Illumina and Oxford Nanopore Technologies.
I have found these ...
- 349
2
votes
1
answer
148
views
What is the meaning of these misaligned reads in a sequencing run?
I am analyzing some SARS-CoV-2 sequencing runs abd often find read alignments like the one in the image.
...
- 349
2
votes
2
answers
872
views
Marking optical or PCR duplicates with picard vs. samtools flagstat
I am trying to identify the best library prep method for noninvasive prenatal test samples, to be sequenced with Illumina Novaseq.
One metric that I am evaluating is the number/ percentage of PCR or ...
- 243
3
votes
1
answer
96
views
Overlap in Paired-end Reads for Sequencing?
Regarding Sequencing:
Do/can the paired-end reads have overlaps?
- 107