Question: BWA mem and BWA aln lead to different results
0
gravatar for zhangdezhi008
11 weeks 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 11 weeks ago by dariober8.4k • written 11 weeks ago by zhangdezhi0080

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

ADD REPLYlink written 11 weeks ago by genomax40k
2
gravatar for dariober
11 weeks ago by
dariober8.4k
Glasgow - UK
dariober8.4k 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 11 weeks ago • written 11 weeks ago by dariober8.4k

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 10 weeks 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: 1575 users visited in the last hour