gargammel: a sequence simulator for ancient DNA.

Abstract

Ancient DNA has emerged as a remarkable tool to infer the history of extinct species and past populations. However, many of its characteristics, such as extensive fragmentation, damage and contamination, can influence downstream analyses. To help investigators measure how these could impact their analyses in silico , we have developed gargammel, a package that simulates ancient DNA fragments given a set of known reference genomes. Our package simulates the entire molecular process from post-mortem DNA fragmentation and DNA damage to experimental sequencing errors, and reproduces most common bias observed in ancient DNA datasets. The package is publicly available on github: https://grenaud.github.io/gargammel/ and released under the GPL. gabriel.renaud@snm.ku.dk. Supplementary data are available at Bioinformatics online.