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

macs2: set default qvalue to 0.99

parent 5dc7c7f8
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
src/nf_modules/macs2/main.nf
+ 2
− 0
View file @ 6ddf33f5
...@@ -30,6 +30,7 @@ macs2 callpeak \ ...@@ -30,6 +30,7 @@ macs2 callpeak \
--call-summits \ --call-summits \
--control ${bam_control} \ --control ${bam_control} \
--keep-dup all \ --keep-dup all \
--qvalue 0.99 \
--name ${bam_ip.simpleName} 2> \ --name ${bam_ip.simpleName} 2> \
${bam_ip.simpleName}_macs2_report.txt ${bam_ip.simpleName}_macs2_report.txt
...@@ -68,6 +69,7 @@ awk '{print \$1"\t"\$2"\t"\$3"\t.\t+\t"\$4}' ${bg_control} > \ ...@@ -68,6 +69,7 @@ awk '{print \$1"\t"\$2"\t"\$3"\t.\t+\t"\$4}' ${bg_control} > \
macs2 callpeak \ macs2 callpeak \
${params.peak_calling_bg} \ ${params.peak_calling_bg} \
--treatment ${bg_ip.simpleName}.bed \ --treatment ${bg_ip.simpleName}.bed \
--qvalue 0.99 \
--call-summits \ --call-summits \
--control ${bg_control.simpleName}.bed \ --control ${bg_control.simpleName}.bed \
--keep-dup all \ --keep-dup all \
......
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