coverage.target: Calculates read coverage
In TEQC: Quality control for target capture experiments

Description Usage Arguments Value Author(s) See Also Examples

Calculates average coverage over all target bases, average coverage for each target separately, and per-base coverage for all and for targeted bases

1	coverage.target(reads, targets, Offset = 0, perTarget = TRUE, perBase = TRUE)

`reads`	`RangedData` table containing positions of sequenced reads, i.e. output from `get.reads`
`targets`	`RangedData` table containing positions of target regions, i.e. output from `get.targets`
`Offset`	integer; add `Offset` bases on both sides to targeted regions and potentially collapse resulting overlapping target regions
`perTarget`	if TRUE, coverage average and standard deviation per target are calculated and returned
`perBase`	if TRUE, the per-base coverages i) only for targeted bases and ii) for all sequenced and/or targeted bases, are returned

A list is returned with elements

`avgTargetCoverage`	average coverage over all target bases
`targetCoverageSD`	standard deviation of coverage of all target bases
`targetCoverageQuantiles`	0% (minium), 25%, 50% (median), 75% and 100% (maximum) quantiles of coverage of all target bases
`targetCoverages`	Input `RangedData` table `targets` with two additional 'values' columns `avgCoverage` and `coverageSD`. The former contains the average coverage for each target, the latter the respective coverage standard deviation. Only returned if `perTarget` equals `TRUE`.
`coverageAll`	`RleList` containing a `Rle` vector for each chromosome with coverages for all bases that are sequenced and/or within a targeted; only returned if `perBase` equals `TRUE`
`coverageTarget`	`RleList` containing a `Rle` vector for each chromosome with coverages for target bases only; only returned if `perBase` equals `TRUE`

Manuela Hummel m.hummel@dkfz.de

covered.k, coverage.hist, coverage.uniformity, coverage.plot, coverage.targetlength.plot

## 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)

## total average, per-base and per-target coverages
Coverage <- coverage.target(reads, targets)

Loading required package: BiocGenerics
Loading required package: parallel

Attaching package: 'BiocGenerics'

The following objects are masked from 'package:parallel':

    clusterApply, clusterApplyLB, clusterCall, clusterEvalQ,
    clusterExport, clusterMap, parApply, parCapply, parLapply,
    parLapplyLB, parRapply, parSapply, parSapplyLB

The following objects are masked from 'package:stats':

    IQR, mad, sd, var, xtabs

The following objects are masked from 'package:base':

    Filter, Find, Map, Position, Reduce, anyDuplicated, append,
    as.data.frame, cbind, colMeans, colSums, colnames, do.call,
    duplicated, eval, evalq, get, grep, grepl, intersect, is.unsorted,
    lapply, lengths, mapply, match, mget, order, paste, pmax, pmax.int,
    pmin, pmin.int, rank, rbind, rowMeans, rowSums, rownames, sapply,
    setdiff, sort, table, tapply, union, unique, unsplit, which,
    which.max, which.min

Loading required package: IRanges
Loading required package: S4Vectors
Loading required package: stats4

Attaching package: 'S4Vectors'

The following object is masked from 'package:base':

    expand.grid

Loading required package: Rsamtools
Loading required package: GenomeInfoDb
Loading required package: GenomicRanges
Loading required package: Biostrings
Loading required package: XVector

Attaching package: 'Biostrings'

The following object is masked from 'package:base':

    strsplit

Loading required package: hwriter
[1] "read 19546 sequenced reads"
[1] "read 50 (non-overlapping) target regions"