The sequenced reads in the FASTQ files were aligned to the human genome assembly (hg19) using the Burrows-Wheeler aligner (BWA). The initial alignment BAM files were subjected to sorting (SAMtools), removing duplicated reads (Picard), locally realigning reads around potential small insertion/deletion and recalibrating base quality score [Genome Analysis Toolkit (GATK)]. We used MuTect to generate high-confidence predictions on mutation calls. Variant Effector Predictor was used to annotate the called mutations.
Do you have any questions about this protocol?
Post your question to gather feedback from the community. We will also invite the authors of this article to respond.