diff --git a/sequana/rules/peak_calling/macs2/macs2.rules b/sequana/rules/peak_calling/macs2/macs2.rules index 00487181a..1f1c50666 100644 --- a/sequana/rules/peak_calling/macs2/macs2.rules +++ b/sequana/rules/peak_calling/macs2/macs2.rules @@ -16,23 +16,23 @@ rule macs2: """ input: - bam = __macs2__input_treatment, - ctl = __macs2__input_control + bam = __macs2__input_treatment params: gsize = config["macs2"]["genomeSize"], readLength = config["macs2"]["readLength"], qval = config["macs2"]["qval"], options = config["macs2"]['options'], - prefix = __macs2__output + prefix = __macs2__output, + ctl = __macs2__input_control log: __macs2__log output: touch(__macs2__output_done) threads: 1 run: - if input.ctl == "None": + if input.ctl == None: shell(""" macs2 callpeak -t {input.bam} -f BAM -g {params.gsize} -n {params.prefix} --extsize={params.readLength} -q {params.pval} {params.options} 2> {log}""") else: - shell(""" macs2 callpeak -t {input.bam} -c {input.ctl} -f BAM -g {params.gsize} -n {params.prefix} + shell(""" macs2 callpeak -t {input.bam} -c {params.ctl} -f BAM -g {params.gsize} -n {params.prefix} --extsize={params.readLength} -q {params.pval} {params.options} 2> {log}""")