Error in MACS2 paired end peakcalling

Dear all,

I am having an issue running MACS2 using ‘BAMPE’ option. My current macs2 version is macs2 2.7.1. I tried with the macs2 2.1.2 version as well. In both cases, I got the following error messages. When I tried to use ‘BAM’ option, then the program was running. I did use MACS2 before and haven’t seen this error before. Kindly suggest!

INFO @ Wed, 24 Mar 2021 15:36:06: #1 fragments after filtering in control: 22
INFO @ Wed, 24 Mar 2021 15:36:06: #1 Redundant rate of control: 1.00
INFO @ Wed, 24 Mar 2021 15:36:06: #1 finished!
INFO @ Wed, 24 Mar 2021 15:36:06: #2 Build Peak Model…
INFO @ Wed, 24 Mar 2021 15:36:06: #2 Skipped…
INFO @ Wed, 24 Mar 2021 15:36:06: #3 Call peaks…
Traceback (most recent call last):
File “/Users/DSaha/anaconda2/bin/macs2”, line 622, in
main()
File “/Users/DSaha/anaconda2/bin/macs2”, line 57, in main
run( args )
File “/Users/DSaha/anaconda2/lib/python2.7/site-packages/MACS2/callpeak_cmd.py”, line 264, in run
peakdetect.call_peaks()
File “MACS2/PeakDetect.pyx”, line 116, in MACS2.PeakDetect.PeakDetect.call_peaks (MACS2/PeakDetect.c:1815)
File “MACS2/PeakDetect.pyx”, line 164, in MACS2.PeakDetect.PeakDetect.__call_peaks_w_control (MACS2/PeakDetect.c:2070)
ZeroDivisionError: float division

## FYI - below was my command -
macs2 callpeak -t /Volumes/ThunderBay/JL_DS_folder/ChIPseq_analysis_DS/Histonemodification/H3K4me_analysis/alignment/Bamfiles/sorted_H3K4me_WT_2i_rep1.bam -c /Volumes/ThunderBay/JL_DS_folder/ChIPseq_analysis_DS/WT\ set/Input_2i/2_Input_2i/30M_Input-2_2i.sort.bam --broad -f BAMPE -g mm --broad-cutoff 0.1 -q 0.05 --extsize 200 --nomodel -n WT_H3K4me1_2i_rep1 -B --outdir /Volumes/ThunderBay/JL_DS_folder/ChIPseq_analysis_DS/Histonemodification/H3K4me_analysis/macs2_2021/

Best
-Dhurjhoti