I am running diffbind using narrowPeaks output from macs2
When using dba.count, I am taking each summit and counting reads +/-250 from the summit so I'm expecting to see in the diffbind report that the width column only contains "501"
bindingmatrix <- dba.count(dbaobject, summits=250) # summit +/- 250 bases
However, around 2% of the rows in the report contain a range of widths up to ~1400
Can anyone explain why this is?