diff --git a/macs2.wdl b/macs2.wdl new file mode 100644 index 0000000000000000000000000000000000000000..999f3f0fdb570289d0614d605741218aebc91984 --- /dev/null +++ b/macs2.wdl @@ -0,0 +1,24 @@ +task peakCalling { + String? preCommand + File bamFile + String outDir + String sampleName + Int? threads + Int? memory + + + command { + set -e -o pipefail + ${preCommand} + macs2 callpeaks --treatment ${bamFile} --outdir ${outDir} --name ${sampleName} + } + + output { + File peakFile = outDir + "/" + sampleName + "/macs2/" + sampleName + "_peaks.narrowPeak" + } + + runtime { + cpu: select_first([threads,1]) + memory: select_first([memory,8]) + } +} \ No newline at end of file