Hello, would it be possible+reasonable to enter csaw peak calling with RleList object such as those generated by the GenomicAlignments' coverage method? I'm not 100% sure but. I guess it'd be more performant to store and load the RleList objects compared to bam files. regards, Arne [[alternative HTML version deleted]]