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

debugged spades

parent 47f2f310
No related branches found
No related tags found
2 merge requests!9Changes for virus assembly pipeline,!10Extra tasks required for assembly.
task sample {
File sequenceFile
String? outFilePath = "subsampledReads"
String? outFilePath = "subsampledReads.fq.gz"
String? preCommand
Int? seed
Boolean? twoPassMode
Float? fraction
Int? number
Boolean? zip
Boolean? zip = true
command {
set -e -o pipefail
......@@ -20,6 +20,6 @@ task sample {
${"> " + outFilePath}
}
output {
File subsampledReads=select_first([outFilePath])
File subsampledReads= select_first([outFilePath])
}
}
\ No newline at end of file
......@@ -39,7 +39,7 @@ task spades {
${true="--rna" false="" rna} \
${true="--plasmid" false="" plasmid} \
${true="--iontorrent" false="" ionTorrent} \
${"--12 " + interlacedReads }
${"--12 " + interlacedReads } \
${"-1 " + read1 } \
${"-2 " + read2 } \
${"-s " + singleRead } \
......@@ -65,12 +65,12 @@ task spades {
Array[File] correctedReads = glob(outputDir + "/corrected/*.fastq*")
File scaffolds = outputDir + "/scaffolds.fasta"
File contigs = outputDir + "/contigs.fasta"
File assemblyGraphGfa = outputDir + "/assembly_graph.gfa"
File assemblyGraphFastg = outputDir + "/assembly_graph.fastq"
File assemblyGraphWithScaffoldsGfa = outputDir + "/assembly_graph_with_scaffolds.gfa"
File assemblyGraphFastg = outputDir + "/assembly_graph.fastg"
File contigsPaths = outputDir + "/contigs.paths"
File scaffoldsPaths = outputDir + "/scaffolds.paths"
File params = outputDir + "/params.txt"
File log = outputDir + "/spades/log"
File log = outputDir + "/spades.log"
}
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