Update long_read.smk

IKIM-Essen · Feb 16, 2022 · 63ae246 · 63ae246
1 parent 095d223
commit 63ae246
Showing 1 changed file with 15 additions and 18 deletions.
diff --git a/workflow/rules/long_read.smk b/workflow/rules/long_read.smk
@@ -121,24 +121,21 @@ rule canu_correct:
         "../envs/canu.yaml"
     threads: 16
     shell:
-        """
-        ( if [ -d {params.outdir} ]; then rm -Rf {params.outdir}; fi &&
-        canu -correct -nanopore {input} -p {wildcards.sample} -d {params.outdir} genomeSize=30k minOverlapLength=10 minReadLength=200 \
-        useGrid=false {params.for_testing} \
-        corMMapMerSize=10 corOutCoverage=50000 corMinCoverage=0 maxInputCoverage=20000 \
-        corOverlapper=minimap utgOverlapper=minimap obtOverlapper=minimap \
-        corConcurrency={params.concurrency} \
-        cormhapConcurrency={params.concurrency} cormhapThreads={params.concurrency} \
-        cormmapConcurrency={params.concurrency} cormmapThreads={params.concurrency} \
-        obtmmapConcurrency={params.concurrency} obtmmapThreads={params.concurrency} \
-        utgmmapConcurrency={params.concurrency} utgmmapThreads={params.concurrency} \
-        redConcurrency={params.concurrency} redThreads={params.concurrency} \
-        ovbConcurrency={params.concurrency} \
-        ovsConcurrency={params.concurrency} \
-        oeaConcurrency={params.concurrency}
-        )
-        > {log} 2>&1
-        """
+        "( if [ -d {params.outdir} ]; then rm -Rf {params.outdir}; fi &&"
+        " canu -correct -nanopore {input} -p {wildcards.sample} -d {params.outdir} genomeSize=30k minOverlapLength=10 minReadLength=200"
+        " useGrid=false {params.for_testing}"
+        " corMMapMerSize=10 corOutCoverage=50000 corMinCoverage=0 maxInputCoverage=20000"
+        " corOverlapper=minimap utgOverlapper=minimap obtOverlapper=minimap"
+        " corConcurrency={params.concurrency}"
+        " cormhapConcurrency={params.concurrency} cormhapThreads={params.concurrency}"
+        " cormmapConcurrency={params.concurrency} cormmapThreads={params.concurrency}"
+        " obtmmapConcurrency={params.concurrency} obtmmapThreads={params.concurrency}"
+        " utgmmapConcurrency={params.concurrency} utgmmapThreads={params.concurrency}"
+        " redConcurrency={params.concurrency} redThreads={params.concurrency}"
+        " ovbConcurrency={params.concurrency}"
+        " ovsConcurrency={params.concurrency}"
+        " oeaConcurrency={params.concurrency})"
+        "> {log} 2>&1"
 
 
 # rule medaka_consensus_reference: