Generate mappable regions of the genome for a given K-mer length.
Knowing where reads can uniquely map in the genome is useful for nascent RNA assays, both in statistical calculations and to make predictions. kmap takes two inputs: a genome and the read length, K. The output is automatically cached on disk.
library(kmap)
gr <- mappable("hg38", kmer = 36)
Install the latest version of kmer using:
BiocManager::install("coregenomics/kmap")
If the above command fails, install Bioconductor.
install.packages("BiocManager")