how to determine the percentage of reads after sequencing?
After sequencing I got two Fastq files with forward and reverse reads. I would like to know what percentage of reads belong to the virus I'm interested in. I mapped (bowtie2) these reads to the reference genome I used in the assembly and got 2%. I tried another method, did a metagenomic analysis (Kaiju), after which I got 30%. Why are these numbers so different? and which method should I choose and why.
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Assembly at lineage level from shotgun data ?
written by v.berriosfarias •
-
mapped paired end reads extraction from .bam file
written by t.ali •
-
Reads mapped on wrong strand
written by qstefano •
-
How to obtain my forward and reverse reads from sam / bam file
written by Bioinfo •
-
Low value of mapped reads against a metagenomic co-assembly
written by Biok •
-
How to make a paired-end reverse file
written by fb_marinho88 •
-
joining paired end illumina data for metagenomic analysis
written by cvu
-
How to Separate a SAM File Containing Forward and Reverse Reads into Two Separate Files
written by kdevakan •
-
Bowtie2: missing forward reads after alignment
written by nikkihathi •
-
Metagenome Assembly Using Combined Reads Through Flash And Using Abyss For Assembly
written by Sanjiv Kumar •
Is this a metagenomic dataset? What do you mean by
the reference genome I used in the assembly?If you have a genome available for the virus you are interested in then aligning the data to it is going to be a reasonably good way of identifying those reads. Depending on where the sample comes from there will be a chance of some misassignments when doing alignment to just viral genome (false positive and negative).