I am trying to run a macs2 peak analysis on bam files generated from unpaired reads. Format = AUTO, with no additional options invoked. The warning I receive is "Too few paired peaks (0) so I can not build the model! Broader your MFOLD range parameter may erase this error..."
Re running the script with --nomodel and a suggested --extsize provides strange data that does not line up with what I can see on IGV.
I suspect the issue is that macs is trying to build a model with paired reads when I haven't provided paired data - is this possible? Should I be adjusting my usage parameters?