getCoverage: Get coverage count for positions of interest

In JEFworks/HoneyBADGER: HMM-integrated Bayesian approach for detecting CNV and LOH events using single-cell RNA-seq data

Description

Get coverage count for positions of interest

Usage

getCoverage(gr, bamFile, indexFile, verbose = FALSE)

Arguments

gr	GenomicRanges object for positions of interest
bamFile	bam file
indexFile	bai index file
verbose	Boolean of whether or not to print progress and info

Value

totCount Total coverage count information for each position of interest

Examples

## Not run: 
# Sites of interest (chr1:4600000, chr2:2000)
gr <- GRanges(c('chr1', 'chr2'), IRanges(start=c(4600000, 2000), width=1))
# we can get the coverage at these SNP sites from our bams
path <- '../data-raw/bams/'
files <- list.files(path = path)
files <- files[grepl('.bam$', files)]
cov <- do.call(cbind, lapply(files, function(f) {
    bamFile <- paste0(path, f)
    indexFile <- paste0(path, paste0(f, '.bai'))
    getCoverage(gr, bamFile, indexFile)
}))
colnames(cov) <- files

## End(Not run)

JEFworks/HoneyBADGER documentation built on July 24, 2021, 3:01 p.m.