The paired-end contemporary DNA sequences were mapped to the human reference genome (GRCh37/hg19) or the revised Cambridge reference mitochondrial genome (32) using BWA version 0.6.2 (33). Duplicate sequences were then removed from Alignment/Map (SAM) files using SAMtools (35, 36). Using the mpileup command, with a maximum depth parameter of 1000, variant call format (VCF) files were generated for each chromosome separately. Using an awk command, indel variations were excluded. The VCF files of individual modern genomes were merged using VCFtools (37) after zipping and indexing using tabix.

Note: The content above has been extracted from a research article, so it may not display correctly.

Please log in to submit your questions online.
Your question will be posted on the Bio-101 website. We will send your questions to the authors of this protocol and Bio-protocol community members who are experienced with this method. you will be informed using the email address associated with your Bio-protocol account.

We use cookies on this site to enhance your user experience. By using our website, you are agreeing to allow the storage of cookies on your computer.