https://github.com/nf-core/methylseq/issues/632 Regarding this issue I opened (and the corresponding PR here https://github.com/nf-core/methylseq/pull/633), any idea from your side whether there is a specific reason why Preseq is running on deduplicated instead of raw aligned BAM files? Also happy to help if needed!