| read_bam_file {segmenter} | R Documentation |
bam filesCount the reads in each range of the GRanges object
read_bam_file(file, features, ...)
file |
A string. The path to the file. |
features |
A |
... |
Other arguments passed to |
A matrix
# locate the bam file
bam_file <- system.file("extdata", "randomBam.bam", package = "bamsignals")
# load a granges object
rand_anno <- system.file("extdata",
"randomAnnot.Rdata",
package = "bamsignals")
features <- GenomicRanges::promoters(get(load(rand_anno)))
# count reads in ranges
read_bam_file(bam_file, features)