Alignment and Calling

Overall description of VarSome's pipelines

VarSome Clinical uses Sentieon's implementation of bwa-mem for aligning reads to the genome and Sentieon's implementation of GATK for variant calling [1].

  • Sentieon: Sentieon (bwa-mem) aligner & Sentieon DNAscope caller
  • Tumor/somatic samples: Sentieon (bwa-mem) aligner & Sentieon TNhaplotyper2 caller (see below)

Default Workflow

  1. Adapter trimming is normally done before the FASTQ file is generated. The FASTQ files VarSome gets are therefore already trimmed.
  2. Read mapping
  3. Primer clipping
  4. Deduplication
  5. Indel realignment
  6. Base-quality recalibration
  7. Variant calling

Joint calling

Please see our article on joint calling.

Tumor Samples

 We use Sentieon's TNhaplotyper2 tool for both matched tumor-normal pairs and single tumor samples.

Related articles