Commit e60c44c9 authored by Kenzo-Hugo Hillion's avatar Kenzo-Hugo Hillion
Browse files

update paired metaphlan with new outputs

parent 6cad6524
......@@ -18,13 +18,13 @@ rule metaphlan3_paired:
r1 = __metaphlan3_input_r1,
r2 = __metaphlan3_input_r2
output:
__metaphlan3_output
profile = __metaphlan3_output,
bowtie2out = __metaphlan3_output_bowtie2out,
sams = __metaphlan3_output_sams
params:
exec_command = __metaphlan3_exec_command,
modules = __metaphlan3_modules,
input_type = __metaphlan3_input_type,
bowtie2out = "{output_dir}/{sample}.bowtie2.bz2".format(output_dir=__metaphlan3_output_dir, sample="{sample}"),
sams = "{output_dir}/{sample}.sam.bz2".format(output_dir=__metaphlan3_output_dir, sample="{sample}"),
options = __metaphlan3_options
threads:
__metaphlan3_threads
......@@ -32,5 +32,5 @@ rule metaphlan3_paired:
command = []
if params.modules:
command.append("module load {params.modules}")
command.append("{params.exec_command} --nproc {threads} --input_type {params.input_type} -s {params.sams} --bowtie2out {params.bowtie2out} {params.options} {input.r1},{input.r2} {output}")
command.append("{params.exec_command} --nproc {threads} --input_type {params.input_type} -s {output.sams} --bowtie2out {output.bowtie2out} {params.options} {input.r1},{input.r2} {output.profile}")
shell(" && ".join(command))
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment