footprintR
provides tools for working with single-molecule footprinting data
in R. These include functions for reading collapsed or read-level data from
files generated by modkit and
representation of such data as R objects, functions to manipulate such objects
and for visualization..
Current contributors include:
footprintR
can be installed from GitHub via the
BiocManager
package:
if (!requireNamespace("BiocManager", quietly = TRUE))
install.packages("BiocManager")
BiocManager::install("fmicompbio/footprintR")
Here is a minimal example to for using footprintR
:
library(footprintR)
infile <- system.file("extdata", "modkit_pileup_1.bed.gz", package = "footprintR")
se <- readBedMethyl(bmfile)
se
assayNames(se)
The return value se
is a SummarizedExperiment
with positions in rows and
samples in columns, and assays Nmod
and Nvalid
with the counts of modified
and total (valid) bases.
The results can be conveniently visualized:
plotRegion(se)