BWA-MEM Vs BWA-ALN
0
1
Entering edit mode
6.5 years ago
SOHAIL ▴ 400

Hi everyone,

I'm sorry if it is answered anywhere else, I have raw dna-seq paired-end sequence reads data from different libraries of read length for different samples. Ranging from 50-90, 50-101, 50-150 (after quality trimming). In FASTQC check i noticed that about 10% of reads are lesser than 70 bp for each sample.

As BWA-mem is pretty mature now but the bwa docs say that bwa mem is preferable for longer reads ( > 70 bp). and BWA-aln is good for shorter reads, but i guess most of the reads from my library are not shorter.

I want to use one mapping algorithm for consistency reasons.Can i still use bwa-mem ignoring small fraction of smaller read lengths? or i should use BWA-aln ?

Thanks

sohail

next-gen genome alignment • 9.4k views
ADD COMMENT
2
Entering edit mode

bwa-mem could be used for short reads as well, but this blog suggests that bwa-aln performs better for 50bp reads but over all mem has a good accuracy at different mismatch levels. Would it be useful if you align the 10% of the short reads with mem and aln and compare the stats to see if it makes a big difference.

ADD REPLY
0
Entering edit mode

yes, i think i should try to compare..

ADD REPLY

Login before adding your answer.

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