Skip to content
Snippets Groups Projects
Commit 1dbfac51 authored by Ruben Vorderman's avatar Ruben Vorderman
Browse files

macs2 and htseq.wdl to 1.0

parent 0e6bde49
No related branches found
No related tags found
1 merge request!34Move all files to version 1.0
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 {
......
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 {
......
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