ENCODE-DCC / atac-seq-pipeline

ENCODE ATAC-seq pipeline
MIT License
389 stars 174 forks source link

Discrepancy between signal p-value and fold change bigwig tracks around the open regions #439

Open AayushGrover opened 8 months ago

AayushGrover commented 8 months ago

Hello,

I tried correlating two bigwig tracks (fold change over control and signal p-value) around the ATAC-seq peaks from the GM12878 ATAC-seq ENCODE experiment and similarly for K562 ATAC-seq ENCODE experiment. To be precise, the following are the accession codes of the files:

However, the rank correlation is only around 0.8. I would expect this to be 0.95 and above. Could you help me understand why this might be the case, especially since both tracks are created from the same bam files and hence, the same control?

FC vs P-val (gm12878) FC vs P-val (k562)