dmPlotBasic
$ dmtools stats -i sample1.methratio.dm -o prefix1
# user also can obtained coverage and levels of CG:
$ dmtools stats -i sample1.methratio.dm -o prefix1.cg --context 1
$ python3 dmPlotBasic.py -c prefix1.cover.txt prefix2.cover.txt -o tt.pdf
$ dmtools bodystats -i prefix1.dm --gff gene.gff -o prefix1.bodym \
--printcoverage 1
$ python3 dmPlotBasic.py -f prefix1.bodym.cover.cg prefix2.bodym.cover.cg \
-c prefix1.cover.txt prefix2.cover.txt -o tt.pdf