Fastq files from the sequencer were redirected to a custom pipeline for HaloPlex Target Enrichment System on the DNA nexus platform and/or to Agilent Surecall software.
Briefly, reads were aligned to the human reference genome (GRCh37/hg19) (http://hgdownload.cse.ucsc.edu/) with Burrows-Wheeler Aligner (BWA) [28 (link)] and variants were called using at least 2 of the 3 following variant callers: Genome Analysis Toolkit (GATK) [29 (link)–31 ], Freebayes [32 ] (both within the DNA nexus platform), and Base Alignment Quality (BAQ) Single Nucleotide Polymorphism (SNP) caller (within SureCall tool).
Variants passing quality filters were annotated separately against NCBI RefGene (http://www.ncbi.nlm.nih.gov) and ENSEMBL Variant Effect Predictor ver.72 (http://www.ensembl.org/info/docs/tools/vep).
Free full text: Click here