Bowtie2 mapping algorithm
Hi all,
Could you help me to understand the following problem? I am mapping reads to a viral reference genome, including different genotypes (someone very similar).
This is my actual strategy:
- First mapping allow me to find the best genome.
- Second mapping is against the choosen one, in order to increment the reads mapped.
This works, however sometime I get weird results.
For instance the total reads count is always greater but in some regions there are less reads!
How can be possible? Bowtie2 algorithm takes into account the reference length?
Many thanks, Diego
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Mapping to mtDNA and then align the unmapped
written by DanielEB_fisk •
-
Artificial reads - remove multiple mapped reads against reference genome, and only keep reads that …
written by Jalil Sharif •
-
Removal of host sequences without reference genome
written by emiliomastriani •
-
Analyzing differential expression of host genes in RNA-Virus infected cells?
written by DIG •
-
Bowtie2 reference genome
written by Ada •
-
How to visualize a pair reads mapped genome.
written by K.Gee •
-
mapping unmapped reads to viral genome using blast+
written by Shaway •
-
Viral de-novo Variant Call
written by Ribas •
-
Whole Genome Sequencing (WGS) Alignment Quality check softwares/tools
written by ravi.uhdnis
-
What Better Way To Get The Paired Reads Aligned Against The Reference Genome?
written by thiagomafra
You should include the commands you use, and illustrate your post with a particular example of the puzzling coverage issue. Did you check the mapping sam flags? Maybe it seems less reads because there are less supplementary / secondary mapping in some regions?
I am running bowtie2 with --very-sensitive-local alignment. I do not think is regarding supplementary/secondary reads since keeping only primary one.