We aligned the Illumina short-read sequences to the GRCh37+decoy reference genome with BWA-mem56 using the default parameters. Putative PCR duplicates were flagged using SAMBLASTER57. The SAM outputs were converted to BAM format, and sorted by chromosomal coordinates using Sambamba58, and all BAM files for the same samples were merged.