Sex determination
We calculated the depth of coverage of each sample in windows of 10kb using Mosdepth (47). Next, we normalized each window coverage with the mean coverage of each individual to exclude the influence of sequence effort and plotted the density of normalized mean coverage in the autosomes and the sex chromosome (chrZ) to genetically identify the sex of each sample (half of the average coverage is expected in the chrZ for females).