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).