macs3-project/MACS

Q: using macs3 to detect peaks on long-read readings

Opened this issue · 0 comments

I wanted to test macs3 for genomic peaks detection. For that, a number of reads obtained from Pacbio, previously aligned and subsequently sorted and indexed has been provided for macs3.
The line of code is that.
macs3 callpeak -t /home/aligment/1A_.bam -g hs -f BAM --outdir /home/aligment/-n 1A_peaks

Traceback (most recent call last):
File "/conda/envs/macs3_env/bin/macs3", line 1028, in
main()
File "/conda/envs/macs3_env/bin/macs3", line 49, in main
run( args )
File "/conda/envs/macs3_env/lib/python3.12/site-packages/MACS3/Commands/callpeak_cmd.py", line 70, in run
(treat, control) = load_tag_files_options (options)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/conda/envs/macs3_env/lib/python3.12/site-packages/MACS3/Commands/callpeak_cmd.py", line 364, in load_tag_files_options
ttsize = tp.tsize()
^^^^^^^^^^
File "MACS3/IO/Parser.pyx", line 1131, in MACS3.IO.Parser.BAMParser.tsize
File "MACS3/IO/Parser.pyx", line 1160, in MACS3.IO.Parser.BAMParser.tsize
struct.error: unpack requires a buffer of 4 bytes