Merge pull request #217 from biowdl/BIOWDL-457

Add necessary options for RNA genotyping to haplotypecaller.

CHANGELOG.md

@@ -11,6 +11,9 @@ that users understand how the changes affect the new version.
 
 version 3.2.0-develop
 ---------------------------
++ GATK.HaplotypeCaller: Add `--dont-use-soft-clipped-bases` and 
+  `--standard-min-confidence-threshold-for-calling` options. These are 
+  required for RNA seq variant calling according to GATK best practices.
 + Samtools: Fix quotations in sort command.
 + Samtools SortByName is now called Sort.
 + Generalize sort task to now also sort by position, instead of just read name.

gatk.wdl

@@ -910,6 +910,8 @@ task HaplotypeCaller {
         String? outputMode
         Boolean gvcf = false
         String emitRefConfidence = if gvcf then "GVCF" else "NONE"
+        Boolean dontUseSoftClippedBases = false
+        Float? standardMinConfidenceThresholdForCalling
 
         String memory = "12G"
         String javaXmx = "4G"
@@ -931,7 +933,9 @@ task HaplotypeCaller {
         ~{"--pedigree " + pedigree} \
         ~{"--contamination-fraction-per-sample-file " + contamination} \
         ~{"--output-mode " + outputMode} \
-        --emit-ref-confidence ~{emitRefConfidence}
+        --emit-ref-confidence ~{emitRefConfidence} \
+        ~{true="--dont-use-soft-clipped-bases" false="" dontUseSoftClippedBases} \
+        ~{"--standard-min-confidence-threshold-for-calling " + standardMinConfidenceThresholdForCalling}
     }
 
     output {
@@ -962,6 +966,8 @@ task HaplotypeCaller {
                      category: "advanced"}
         emitRefConfidence: {description: "Whether to include reference calls. Three modes: 'NONE', 'BP_RESOLUTION' and 'GVCF'",
                             category: "advanced"}
+        dontUseSoftClippedBases: {description: "Do not use soft-clipped bases. Should be 'true' for RNA variant calling.", category: "common"}
+        standardMinConfidenceThresholdForCalling: {description: "Confidence threshold used for calling variants.", category: "advanced"}
         dbsnpVCF: {description: "A dbSNP VCF.", category: "common"}
         dbsnpVCFIndex: {description: "The index for the dbSNP VCF.", category: "common"}
         pedigree: {description: "Pedigree file for determining the population \"founders\"", category: "common"}