I have a WIG file containing read coverage data from ChIP-seq peak calling analysis.
For a given set of coordinates is there a simple fast way of determining the 'summit' coordinate, i.e the coordinate with the highest read pileup.
My own method would be slow as i would run through the (big) WIG file, extract the region of interest and find the coordinate with the highest score. But surely there must be a quicker way of doing this? If not fair enough.
The peak caller does report the summit from its own data, but i want to check out different regions.