Merge pull request #64 from biowdl/BIOPET-1080

Added CaseControl tool

biopet/sampleconfig.wdl

 version 1.0
 
+import "../common.wdl" as common
+
 task SampleConfig {
     input {
         File? toolJar
@@ -77,3 +79,41 @@ task SampleConfigCromwellArrays {
         memory: ceil(memory * memoryMultiplier)
     }
 }
+
+task CaseControl {
+    input {
+        File? toolJar
+        String? preCommand
+        Array[File]+ inputFiles
+        Array[File]+ sampleConfigs
+        String outputPath
+        String controlTag = "control"
+
+        Int memory = 4
+        Float memoryMultiplier = 1.5
+    }
+
+    String toolCommand = if defined(toolJar)
+        then "java -Xmx" + memory + "G -jar " + toolJar
+        else "biopet-sampleconfig -Xmx" + memory + "G"
+
+    command {
+        set -e -o pipefail
+        ~{preCommand}
+        mkdir -p $(dirname ~{outputPath})
+        ~{toolCommand} CromwellArrays \
+        -i ~{sep="-i " inputFiles} \
+        -s ~{sep="-s " sampleConfigs} \
+        ~{"-o " + outputPath} \
+        ~{"--controlTag " + controlTag}
+    }
+
+    output {
+        File outputFile = outputPath
+        Array[CaseControl] caseControls = read_json(outputFile)
+    }
+
+    runtime {
+        memory: ceil(memory * memoryMultiplier)
+    }
+}

common.wdl

@@ -163,3 +163,10 @@ struct FastqPair {
     File? R2
     File? R2_md5
 }
+
+struct CaseControl {
+    String inputName
+    IndexedBamFile inputBam
+    String controlName
+    IndexedBamFile controlBam
+}