diff --git a/src/pipelines/nf_modules/peak_calling.nf b/src/pipelines/nf_modules/peak_calling.nf index 07f0b942d99b982616c04d733f9392f13ecb2d69..3e6645c338c316cd422f0f9788146489e0cc6836 100644 --- a/src/pipelines/nf_modules/peak_calling.nf +++ b/src/pipelines/nf_modules/peak_calling.nf @@ -86,7 +86,7 @@ process bedtools_slop { ${new_id}.enlarged.tmp.bed sort -k 1,1 -k2,2n ${new_id}.enlarged.tmp.bed > ${new_id}.enlarged.sort.bed bedtools merge -i ${new_id}.enlarged.sort.bed \ - -d 6 -c 5,6 -o mean,distinct | awk '{print \$1"\\t"\$2"\\t"\$3"\\tpeak_"NR"\\t"\$4"\\t"\$5}' > ${new_id}.merged.bed + -c 5,6 -o mean,distinct | awk '{print \$1"\\t"\$2"\\t"\$3"\\tpeak_"NR"\\t"\$4"\\t"\$5}' > ${new_id}.merged.bed """ } @@ -135,7 +135,7 @@ process bedtools_merge { """ cat *.bed | sort -k1,1 -k2,2n > tmp.bed bedtools merge -i tmp.bed \ - -d 6 -c 1,5,6 -o count,mean,distinct | awk '{if(\$4 >= 2){print \$1"\\t"\$2"\\t"\$3"\\tpeak_"NR"\\t"\$5"\\t"\$6}}' > ${sp}.rep_merged.bed + -c 1,5,6 -o count,mean,distinct | awk '{if(\$4 >= 2){print \$1"\\t"\$2"\\t"\$3"\\tpeak_"NR"\\t"\$5"\\t"\$6}}' > ${sp}.rep_merged.bed """ }