Thanks Bjoern,
I read that wigCorrelate, finds the overlapping features between two wig and calculate base-to-base Pearson's Coefficient. When I correlated the data, some correlation was way too low then expected. So, now, I want create 1000 bp windowed wigs from bam files with igvtools count and then correlate them. Is there a way to create 1000 bp windowed wig from bigWig? I will check out deepTools.
L