text: **Imputation of ancient genomes.**We follow the recommended pipeline 
  of GLIMPSE73 and first call genotype likelihoods for each genome in 
  the 1000GP, segregating sites using bcftools mpileup with filter -q 20, 
  -Q 20 and -C 50. We subsequently impute each genome separately using 
  GLIMPSE v. 1.1.1 using the 1000GP phase 3 reference panel74 downloaded 
  from https://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20130502/. 
  These imputed genomes are merged into a single VCF (variant call 
  format) for further downstream processing. 
  when migration rates are high (and, therefore,*F*ST is low), indicating 
  that Twigstats is better powered to detect such scenarios of continued 
  migration. Encouragingly, a model that involves the two immediately 
  adjacent populations is selected in all replicates as the ‘best’ model 
  (highest qpAdm*P*value) using Twigstats, whereas this is the case in 
  only 80% (migration rate of 0.001) and 30% (migration rate of 0.005) 
  of all replicates using regular qpAdm. 