Solving the differential peak calling problem in ChIP-seq data