mapDamage: testing for damage patterns in ancient DNA sequences

Research output: Contribution to journalJournal articleResearchpeer-review

Ancient DNA extracts consist of a mixture of contaminant DNA molecules, most often originating from environmental microbes, and endogenous fragments exhibiting substantial levels of DNA damage. The latter introduce specific nucleotide misincorporations and DNA fragmentation signatures in sequencing reads that could be advantageously used to argue for sequence validity. mapDamage is a Perl script that computes nucleotide misincorporation and fragmentation patterns using next-generation sequencing reads mapped against a reference genome. The Perl script outputs are further automatically processed in embedded R script in order to detect typical patterns of genuine ancient DNA sequences. Availability and implementation: The Perl script mapDamage is freely available with documentation and example files at https://geogenetics.ku.dk/all_literature/mapdamage/. The script requires prior installation of the SAMtools suite and R environment and has been validated on both GNU/Linux and MacOSX operating systems.
Original languageEnglish
JournalBioinformatics
Volume27
Issue number15
Pages (from-to)2153-2155
Number of pages3
ISSN1367-4803
DOIs
Publication statusPublished - 1 Aug 2011

ID: 33910533