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

output as sorted and indexed bam

parent 58b232cb
No related branches found
No related tags found
2 merge requests!9Changes for virus assembly pipeline,!10Extra tasks required for assembly.
......@@ -2,7 +2,7 @@ task mem {
File read1
File referenceFile
File? read2
String? outputFile
String? outputFile = "aligned.bam"
String? preCommand
Int? threads = 1
String? memory = "4G"
......@@ -78,10 +78,12 @@ task mem {
${true="-Y" false="" Y } \
${true="-M" false="" M } \
${"-I " + I } \
${referenceFile} ${read1} ${read2}
${referenceFile} ${read1} ${read2} \
| picard SortSam CREATE_INDEX=TRUE TMP_DIR=null \
INPUT=/dev/stdin OUTPUT=${outputFile}
}
output {
File samFile = if defined(outputFile) then select_first([outputFile]) else stdout()
File samFile = select_first([outputFile])
}
runtime {
cpu: select_first([threads])
......
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