# RBedMethyl
<!-- badges: start -->
[](https://lifecycle.r-lib.org/articles/stages.html#stable)
[](https://github.com/CMG-UA/RBedMethyl/actions/workflows/R-CMD-check.yaml)
[](https://codecov.io/gh/CMG-UA/RBedMethyl)
<!-- badges: end -->
Disk-backed access to nanoporetech modkit bedMethyl files for ONT-scale workflows.
## Disclaimer
I am not affiliated with nanoporetech. This package has been tested with modkit
pileup v0.6.1 and commit
https://github.com/nanoporetech/modkit/tree/481e3c9e7930f3f499eadf1ef441606f33e6881c.
## Installation
### Bioconductor (when available)
```r
if (!requireNamespace("BiocManager", quietly = TRUE)) {
install.packages("BiocManager")
}
BiocManager::install("RBedMethyl")
```
### From source (local checkout)
```r
install.packages("devtools")
devtools::install(".")
```
## Usage
```r
library(RBedMethyl)
lines <- c(
paste("chr1", 0, 1, "m", 0, "+", 0, 1, 0, 10, 0.5, 5, 5, 0, 0, 0, 0, 0, sep = "\t"),
paste("chr1", 10, 11, "m", 0, "+", 10, 11, 0, 20, 0.25, 5, 15, 0, 0, 0, 0, 0, sep = "\t")
)
bed <- tempfile(fileext = ".bed")
writeLines(lines, bed)
# Show available fields
bedMethylFields()
# Read coverage and mod_reads
bm <- readBedMethyl(bed, mod = "m", fields = c("coverage", "mod_reads"))
# Per-site methylation fraction
beta(bm)
# Subset by region and summarize
regions <- GenomicRanges::GRanges(
seqnames = "chr1",
ranges = IRanges::IRanges(start = 1, end = 12)
)
bm_region <- subsetByRegion(bm, regions)
summarizeByRegion(bm, regions)
# Bracket subsetting with GRanges
bm[regions]
```
