diff --git a/spades.wdl b/spades.wdl index 960446b966faf532078ee17298bbd910b62dc284..f1c94e8005d914414370016b885cc4664b10a924 100644 --- a/spades.wdl +++ b/spades.wdl @@ -24,14 +24,14 @@ task spades { Boolean? disableGzipOutput Boolean? disableRepeatResolution File? dataset - Int? threads + Int? threads = 1 Float? memoryGb File? tmpDir String? k Float? covCutoff Int? phredOffset } - Int finalThreads = select_first([threads,1]) + Int finalThreads = select_first([threads]) Float totalMemory = select_first([memoryGb, finalThreads * 16.0]) Int finalMemory = ceil(totalMemory) Int clusterMemory = ceil(totalMemory / finalThreads) diff --git a/unicycler.wdl b/unicycler.wdl index f7aa36a48067bcd2f880a843d1a4aa401c7f5078..9b442afd36fabe7ff422736cd69e118382bbb521 100644 --- a/unicycler.wdl +++ b/unicycler.wdl @@ -1,94 +1,96 @@ +version 1.0 + task unicycler { - String? preCommand - File? short1 - File? short2 - File? unpaired - File? long - String out - Int? verbosity - Int? minFastaLength - Int? keep - Boolean? vcf - Int? threads - Int? memory - Int finalThreads = select_first(threads, 1) - Int finalMemory = select_first(memory, 4) - String? mode - Float? minBridgeQual - Int? linearSeqs - File? spadesPath - Boolean? noCorrect - Float? minKmerFrac - Float? maxKmerFrac - Int? kmerCount - Float? depthFilter - Boolean? noMiniasm - File? raconPath - File? existingLongReadAssembly - Boolean? noRotate - File? startGenes - Float? startGeneId - Float? startGeneCov - String? makeblastdbPath - File? tblastnPath - Boolean? noPilon - File? bowtie2Path - File? bowtie2buildPath - File? samtoolsPath - File? pilonPath - File? javaPath - Int? minPolishSize - File? bcftoolsPath - Int? minComponentSize - Int? minDeadEndSize - File? contamination - String? scores - String? lowScore + input { + String? preCommand + File? short1 + File? short2 + File? unpaired + File? long + String out + Int? verbosity + Int? minFastaLength + Int? keep + Boolean? vcf + Int? threads = 1 + Int? memory = 4 + String? mode + Float? minBridgeQual + Int? linearSeqs + File? spadesPath + Boolean? noCorrect + Float? minKmerFrac + Float? maxKmerFrac + Int? kmerCount + Float? depthFilter + Boolean? noMiniasm + File? raconPath + File? existingLongReadAssembly + Boolean? noRotate + File? startGenes + Float? startGeneId + Float? startGeneCov + String? makeblastdbPath + File? tblastnPath + Boolean? noPilon + File? bowtie2Path + File? bowtie2buildPath + File? samtoolsPath + File? pilonPath + File? javaPath + Int? minPolishSize + File? bcftoolsPath + Int? minComponentSize + Int? minDeadEndSize + File? contamination + String? scores + String? lowScore + } command { set -e -o pipefail - mkdir -p ${out} - ${preCommand} + mkdir -p ~{out} + ~{preCommand} unicycler \ - ${"--short1 " + short1} \ - ${"--short2 " + short2} \ - ${"--unpaired " + unpaired} \ - ${"--long " + long} \ - --out ${out} \ - ${"--min_fasta_length " + minFastaLength} \ - ${"--keep " + keep } \ - ${true="--vcf" false="" vcf } \ - ${"--threads " + finalThreads } \ - ${"--mode " + mode } \ - ${"--min_bridge_qual " + minBridgeQual } \ - ${"--linear_seqs " + linearSeqs } \ - ${"--spades_path " + spadesPath } \ - ${true="--no_correct" false="" noCorrect } \ - ${"--min_kmer_frac " + minKmerFrac } \ - ${"--max_kmer_frac " + maxKmerFrac } \ - ${"--kmer_count " + kmerCount } \ - ${"--depth_filter " + depthFilter } \ - ${true="--no_miniasm" false="" noMiniasm } \ - ${"--racon_path " + raconPath } \ - ${"--existing_long_read_assembly " + existingLongReadAssembly } \ - ${true="--no_rotate" false="" noRotate } \ - ${"--start_genes " + startGenes } \ - ${"--start_gene_id " + startGeneId } \ - ${"--start_gene_cov " + startGeneCov } \ - ${"--makeblastdb_path " + makeblastdbPath } \ - ${"--tblastn_path " + tblastnPath } \ - ${true="--no_pilon" false="" noPilon } \ - ${"--bowtie2_path " + bowtie2Path } \ - ${"--bowtie2_build_path " + bowtie2buildPath } \ - ${"--samtools_path " + samtoolsPath } \ - ${"--pilon_path " + pilonPath } \ - ${"--java_path " + javaPath } \ - ${"--min_polish_size " + minPolishSize } \ - ${"--bcftools_path " + bcftoolsPath } \ - ${"--min_component_size " + minComponentSize } \ - ${"--min_dead_end_size " + minDeadEndSize } \ - ${"--contamination " + contamination } \ - ${"--scores " + scores } \ - ${"--low_score " + lowScore } + ~{"--short1 " + short1} \ + ~{"--short2 " + short2} \ + ~{"--unpaired " + unpaired} \ + ~{"--long " + long} \ + --out ~{out} \ + ~{"--min_fasta_length " + minFastaLength} \ + ~{"--keep " + keep } \ + ~{true="--vcf" false="" vcf } \ + ~{"--threads " + threads } \ + ~{"--mode " + mode } \ + ~{"--min_bridge_qual " + minBridgeQual } \ + ~{"--linear_seqs " + linearSeqs } \ + ~{"--spades_path " + spadesPath } \ + ~{true="--no_correct" false="" noCorrect } \ + ~{"--min_kmer_frac " + minKmerFrac } \ + ~{"--max_kmer_frac " + maxKmerFrac } \ + ~{"--kmer_count " + kmerCount } \ + ~{"--depth_filter " + depthFilter } \ + ~{true="--no_miniasm" false="" noMiniasm } \ + ~{"--racon_path " + raconPath } \ + ~{"--existing_long_read_assembly " + existingLongReadAssembly } \ + ~{true="--no_rotate" false="" noRotate } \ + ~{"--start_genes " + startGenes } \ + ~{"--start_gene_id " + startGeneId } \ + ~{"--start_gene_cov " + startGeneCov } \ + ~{"--makeblastdb_path " + makeblastdbPath } \ + ~{"--tblastn_path " + tblastnPath } \ + ~{true="--no_pilon" false="" noPilon } \ + ~{"--bowtie2_path " + bowtie2Path } \ + ~{"--bowtie2_build_path " + bowtie2buildPath } \ + ~{"--samtools_path " + samtoolsPath } \ + ~{"--pilon_path " + pilonPath } \ + ~{"--java_path " + javaPath } \ + ~{"--min_polish_size " + minPolishSize } \ + ~{"--bcftools_path " + bcftoolsPath } \ + ~{"--min_component_size " + minComponentSize } \ + ~{"--min_dead_end_size " + minDeadEndSize } \ + ~{"--contamination " + contamination } \ + ~{"--scores " + scores } \ + ~{"--low_score " + lowScore } } output { File assemblyFasta = out + "/assembly.fasta" @@ -96,7 +98,7 @@ task unicycler { File log = out + "/unicycler.log" } runtime { - cpu: finalThreads - memory: finalMemory + cpu: select_first([threads]) + memory: select_first([memory]) } } \ No newline at end of file