Running realignment and recalibration, GATK

This function executes the docker container snv.1 where GATK software is used to do INDEL realignment and quality recalibration. This analysis is required only to run mutect1. The bwa index has to be prepared with bwaIndex

gatkDNA(
  group = c("sudo", "docker"),
  bam.folder = getwd(),
  scratch.folder = "/data/scratch",
  gatk.filename,
  genome.folder,
  threads = 1
)

Arguments

group,

a character string. Two options: "sudo" or "docker", depending to which group the user belongs

bam.folder,

a character string indicating where bam files generated with bwa.R are located. In this folder should be loacted also the GATK file GenomeAnalysisTK-X.X-0.tar.bz2.

scratch.folder,

a character string indicating the scratch folder where docker container will be mounted

gatk.filename,

a character string for GenomeAnalysisTK-X.X-0.tar.bz2.

genome.folder,

a character string indicating the folder where the indexed reference genome for bwa is located

threads,

a number indicating the number of cores to be used from the application

Value

three files: dedup_reads.bam, which is sorted and duplicates marked bam file, dedup_reads.bai, which is the index of the dedup_reads.bam, and dedup_reads.stats, which provides mapping statistics

Examples