Question: BWA mem and BWA aln lead to different results
0
gravatar for zhangdezhi008
13 months ago by
zhangdezhi0080 wrote:

Hi, all,

I tried two uncover the genetic structure for a bird population using resequencing data, with PAIRED-END 100 bp reads. I used two different mapping methods BWA mem and BWA aln, and then used SAMtools to get vcf file. After filtering, we get two different SNP datasets when using BWA mem and BWA aln, and of course, this leads to two diffrent genetic structures for the same population! The problem is the mapping methods BWA mem and BWA aln, which one should I trust?

Any suggestions would help a lot! My best wishes, Dezhi

ADD COMMENTlink modified 13 months ago by dariober9.7k • written 13 months ago by zhangdezhi0080

Take a look at associated links: C: BWA aln vs. mem for 100bp reads

ADD REPLYlink written 13 months ago by genomax59k
2
gravatar for dariober
13 months ago by
dariober9.7k
Glasgow - UK
dariober9.7k wrote:

this leads to two diffrent genetic structures for the same population!

Provided that there are no bugs in your analysis pipeline, I think the question is not so much which aligner is better but rather on whether there is an actual population structure or whether the two different structures are actually the same.

I mean, bwa mem seems to be the standard now and it has superseded bwa aln (so use mem I would suggest). Still, bwa aln should give good alignments and if your final results depend on the choice between two good aligners, then these results are not very robust and should be taken with care,

ADD COMMENTlink modified 13 months ago • written 13 months ago by dariober9.7k

Thanks for your reply very much! According to the previous study, this species should have population structure with two lineages based on sanger sequences, which is consistent with the result generated from BWA-aln. But as you said, BWA-mem should be the better solution for mapping. Confused still...

ADD REPLYlink written 13 months ago by zhangdezhi0080
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 579 users visited in the last hour