This is an HPC workflow to call SNPs and INDELs in the mitochondrial genome from matched tumor-normal whole genome sequencing (WGS) data from 25 childhood acute lymphoblastic leukemia cases. This workflow is largely based on GATK's best practice WDL. Check out GATK's doc for more details.
Step 2: Align unmapped bam, mark duplicates, collect coverage information, and call SNPs and INDELs from normal bams
- 2.1.AlignToMt.sh
- 2.2.AlignToShiftedMt.sh
- 2.3.1.CollectWgsMetrics.sh
- 2.3.2.CollectWgsMetrics.R
- 2.3.3.RunCollectWgsMetricsR.sh
- 2.4.CallMt.Germline.sh
- 2.5.CallShiftedMt.Germline.sh
- 2.6.LiftoverAndCombineVcfs.sh
- 2.7.MergeStats.sh
- 2.8.InitialFilter.sh
- 2.9.SplitMultiAllelicsAndRemoveNonPassSites.sh