Skip to content
Snippets Groups Projects
Commit 8b0d6891 authored by Cats's avatar Cats
Browse files

unpaired vardict fix, split strelka/manta configuration and run to fix retries

parent 1bc3f6cc
No related branches found
No related tags found
1 merge request!45Changes for somatic variantcalling
version 1.0
task Somatic {
task ConfigureSomatic {
input {
File tumorBam
File tumorIndex
......@@ -14,9 +14,6 @@ task Somatic {
Boolean exome = false
String? preCommand
String? installDir
Int cores = 1
Int memory = 4
}
String toolCommand = if defined(installDir)
......@@ -33,7 +30,22 @@ task Somatic {
~{"--callRegions " + callRegions} \
--runDir ~{runDir} \
~{true="--exome" false="" exome}
}
output {
String runDirectory = runDir
}
}
task RunSomatic {
input {
String runDir
Int cores = 1
Int memory = 4
Boolean paired = true
}
command {
~{runDir}/runWorkflow.py \
-m local \
-j ~{cores} \
......@@ -46,10 +58,10 @@ task Somatic {
"/results/variants/candidateSmallIndels.vcf.gz.tbi"
File candidateSV = runDir + "/results/variants/candidateSV.vcf.gz"
File candidateSVindex = runDir + "/results/variants/candidateSV.vcf.gz.tbi"
File tumorSV = if defined(normalBam)
File tumorSV = if paired
then runDir + "/results/variants/somaticSV.vcf.gz"
else runDir + "/results/variants/tumorSV.vcf.gz"
File tumorSVindex = if defined(normalBam)
File tumorSVindex = if paired
then runDir + "/results/variants/somaticSV.vcf.gz.tbi"
else runDir + "/results/variants/tumorSV.vcf.gz.tbi"
File? diploidSV = "/results/variants/diploidSV.vcf.gz"
......
version 1.0
task Germline {
task ConfigureGermline {
input {
String? preCommand
String? installDir
......@@ -13,9 +13,6 @@ task Germline {
File? callRegionsIndex
Boolean exome = false
Boolean rna = false
Int cores = 1
Int memory = 4
}
String toolCommand = if defined(installDir)
......@@ -32,26 +29,14 @@ task Germline {
~{"--callRegions " + callRegions} \
~{true="--exome" false="" exome} \
~{true="--rna" false="" rna}
~{runDir}/runWorkflow.py \
-m local \
-j ~{cores} \
-g ~{memory}
}
output {
File variants = runDir + "/results/variants/variants.vcf.gz"
File variantsIndex = runDir + "/results/variants/variants.vcf.gz.tbi"
}
runtime {
cpu: cores
memory: memory
String runDirectory = runDir
}
}
task Somatic {
task ConfigureSomatic {
input {
String? preCommand
String? installDir
......@@ -67,9 +52,6 @@ task Somatic {
File? indelCandidates
File? indelCandidatesIndex
Boolean exome = false
Int cores = 1
Int memory = 4
}
String toolCommand = if defined(installDir)
......@@ -87,7 +69,22 @@ task Somatic {
~{"--callRegions " + callRegions} \
~{"--indelCandidates " + indelCandidates} \
~{true="--exome" false="" exome} \
}
output {
String runDirectory = runDir
}
}
task Run {
input {
String runDir
Int cores = 1
Int memory = 4
Boolean somatic = true
}
command {
~{runDir}/runWorkflow.py \
-m local \
-j ~{cores} \
......@@ -95,10 +92,14 @@ task Somatic {
}
output {
File indelsVcf = runDir + "/results/variants/somatic.indels.vcf.gz"
File indelsIndex = runDir + "/results/variants/somatic.indels.vcf.gz.tbi"
File snvVcf = runDir + "/results/variants/somatic.snvs.vcf.gz"
File snvIndex = runDir + "/results/variants/somatic.snvs.vcf.gz.tbi"
File? indelsVcf = runDir + "/results/variants/somatic.indels.vcf.gz"
File? indelsIndex = runDir + "/results/variants/somatic.indels.vcf.gz.tbi"
File variants = if somatic
then runDir + "/results/variants/somatic.snvs.vcf.gz"
else runDir + "/results/variants/variants.vcf.gz"
File variantsIndex = if somatic
then runDir + "/results/variants/somatic.snvs.vcf.gz.tbi"
else runDir + "/results/variants/variants.vcf.gz.tbi"
}
runtime {
......
......@@ -39,7 +39,7 @@ task VarDict {
-G ~{refFasta} \
-N ~{tumorSampleName} \
-b "~{tumorBam}~{"|" + normalBam}" \
~{true="-z" false="" defined(normalBam)} \
~{true="" false="-z" defined(normalBam)} \
-c ~{chromosomeColumn} \
-S ~{startColumn} \
-E ~{endColumn} \
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment