cp1 <- callpeak(CHIP, CTRL, gsize = 5.2e7, store_bdg = TRUE, name = "run_callpeak_narrow0", outdir = tempdir(), cutoff_analysis = TRUE)
