struct.error in MACS2 callpeak

0

Hi all,

I am trying to callpeaks on sorted and subsampled BAM files from the results of my ChIP-seq program. When I do callpeaks, I get the following error after the read treatment tags step, and I have no idea what it means/how to solve it:

INFO  @ Sun, 23 May 2021 10:20:55: #1 read treatment tags... 
struct.error: unpack requires a buffer of 4 bytes
Exception ignored in: 'MACS2.IO.Parser.BAMParser.tsize'
struct.error: unpack requires a buffer of 4 bytes
Traceback (most recent call last):
  File "/dartfs-hpc/rc/home/2/f0036r2/.conda/envs/macs2_use/bin/macs2", line 653, in <module>
    main()
  File "/dartfs-hpc/rc/home/2/f0036r2/.conda/envs/macs2_use/bin/macs2", line 51, in main
    run( args )
  File "/dartfs-hpc/rc/home/2/f0036r2/.conda/envs/macs2_use/lib/python3.6/site-packages/MACS2/callpeak_cmd.py", line 65, in run
    else:       (treat, control) = load_tag_files_options  (options)
  File "/dartfs-hpc/rc/home/2/f0036r2/.conda/envs/macs2_use/lib/python3.6/site-packages/MACS2/callpeak_cmd.py", line 391, in load_tag_files_options
    treat = tp.build_fwtrack()
  File "MACS2/IO/Parser.pyx", line 1169, in MACS2.IO.Parser.BAMParser.build_fwtrack
  File "MACS2/IO/Parser.pyx", line 1181, in MACS2.IO.Parser.BAMParser.build_fwtrack
  File "MACS2/IO/Parser.pyx", line 1166, in MACS2.IO.Parser.BAMParser.get_references
struct.error: unpack requires a buffer of 4 bytes

Thanks in advance!


ChIP-seq


MACS2

• 29 views



Source link