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

spades to 1.0

parent 8c40b9e5
No related branches found
No related tags found
1 merge request!34Move all files to version 1.0
version 1.0
task spades {
String outputDir
String? preCommand
File read1
File? read2
File? interlacedReads
File? sangerReads
File? pacbioReads
File? nanoporeReads
File? tslrContigs
File? trustedContigs
File? untrustedContigs
Boolean? singleCell
Boolean? metagenomic
Boolean? rna
Boolean? plasmid
Boolean? ionTorrent
Boolean? onlyErrorCorrection
Boolean? onlyAssembler
Boolean? careful
Boolean? disableGzipOutput
Boolean? disableRepeatResolution
File? dataset
Int? threads
Float? memoryGb
File? tmpDir
String? k
Float? covCutoff
Int? phredOffset
input {
String outputDir
String? preCommand
File read1
File? read2
File? interlacedReads
File? sangerReads
File? pacbioReads
File? nanoporeReads
File? tslrContigs
File? trustedContigs
File? untrustedContigs
Boolean? singleCell
Boolean? metagenomic
Boolean? rna
Boolean? plasmid
Boolean? ionTorrent
Boolean? onlyErrorCorrection
Boolean? onlyAssembler
Boolean? careful
Boolean? disableGzipOutput
Boolean? disableRepeatResolution
File? dataset
Int? threads
Float? memoryGb
File? tmpDir
String? k
Float? covCutoff
Int? phredOffset
}
Int finalThreads = select_first([threads,1])
Float totalMemory = select_first([memoryGb, finalThreads * 16.0])
Int finalMemory = ceil(totalMemory)
......@@ -34,34 +38,34 @@ task spades {
command {
set -e -o pipefail
${preCommand}
~{preCommand}
spades.py \
${"-o " + outputDir} \
${true="--sc" false="" singleCell} \
${true="--meta" false="" metagenomic} \
${true="--rna" false="" rna} \
${true="--plasmid" false="" plasmid} \
${true="--iontorrent" false="" ionTorrent} \
${"--12 " + interlacedReads } \
${true="-1" false="-s" defined(read2)} ${read1} \
${"-2 " + read2 } \
${"--sanger " + sangerReads } \
${"--pacbio " + pacbioReads } \
${"--nanopore " + nanoporeReads } \
${"--tslr " + tslrContigs } \
${"--trusted-contigs " + trustedContigs } \
${"--untrusted-contigs " + untrustedContigs } \
${true="--only-error-correction" false="" onlyErrorCorrection } \
${true="--only-assembler" false="" onlyAssembler } \
${true="--careful" false="" careful } \
${true="--disable-gzip-output" false="" disableGzipOutput} \
${true="--disable-rr" false="" disableRepeatResolution } \
${"--dataset " + dataset } \
${"--threads " + finalThreads} \
${"--memory " + finalMemory } \
${"-k " + k } \
${"--cov-cutoff " + covCutoff } \
${"--phred-offset " + phredOffset }
~{"-o " + outputDir} \
~{true="--sc" false="" singleCell} \
~{true="--meta" false="" metagenomic} \
~{true="--rna" false="" rna} \
~{true="--plasmid" false="" plasmid} \
~{true="--iontorrent" false="" ionTorrent} \
~{"--12 " + interlacedReads } \
~{true="-1" false="-s" defined(read2)} ~{read1} \
~{"-2 " + read2 } \
~{"--sanger " + sangerReads } \
~{"--pacbio " + pacbioReads } \
~{"--nanopore " + nanoporeReads } \
~{"--tslr " + tslrContigs } \
~{"--trusted-contigs " + trustedContigs } \
~{"--untrusted-contigs " + untrustedContigs } \
~{true="--only-error-correction" false="" onlyErrorCorrection } \
~{true="--only-assembler" false="" onlyAssembler } \
~{true="--careful" false="" careful } \
~{true="--disable-gzip-output" false="" disableGzipOutput} \
~{true="--disable-rr" false="" disableRepeatResolution } \
~{"--dataset " + dataset } \
~{"--threads " + finalThreads} \
~{"--memory " + finalMemory } \
~{"-k " + k } \
~{"--cov-cutoff " + covCutoff } \
~{"--phred-offset " + phredOffset }
}
output {
Array[File] correctedReads = glob(outputDir + "/corrected/*.fastq*")
......
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