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 )
a character string. Two options:
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.
a character string indicating the scratch folder where docker container will be mounted
a character string for GenomeAnalysisTK-X.X-0.tar.bz2.
a character string indicating the folder where the indexed reference genome for bwa is located
a number indicating the number of cores to be used from the application
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