In the sequencing stage, reads were aligned to human reference genome 19 using BWA software (Li and Durbin, 2009). The per-base quality scores were recalibrated, and reads were locally realigned using the Genome Analysis Toolkit (DePristo et al., 2011). Syzygy software (Rivas et al., 2011) was used to call variants from reads generated from pooled DNA samples.
