Skip to content
Snippets Groups Projects
Commit d8f4e72b authored by Valentina Galata's avatar Valentina Galata
Browse files

assembly w/ operams: added perl for execution (issue #40), minor change of rule name

parent e2acad95
No related branches found
No related tags found
1 merge request!76Merge "cleanup" branch with "master" branch
......@@ -145,7 +145,7 @@ rule assemble_hy_metaspades:
"cd $(dirname {output}) && ln -sf contigs.fasta $(basename {output}) && "
"date) 2> {log.err} > {log.out}"
rule assemble_hy_operams:
rule assembly_hy_operams:
input:
lr=os.path.join(RESULTS_DIR, "preproc/metag/lr/lr.fastq.gz"),
r1=os.path.join(RESULTS_DIR, "preproc/metag/sr/R1.fastp.fastq.gz"),
......@@ -166,7 +166,7 @@ rule assemble_hy_operams:
shell:
"(date && "
"zcat {input.lr} > {output.lr} && "
"{config[operams][bin]} --short-read1 {input.r1} --short-read2 {input.r2} --long-read {output.lr} --contig-file {input.asm} "
"perl {config[operams][bin]} --short-read1 {input.r1} --short-read2 {input.r2} --long-read {output.lr} --contig-file {input.asm} "
"--out-dir $(dirname {output.asm}) --long-read-mapper minimap2 --num-processors {threads} && "
"cd $(dirname {output.asm}) && ln -sf contigs.fasta $(basename {output.asm}) && "
"date) 2> {log.err} > {log.out}"
\ No newline at end of file
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