What could make Bowtie2 miss legitimate hits?
0
1
Entering edit mode
18 months ago
iddo.nadav ▴ 10

Hi,

We've noticed some strange behavior by Bowtie2: In specific cases - it simply misses a proper hit for a read. We noticed this after investigating unmapped reads - we were surprised to see that they have a perfect end-to-end match to the genome we were aligning them against. This raised the question - why did Bowtie fail to catch this hit.

Repeating the process showed that these were the same reads that bowtie was failing with (so this is not a random error rate). These reads are trimmed, contain no adaptor, have high quality, and proper length, and - as mentioned above - have a perfect match in the genome.

Any idea why Bowtie2 is missing their alignemnt?

Thanks, Iddo

alignment • 279 views
ADD COMMENT
0
Entering edit mode

One of a set of PCR duplicates? Can you filter them out with picard tools? What does the bam flag say? Did you try a different aligner, e.g. BWA mem?

ADD REPLY
0
Entering edit mode

Yes, we did. Surprisingly BWA was able to map these reads. How would you interpret this?

ADD REPLY

Login before adding your answer.

Traffic: 2007 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6