Genotyping of newly sequenced as well as previously reported modern individuals was carried out using SAMtools/bcftools (35), followed by filtering, as previously described (44). Briefly, per-individual diploid genotypes were called using SAMtools mpileup (-C50 option) and bcftools call with the consensus caller (-c option). Calls from each genome were then filtered, excluding calls with low (six reads or one-third of the average sequencing depth, whichever was higher) or high (>2 times the average depth) coverage. We further filtered variants called within 5 base pairs of each other, with a Phred posterior probability of < 30 or a strand bias or end distance bias P value of <104, or with deviations from Hardy-Weinberg equilibrium with a P value of <104 across all samples.

For population genetic analyses, those genotypes were merged with a reference panel of 2286 modern individuals from worldwide populations, genotyped at 593,610 SNPs using the Affymetrix Human Origins array (26, 27). All ancient individuals were represented by pseudo-haploid genotypes obtained by sampling a random allele at each SNP position. For the two ancient samples with higher coverage, KP1 (6.9×) and MH8 (6.8×), we also carried out diploid genotyping, as described above, to be used in the haplotype sharing analyses.

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.