Hybrid selection was performed using custom targeted Nimblegen Roche probes. Libraries were sequenced on HiSeq2000 sequencers (Illumina Inc, San Diego CA). Reads were aligned with BWA (Li et al 2009) and SNVs were called with GATK (McKenna et al 2010).