Skip to content
Snippets Groups Projects
Commit a8a06cd7 authored by DavyCats's avatar DavyCats
Browse files

add star task

parent 034e9ac2
No related branches found
No related tags found
1 merge request!2Star
.idea/
star.wdl 0 → 100644
task Star {
String? preCommand
File inputR1
File? inputR2
String genomeDir
String outFileNamePrefix
String? outSAMtype = "BAM SortedByCoordinate"
String? readFilesCommand = "zcat"
Int? runThreadN
String? outStd
String? twopassMode
String? outSAMattrRGline
Map[String, String] samOutputNames = {"BAM SortedByCoordinate": "sortedByCoord.out.bam"} #needs to be extended for all possible values
command {
set -e -o pipefail
mkdir -p ${sub(outFileNamePrefix, basename(outFileNamePrefix) + "$", "")}
${preCommand}
STAR \
--readFilesIn ${inputR1} ${inputR2} \
--outFileNamePrefix ${outFileNamePrefix} \
--genomeDir ${genomeDir} \
${"--readFilesCommand " + readFilesCommand} \
${"--outSAMtype " + outSAMtype} \
${"--runThreadN " + runThreadN} \
${"--outStd " + outStd} \
${"--twopassMode " + twopassMode} \
${"--outSAMattrRGline " + outSAMattrRGline}
}
output {
File bamFile = outFileNamePrefix + "Aligned." + samOutputNames["${outSAMtype}"]
}
}
\ No newline at end of file
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