Skip to content
Snippets Groups Projects
Commit c0bbdd8c authored by pjvan_thof's avatar pjvan_thof
Browse files

Added gatk as a non jar command

parent ad685f61
No related branches found
No related tags found
No related merge requests found
# Apply Base Quality Score Recalibration (BQSR) model # Apply Base Quality Score Recalibration (BQSR) model
task ApplyBQSR { task ApplyBQSR {
String? preCommand String? preCommand
File gatkJar File? gatkJar
File inputBam File inputBam
File inputBamIndex File inputBamIndex
String outputBamPath String outputBamPath
...@@ -16,10 +16,15 @@ task ApplyBQSR { ...@@ -16,10 +16,15 @@ task ApplyBQSR {
Float? memoryMultiplier Float? memoryMultiplier
Int mem = ceil(select_first([memory, 4.0])) Int mem = ceil(select_first([memory, 4.0]))
String toolCommand = if defined(gatkJar)
then "java -Xmx" + mem + "G -jar " + gatkJar
else "gatk -Xmx" + mem + "G"
command { command {
set -e -o pipefail set -e -o pipefail
${preCommand} ${preCommand}
java ${"-Dsamjdk.compression_level=" + compressionLevel} \ ${toolCommand} \
-Xms${mem}G -jar ${gatkJar} \ -Xms${mem}G -jar ${gatkJar} \
ApplyBQSR \ ApplyBQSR \
--create-output-bam-md5 \ --create-output-bam-md5 \
...@@ -46,7 +51,7 @@ task ApplyBQSR { ...@@ -46,7 +51,7 @@ task ApplyBQSR {
# Generate Base Quality Score Recalibration (BQSR) model # Generate Base Quality Score Recalibration (BQSR) model
task BaseRecalibrator { task BaseRecalibrator {
String? preCommand String? preCommand
File gatkJar File? gatkJar
File inputBam File inputBam
File inputBamIndex File inputBamIndex
String recalibrationReportPath String recalibrationReportPath
...@@ -61,10 +66,15 @@ task BaseRecalibrator { ...@@ -61,10 +66,15 @@ task BaseRecalibrator {
Float? memoryMultiplier Float? memoryMultiplier
Int mem = ceil(select_first([memory, 4.0])) Int mem = ceil(select_first([memory, 4.0]))
String toolCommand = if defined(gatkJar)
then "java -Xmx" + mem + "G -jar " + gatkJar
else "gatk -Xmx" + mem + "G"
command { command {
set -e -o pipefail set -e -o pipefail
${preCommand} ${preCommand}
java -Xms${mem}G -jar ${gatkJar} \ ${toolCommand} \
BaseRecalibrator \ BaseRecalibrator \
-R ${refFasta} \ -R ${refFasta} \
-I ${inputBam} \ -I ${inputBam} \
...@@ -91,7 +101,7 @@ task CombineGVCFs { ...@@ -91,7 +101,7 @@ task CombineGVCFs {
String outputPath String outputPath
String gatkJar String? gatkJar
File refFasta File refFasta
File refFastaIndex File refFastaIndex
...@@ -102,12 +112,17 @@ task CombineGVCFs { ...@@ -102,12 +112,17 @@ task CombineGVCFs {
Float? memoryMultiplier Float? memoryMultiplier
Int mem = ceil(select_first([memory, 4.0])) Int mem = ceil(select_first([memory, 4.0]))
String toolCommand = if defined(gatkJar)
then "java -Xmx" + mem + "G -jar " + gatkJar
else "gatk -Xmx" + mem + "G"
command { command {
set -e -o pipefail set -e -o pipefail
${preCommand} ${preCommand}
if [ ${length(gvcfFiles)} -gt 1 ]; then if [ ${length(gvcfFiles)} -gt 1 ]; then
java ${"-Dsamjdk.compression_level=" + compressionLevel} \ ${toolCommand} \
-Xmx${mem}G -jar ${gatkJar} \ -Xmx${mem}G -jar ${gatkJar} \
CombineGVCFs \ CombineGVCFs \
-R ${refFasta} \ -R ${refFasta} \
...@@ -133,7 +148,7 @@ task CombineGVCFs { ...@@ -133,7 +148,7 @@ task CombineGVCFs {
# Combine multiple recalibration tables from scattered BaseRecalibrator runs # Combine multiple recalibration tables from scattered BaseRecalibrator runs
task GatherBqsrReports { task GatherBqsrReports {
String? preCommand String? preCommand
String gatkJar String? gatkJar
Array[File] inputBQSRreports Array[File] inputBQSRreports
String outputReportPath String outputReportPath
...@@ -141,10 +156,15 @@ task GatherBqsrReports { ...@@ -141,10 +156,15 @@ task GatherBqsrReports {
Float? memoryMultiplier Float? memoryMultiplier
Int mem = ceil(select_first([memory, 4.0])) Int mem = ceil(select_first([memory, 4.0]))
String toolCommand = if defined(gatkJar)
then "java -Xmx" + mem + "G -jar " + gatkJar
else "gatk -Xmx" + mem + "G"
command { command {
set -e -o pipefail set -e -o pipefail
${preCommand} ${preCommand}
java -Xms${mem}G -jar ${gatkJar} \ ${toolCommand} \
GatherBQSRReports \ GatherBQSRReports \
-I ${sep=' -I ' inputBQSRreports} \ -I ${sep=' -I ' inputBQSRreports} \
-O ${outputReportPath} -O ${outputReportPath}
...@@ -167,7 +187,7 @@ task GenotypeGVCFs { ...@@ -167,7 +187,7 @@ task GenotypeGVCFs {
String outputPath String outputPath
String gatkJar String? gatkJar
File refFasta File refFasta
File refFastaIndex File refFastaIndex
...@@ -181,11 +201,16 @@ task GenotypeGVCFs { ...@@ -181,11 +201,16 @@ task GenotypeGVCFs {
Float? memoryMultiplier Float? memoryMultiplier
Int mem = ceil(select_first([memory, 4.0])) Int mem = ceil(select_first([memory, 4.0]))
String toolCommand = if defined(gatkJar)
then "java -Xmx" + mem + "G -jar " + gatkJar
else "gatk -Xmx" + mem + "G"
command { command {
set -e -o pipefail set -e -o pipefail
${preCommand} ${preCommand}
java ${"-Dsamjdk.compression_level=" + compressionLevel} \ ${toolCommand} \
-Xmx${mem}G -jar ${gatkJar} \ -Xmx${mem}G -jar ${gatkJar} \
GenotypeGVCFs \ GenotypeGVCFs \
-R ${refFasta} \ -R ${refFasta} \
...@@ -220,16 +245,21 @@ task HaplotypeCallerGvcf { ...@@ -220,16 +245,21 @@ task HaplotypeCallerGvcf {
File refFastaIndex File refFastaIndex
Float? contamination Float? contamination
Int? compressionLevel Int? compressionLevel
String gatkJar String? gatkJar
Float? memory Float? memory
Float? memoryMultiplier Float? memoryMultiplier
Int mem = ceil(select_first([memory, 4.0])) Int mem = ceil(select_first([memory, 4.0]))
String toolCommand = if defined(gatkJar)
then "java -Xmx" + mem + "G -jar " + gatkJar
else "gatk -Xmx" + mem + "G"
command { command {
set -e -o pipefail set -e -o pipefail
${preCommand} ${preCommand}
java ${"-Dsamjdk.compression_level=" + compressionLevel} \ ${toolCommand} \
-Xmx${mem}G -jar ${gatkJar} \ -Xmx${mem}G -jar ${gatkJar} \
HaplotypeCaller \ HaplotypeCaller \
-R ${refFasta} \ -R ${refFasta} \
......
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