Sequencing reads in FASTQ format were generated from Illumina image data using bcl2fastq program (Illumina). Following the OQFE (original quality functional equivalent) protocol74 (link), sequence reads were mapped to GRCh38 references using BWA MEM75 (link) in an alt-aware manner, read duplicates were marked, and additional per-read tags were added. Single nucleotide variations (SNV) and short insertion and deletions (indels) were identified using a Parabricks accelerated version of DeepVariant v0.10 with a custom WES model and reported in per-sample genome VCF (gVCF)76 (link). These gVCFs were aggregated with GLnexus v1.4.377 (link) into joint-genotyped multi-sample project-level VCF (pVCF), which was converted to bed/bim/fam format using PLINK 1.978 (link) for downstream analyses.