“…SNPs were called on the chromosomal Y. pestis mapped reads using the UnifiedGenotyper in GATK (DePristo et al, 2011). GLZ001 and GLZ002 were analyzed alongside a dataset of modern Y. pestis genomes (n = 233) (Cui et al, 2013;Eroshenko et al, 2017;Kislichkina et al, 2015Kislichkina et al, , 2018aKislichkina et al, , 2018bKutyrev et al, 2018;Morelli et al, 2010;Zhgenti et al, 2015) (as listed in (Keller et al, 2019)). In addition, we included previously published historical Y. pestis genomes with R 10-fold coverage (Bos et al, 2011(Bos et al, , 2016Feldman et al, 2016;Keller et al, 2019;Namouchi et al, 2018;Spyrou et al, 2016Spyrou et al, , 2019, and all previously published Late Neolithic and Bronze Age genomes with R 3-fold coverage (Andrades Valtueñ a et al, 2017;Rasmussen et al, 2015b;Spyrou et al, 2018) Article SNPs were called using the 'EMIT_ALL_SITES' option in GATK, which generated a call for every position in the chromosomal reference genome (CO92).…”