Question: BWA MEM and ancient DNA
gravatar for smilefreak
5.2 years ago by
New Zealand
smilefreak420 wrote:


I have been doing some ancient DNA bioinformatics for a while now and have been using the bwa aln command. 

When using this command with ancient DNA it is recommended to disable seeding. This slows down alignments somewhat but allows the damaged DNA sequences to be aligned, the characteristic damage is a C->T transition at the 3' end of the sequenced fragments.  For instance I use the parameters -l 1024 (to disable seeding) -o 2 (increase gap open slightly) and -n 0.03 (reduce the edit distance slightly). These parameters were taken from this paper

My question:

Are there any specific settings required to use bwa mem with ancient DNA?
Does bwa mem use seeding at all? 

Thanks in advance

alignment adna • 3.2k views
ADD COMMENTlink modified 5.0 years ago by usmanroshan10 • written 5.2 years ago by smilefreak420
gravatar for Istvan Albert
5.2 years ago by
Istvan Albert ♦♦ 82k
University Park, USA
Istvan Albert ♦♦ 82k wrote:

I had really good experiences with bwa mem when using it on dataset with high error rates. It far outperformed other tools. See this:

A: When and why is bwa aln better then bwa mem?

Your problem seems to be very similar in nature in that it introduces mismatches that are similar to sequencing errors. From what I understood from the paper bwa mem was designed to be robust to this type of errors.

ADD COMMENTlink modified 5.2 years ago • written 5.2 years ago by Istvan Albert ♦♦ 82k
gravatar for usmanroshan
5.0 years ago by
usmanroshan10 wrote:

Hi, you may want to try the GPU program MaxSSmap that is designed for aligning divergent reads to genomes. used it to map ancient DNA reads and compared to BWA MEM it produced many more high quality alignments. You may also want to look into NextGenMap which also achieved a higher mapping rate BWA MEM.

ADD COMMENTlink modified 5.0 years ago • written 5.0 years ago by usmanroshan10
Please log in to add an answer.


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