How to Call Mutations From Your Own Whole-Genome Sequencing Data
A working pipeline that turns 30x whole-genome FASTQs into a filtered and benchmarked VCF, plus what changes when the variants are somatic.
8 posts
A working pipeline that turns 30x whole-genome FASTQs into a filtered and benchmarked VCF, plus what changes when the variants are somatic.
A working pipeline from FASTQ to annotated VCF for one human genome, with DeepVariant and GATK calling and benchmarking against GIAB truth sets.
A working GATK pipeline for calling SNPs and indels from RNA-seq reads, with filters suited to transcriptomic data and a check against genome calls.
A working end-to-end path from raw FASTQ to an annotated, QC'd variant set using GATK, DeepVariant, and VEP, with the flags, file sizes, runtimes, and failure modes that matter.
A step-by-step pipeline from raw FASTQ to a filtered, benchmarked VCF for a single human genome, with the tool choices, parameters, and quality checks we would use.
A working pipeline for calling somatic SNVs and indels with GATK Mutect2, from BAM to annotated VCF. Covers results at low variant allele fraction.
A practical comparison of DeepVariant, GATK HaplotypeCaller, bcftools mpileup, Strelka2, and Clair3 for calling SNPs from your own whole-genome data, with commands, benchmarks, and the failure modes that matter.
A working end-to-end guide to going from FASTQ to a filtered, benchmarked VCF on a personal 30x human genome, with the specific tools, flags, reference files, and failure modes we would use.