wes
Genomic DNA samples were collected and sequence libraries were constructed using the Agilent Sure Select Human Whole Exome V2 Kit (Agilent Technologies, Santa Clara, CA). Prepared libraries were sequenced using the HiSeq2500 System (Illumina, San Diego, CA). Reads obtained from the BWA software package (v. 0.7.15) was mapped with the human reference genome (GRCh37/hg19). Variant calling and variant annotation were performed using the Genome Analysis Toolkit (GATK)