Skip to content

Commit

Permalink
fixe syntax error
Browse files Browse the repository at this point in the history
  • Loading branch information
rlegendre committed Jan 19, 2018
1 parent 0a4ea11 commit 9341af7
Showing 1 changed file with 5 additions and 5 deletions.
10 changes: 5 additions & 5 deletions sequana/rules/peak_calling/macs2/macs2.rules
Expand Up @@ -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}""")

0 comments on commit 9341af7

Please sign in to comment.