I want to estimate the methylation coverage in the gene region. I have an output file from bismark. I would like to get the methylation coverage ( % of methylated cytosines/total cytosines) in the genomic region. I am not sure if methylpy can handle the output from bismark. Does anyone have suggestion?