Integrated snATAC-seq data were split into single bam files for each of the 143 individuals using sinto (https://timoast.github.io/sinto/). Duplicated bam files for the same sample were merged together with samtools.103 For variant calling, the Sentieon germline variant calling pipeline104 was used, namely: