Skip to content
Snippets Groups Projects
Verified Commit c4218fde authored by Laurent Modolo's avatar Laurent Modolo
Browse files

macs3: fix call cmd

parent 79984ea1
No related branches found
No related tags found
No related merge requests found
......@@ -18,7 +18,7 @@ process peak_calling {
script:
/* remove --nomodel option for real dataset */
"""
macs2 callpeak \
macs3 callpeak \
--treatment ${bam_ip} \
--call-summits "True"\
--control ${bam_control} \
......@@ -53,7 +53,7 @@ awk '{print \$1"\t"\$2"\t"\$3"\t.\t+\t"\$4}' ${bg_ip} > \
${bg_ip.simpleName}.bed
awk '{print \$1"\t"\$2"\t"\$3"\t.\t+\t"\$4}' ${bg_control} > \
${bg_control.simpleName}.bed
macs2 callpeak \
macs3 callpeak \
--treatment ${bg_ip.simpleName}.bed \
--call-summits "True"\
--control ${bg_control.simpleName}.bed \
......
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