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

index solved

parent 9d7eaa0d
No related branches found
No related tags found
No related merge requests found
...@@ -27,7 +27,7 @@ task Bowtie { ...@@ -27,7 +27,7 @@ task Bowtie {
Array[File]+ readsUpstream Array[File]+ readsUpstream
Array[File]+? readsDownstream Array[File]+? readsDownstream
String outputPath = "mapped.bam" String outputPath = "mapped.bam"
BowtieIndex index Array[File]+ indexFiles
Int? seedmms Int? seedmms
Int? seedlen Int? seedlen
Int? k Int? k
...@@ -41,7 +41,7 @@ task Bowtie { ...@@ -41,7 +41,7 @@ task Bowtie {
# Image contains bowtie=1.2.2 and picard=2.9.2 # Image contains bowtie=1.2.2 and picard=2.9.2
String dockerImage = "quay.io/biocontainers/mulled-v2-bfe71839265127576d3cd749c056e7b168308d56:1d8bec77b352cdcf3e9ff3d20af238b33ed96eae" String dockerImage = "quay.io/biocontainers/mulled-v2-bfe71839265127576d3cd749c056e7b168308d56:1d8bec77b352cdcf3e9ff3d20af238b33ed96eae"
} }
String indexBasename = sub(index.fasta, "\.fa(sta)?$", "") String indexBasename = sub(indexFiles[0], "(\.rev)?\.[0-9]\.ebwt$", "")
# Assume fastq input with -q flag. # Assume fastq input with -q flag.
...@@ -55,7 +55,8 @@ task Bowtie { ...@@ -55,7 +55,8 @@ task Bowtie {
~{"--threads " + threads} \ ~{"--threads " + threads} \
~{true="--sam" false="" sam} \ ~{true="--sam" false="" sam} \
~{"--sam-RG " + samRG} \ ~{"--sam-RG " + samRG} \
-1 ~{sep="," readsUpstream} \ ~{indexBasename} \
~{true="-1" false="" defined(readsDownstream)} ~{sep="," readsUpstream} \
~{true="-2" false="" defined(readsDownstream)} ~{sep="," readsDownstream} \ ~{true="-2" false="" defined(readsDownstream)} ~{sep="," readsDownstream} \
| picard -Xmx~{picardMemory}G SortSam \ | picard -Xmx~{picardMemory}G SortSam \
INPUT=/dev/stdin \ INPUT=/dev/stdin \
......
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