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

solved indexing issue

parent 0b157958
No related branches found
No related tags found
2 merge requests!9Changes for virus assembly pipeline,!10Extra tasks required for assembly.
...@@ -70,7 +70,6 @@ task mem { ...@@ -70,7 +70,6 @@ task mem {
${true="-p " false="" smartPairing} \ ${true="-p " false="" smartPairing} \
${"-r " + readGroupHeaderLine} \ ${"-r " + readGroupHeaderLine} \
${"-H " + H } \ ${"-H " + H } \
${"-o " + outputFile } \
${true="-j" false="" j } \ ${true="-j" false="" j } \
${true="-5" false="" five } \ ${true="-5" false="" five } \
${true="-q" false="" q } \ ${true="-q" false="" q } \
...@@ -86,10 +85,12 @@ task mem { ...@@ -86,10 +85,12 @@ task mem {
${indexBase} ${read1} ${read2} \ ${indexBase} ${read1} ${read2} \
| picard SortSam CREATE_INDEX=TRUE TMP_DIR=null \ | picard SortSam CREATE_INDEX=TRUE TMP_DIR=null \
INPUT=/dev/stdin SORT_ORDER=coordinate OUTPUT=${outputFile} INPUT=/dev/stdin SORT_ORDER=coordinate OUTPUT=${outputFile}
ln -s $(basename ${outputFile} | sed 's/.bam$/.bai/') ${outputFile}.bai
} }
output { output {
File bamFile = select_first([outputFile]) File bamFile = select_first([outputFile])
File bamIndex = sub(bamFile, ".bam$", ".bai") File bamIndexPicard = sub(bamFile, ".bam$", ".bai")
File bamIndexSamtools = select_first([outputFile]) + ".bai"
} }
runtime { runtime {
cpu: select_first([threads]) 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