Bug: Value too large to convert to unsigned char #705
Description
Describe the bug
Unsigned char issue raised when running macs3.
Here is the output
Command line: callpeak -t fragments.sorted.tsv.gz --name frags_peak -f FRAG --outdir ./ -g mm
ARGUMENTS LIST:
name = frags_peak
format = FRAG
ChIP-seq file = ['fragments.sorted.tsv.gz']
control file = None
effective genome size = 2.65e+09
band width = 300
model fold = [5, 50]
qvalue cutoff = 5.00e-02
The maximum gap between significant sites is assigned as the read length/tag size.
The minimum length of peaks is assigned as the predicted fragment length "d".
Larger dataset will be scaled towards smaller dataset.
Range for calculating regional lambda is: 10000 bps
Broad region calling is off
Paired-End mode is on
INFO @ 26 Jun 2025 18:01:20: [34 MB] #1 read fragment files...
INFO @ 26 Jun 2025 18:01:22: [34 MB] #1 read treatment fragments...
INFO @ 26 Jun 2025 18:01:22: [125 MB] 1000000 fragments parsed
INFO @ 26 Jun 2025 18:01:23: [136 MB] 2000000 fragments parsed
INFO @ 26 Jun 2025 18:01:24: [149 MB] 3000000 fragments parsed
INFO @ 26 Jun 2025 18:01:25: [163 MB] 4000000 fragments parsed
INFO @ 26 Jun 2025 18:01:26: [177 MB] 5000000 fragments parsed
INFO @ 26 Jun 2025 18:01:26: [191 MB] 6000000 fragments parsed
INFO @ 26 Jun 2025 18:01:27: [208 MB] 7000000 fragments parsed
INFO @ 26 Jun 2025 18:01:29: [222 MB] 8000000 fragments parsed
Traceback (most recent call last):
File "/home/ubuntu/miniconda3/envs/spipe/bin/macs3", line 1017, in
main()
File "/home/ubuntu/miniconda3/envs/spipe/bin/macs3", line 55, in main
run(args)
File "/home/ubuntu/miniconda3/envs/spipe/lib/python3.12/site-packages/MACS3/Commands/callpeak_cmd.py", line 71, in run
(treat, control) = load_frag_files_options(options)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/ubuntu/miniconda3/envs/spipe/lib/python3.12/site-packages/MACS3/Commands/callpeak_cmd.py", line 351, in load_frag_files_options
treat = tp.build_petrack()
^^^^^^^^^^^^^^^^^^
File "MACS3/IO/Parser.py", line 1531, in MACS3.IO.Parser.FragParser.build_petrack
File "MACS3/IO/Parser.py", line 1556, in MACS3.IO.Parser.FragParser.build_petrack
OverflowError: value too large to convert to unsigned char
MACS version: macs3 3.0.3