## Coverage

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@ Statistics table

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@ Sensitivity table

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@