Verified Commit a97de040 authored by Aurélien Ginolhac's avatar Aurélien Ginolhac 🚴
Browse files

revert macs2 hard-coded settings

parent 88cce103
......@@ -58,7 +58,7 @@ rule macs2_callpeak_broad:
log:
"logs/macs2/callpeak.{sample}-{control}.broad.log"
params: # --nomodel --extsize 147 added for fake data
lambda w, input:"--broad-cutoff 0.1 -f {bam_format} {gsize} --nomodel --extsize 147 -B --SPMR --keep-dup all {pvalue} {qvalue}".format(
lambda w, input:"--broad-cutoff 0.1 -f {bam_format} {gsize} -B --SPMR --keep-dup all {pvalue} {qvalue}".format(
gsize="{}".format(open(input.gsizepath).read().strip()),
pvalue="-p {}".format(config["params"]["callpeak"]["p-value"]) if config["params"]["callpeak"]["p-value"] else "",
qvalue="-q {}".format(config["params"]["callpeak"]["q-value"]) if config["params"]["callpeak"]["q-value"] else "",
......
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