Bowtie2 and bwa *big* disparity in alignment percentages
1
0
Entering edit mode
2.5 years ago

Dear all, I have aligned my 300*2 data against a 4 Mb reference using both Bowtie2 and bwa-0.7.

The alignment percentage observed in bowtie2 is 80% and with bwa it is 99%. What is fundamentally different causing such a large difference between the aligners.

The coverage once calculated from the mpileup drops tp 1% at 20x. So, the depth is low in both the methods. But bwa has a better horizontal coverage in terms of alignment but depth drops in this too.

Thanks

alignment • 556 views
ADD COMMENT
0
Entering edit mode
2.5 years ago

This is because by default the two aligners are set up to run differently.

To produce equivalent alignments try:

bowtie2 --very-sensitive-local

or if feeling adventourous:

bowtie2 -D 20 -R 3 -N 1 -L 20

Beyond that the changes run deeper than just how many alignments are produced. I have a chapter on this in the Biostar Handbook (not a free resource): https://www.biostarhandbook.com/align/compare-aligners.html

ADD COMMENT
0
Entering edit mode

Sure I'll try this out!

ADD REPLY

Login before adding your answer.

Traffic: 2253 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