diff --git a/htseq.wdl b/htseq.wdl index b634bf5ee0f08128729723eaf71e77536bf401f6..87c4c676fe5ea2485837449b0d802ea52e927fe4 100644 --- a/htseq.wdl +++ b/htseq.wdl @@ -1,25 +1,29 @@ +version 1.0 + task HTSeqCount { - String? preCommand - Array[File] alignmentFiles - File gtfFile - String outputTable - String? format - String? order - String? stranded + input { + String? preCommand + Array[File] alignmentFiles + File gtfFile + String outputTable + String? format + String? order + String? stranded - Int? memory + Int? memory + } command { set -e -o pipefail - mkdir -p ${sub(outputTable, basename(outputTable), "")} - ${preCommand} + mkdir -p ~{sub(outputTable, basename(outputTable), "")} + ~{preCommand} htseq-count \ - -f ${default="bam" format} \ - -r ${default="pos" order} \ - -s ${default="no" stranded} \ - ${sep=" " alignmentFiles} \ - ${gtfFile} \ - > ${outputTable} + -f ~{default="bam" format} \ + -r ~{default="pos" order} \ + -s ~{default="no" stranded} \ + ~{sep=" " alignmentFiles} \ + ~{gtfFile} \ + > ~{outputTable} } output { diff --git a/macs2.wdl b/macs2.wdl index 9880c944a6b0942effc88472f64487b165d15d39..026b418eb1494117afd5149ba5a750e4397d6925 100644 --- a/macs2.wdl +++ b/macs2.wdl @@ -1,20 +1,23 @@ -task PeakCalling { - String? preCommand - Array[File] bamFiles - String outDir - String sampleName - Int? threads - Int? memory - Boolean? nomodel +version 1.0 +task PeakCalling { + input { + String? preCommand + Array[File] bamFiles + String outDir + String sampleName + Int? threads + Int? memory + Boolean? nomodel + } command { set -e -o pipefail - ${preCommand} + ~{preCommand} macs2 callpeak \ - --treatment ${sep = ' ' bamFiles} \ - --outdir ${outDir} \ - --name ${sampleName} \ - ${default=false true='--nomodel' false='' nomodel} + --treatment ~{sep = ' ' bamFiles} \ + --outdir ~{outDir} \ + --name ~{sampleName} \ + ~{default=false true='--nomodel' false='' nomodel} } output {