Hi,

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

http://www.biomedcentral.com/1471-2164/13/178


**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
I had really good experiences with `bwa mem` when using it on dataset with high error rates. It far outperformed other tools. See this:

https://www.biostars.org/p/117225/#117259

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.