Can call single nucleotide peaks to create region peaks ?

[kir1to455]

Hi amaslan, Apologies if this is a confusing statement, what I mean moreso:
If I have the following data (It's not true. It's just my assumption）
chrom start 6mA Coverage 6mA Methylation
chr8 81501 5 3
chr8 81502 2 0
chr8 81503 5 4
chr8 81507 5 0
chr8 81509 2 2
chr8 81511 4 0
chr8 81513 5 0
chr8 81515 4 4
chr8 81517 4 3
chr8 81519 4 3
chr8 81520 2 0

If I want to get a regional peak, what should I do？
In particular, I want to compare with the data of the ChIP-seq peaks (for example CTCF ChIP peaks,H3k27ac ChIP peaks),they can call peaks from MACS2.

Could you give some advice?

Thank you so much for your help.

[amaslan]

We've just released a python package & command-line tool for analysis of DiMeLo-seq data. Check out https://github.com/streetslab/dimelo (documentation here: https://streetslab.github.io/dimelo/). Hopefully one of those functions addresses your question but let me know if not.