peak calling on open chromatin using EM