Using Modkit pileup --modified base C

[serignebasse]

Can someone explain me How Modkit compute Nmod (for 5mC). They do it per reads all they use all reads aligned on same position and count number of read with modified C.

• Nmod - Number of calls passing filters that were classified as a residue with a specified base modification.