Thanks for the insights and ideas Fidel, I will definitely try to make a heatmap instead of just profile plots from now on.
Most of our profile plots are always biased towards the 3' end than 5'. I did try stratifying based on gene length to see the difference, but I still find that the PolII is preferentially higher in the 3' end regardless of gene lengths for my case.
We did use MACS2 (since none of the other broad peak callers were fruitful) for both narrow and broad peaks, and surprisingly we find lot more narrow than broad peaks with this data. I am not sure if this is due to depth, but we have a average of 20X coverage for 23MB genome.
Thanks,
Aarthi