covered.k: Target capture sensitivity

In hummelma/TEQC: Quality control for target capture experiments

Description

Calculates fraction of target bases covered by at least k reads

Usage

covered.k(coverageTarget, k = c(1, 2, 3, 5, 10, 20))

Arguments

coverageTarget	RleList containing Rle vectors of per-target-base coverages for each chromosome, i.e. coverageTarget output from coverage.target
k	integer vector of k-values for which to show fraction of target bases with coverage >= k

Value

Named vector of same length as k giving the corresponding fractions of target bases achieving coverages >= k

Author(s)

Manuela Hummel [e-mail]

See Also

coverage.target, coverage.hist, coverage.uniformity, coverage.plot, coverage.targetlength.plot

Examples

## get reads and targets
exptPath <- system.file("extdata", package="TEQC")
readsfile <- file.path(exptPath, "ExampleSet_Reads.bed")
reads <- get.reads(readsfile, idcol=4, skip=0)
targetsfile <- file.path(exptPath, "ExampleSet_Targets.bed")
targets <- get.targets(targetsfile, skip=0)

## calculate per-base coverages
Coverage <- coverage.target(reads, targets, perBase=TRUE)

covered.k(Coverage$coverageTarget, k=c(1,10,20))

hummelma/TEQC documentation built on March 22, 2021, 9:45 a.m.