This function executes the docker container snvPreprocessing. IMPORTANT: the working folder require the presence of GenomeAnalysisTK.tar.bz2 files

snvPreprocessing(
  group = c("sudo", "docker"),
  bam.folder = getwd(),
  scratch.folder = "/data/scratch",
  genome.folder,
  threads = 8,
  gatk.file
)

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 file generated by bwa is located

scratch.folder,

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

genome.folder,

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

threads,

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

gatk.file,

a character string for GenomeAnalysisTK-X.Y.tar.bz2, this file should be located in the bam.folder

Value

bam and bai after indel relignment and quality recalibration.

Examples

if (FALSE) { #downloading examples 1 million reads of mcf7 exome mixed with 1 million of mouse derived by human exome capturing system("wget http://130.192.119.59/public/hs1m_mm1m_R1.fastq.gz") system("wget http://130.192.119.59/public/hs1m_mm1m_R2.fastq.gz") #required for bwa 61Gb At the present time this is required to run mutect1 system("wget http://130.192.119.59/public/hg19_exome.tar.gz") #this genome index is prepared to be used for indel realignment and quality recalibration #running snvPreprocessing snvPreprocessing(group="docker",bam.folder=getwd(), scratch.folder="/data/scratch", genome.folder="/data/scratch/hg19_exome", threads=24, gatk.file="GenomeAnalysisTK-3.7.tar.bz2") }