Skip to content
Snippets Groups Projects
Commit 70ae489e authored by António Paulo's avatar António Paulo
Browse files

add GATK FilterMutectCalls as a task

parent 00aa5fd1
No related branches found
No related tags found
No related merge requests found
...@@ -318,7 +318,7 @@ task LearnReadOrientationModel { ...@@ -318,7 +318,7 @@ task LearnReadOrientationModel {
} }
output { output {
File artifactPriorTable = "artifact-priors.tar.gz" File artifactPriorsTable = "artifact-priors.tar.gz"
} }
runtime { runtime {
...@@ -420,6 +420,55 @@ task CalculateContamination { ...@@ -420,6 +420,55 @@ task CalculateContamination {
} }
} }
task FilterMutectCalls {
input {
File referenceFasta
File referenceFastaFai
File referenceFastaDict
File unfilteredVcf
File unfilteredVcfIndex
String outputVcf
File? contaminationTable
File? mafTumorSegments
File? artifactPriors
Int uniqueAltReadCount = 4
File mutect2Stats
String? extraArgs
Int memory = 4
Float memoryMultiplier = 1.5
String dockerImage = "quay.io/biocontainers/gatk4:4.1.2.0--1"
}
command {
set -e
gatk --java-options -Xmx~{memory}G \
FilterMutectCalls \
-R ~{referenceFasta} \
-V ~{unfilteredVcf} \
-O ${outputVcf} \
~{"--contamination-table " + contaminationTable} \
~{"--tumor-segmentation " + mafTumorSegments} \
~{"--ob-priors " + artifactPriors} \
~{"--unique-alt-read-count " + uniqueAltReadCount} \
~{"-stats " + mutect2Stats} \
--filtering-stats "filtering.stats" \
--showHidden \
~{extraArgs}
}
output {
File filteredVcf = outputVcf
File filteredVcfIndex = outputVcf + ".tbi"
File filteringStats = "filtering.stats"
}
runtime {
docker: dockerImage
memory: ceil(memory * memoryMultiplier)
}
}
task SplitNCigarReads { task SplitNCigarReads {
input { input {
File inputBam File inputBam
......
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