What could make Bowtie2 miss legitimate hits?

1

Entering edit mode

4.5 years ago

iddo.nadav ▴ 30

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 • 753 views

ADD COMMENT • link 4.5 years ago by iddo.nadav ▴ 30

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 • link 4.5 years ago by Carambakaracho ★ 3.3k

0

Entering edit mode

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

ADD REPLY • link 4.5 years ago by iddo.nadav ▴ 30

Login before adding your answer.