The function 'coverage.target()' calculates the read coverage for each base that is sequenced (and aligned to the reference) and/or located in a target region. Coverage average and standard deviation are calculated for the target as a whole and for each target region separately. In addition to per-base calculation of read coverage, the numbers of reads falling on each target are given as an alternative measure for target coverage. A read is considered on-target, if it overlaps the target by at least one base.
Average on-target coverage
> Coverage <- coverage.target(reads, targets, Offset=@OFFSET@)@AVGCOV@
Coverage per target
> targetcov <- Coverage$targetCoverages > targetcov <- readsPerTarget(reads, targets2, Offset=@OFFSET@)@COVTARG@ Full table
Sensitivity
A high sensitivity of the target capture is shown by a large fraction of target bases that are covered by at least 1, 2, 3... reads.
> covered.k(Coverage$coverageTarget)@SENSITIVITY@
Coverage histogram
Target coverage can be visualized by the coverage histogram. A line is added that shows the cumulative fraction of target bases (right y-axis) with a coverage of at least the corresponding x-axis value. The line represents the results of the previously called function 'covered.k()' for all possible values of k. Additionally, we highlight with dashed lines the target base fraction achieving a coverage of at least 8, in this case.
> coverage.hist(Coverage$coverageTarget, covthreshold=8)@COV_HIST@