healthy = system.file("extdata", "Healthy.fasta", package = "MethTargetedNGS")
reference = system.file("extdata", "Reference.fasta", package = "MethTargetedNGS")
methP <- methAlign(healthy,reference)
entMeth <- methEntropy(methP)
plot(entMeth,type="l")
Run the code above in your browser using DataLab