healthy = system.file("extdata", "Healthy.fasta", package = "MethTargetedNGS")
tumor = system.file("extdata", "Tumor.fasta", package = "MethTargetedNGS")
reference = system.file("extdata", "Reference.fasta", package = "MethTargetedNGS")
healthy = methAlign(healthy,reference)
tumor = methAlign(tumor,reference)
odd_ratio(healthy,tumor)
Run the code above in your browser using DataLab