calculateBamCoverageByInterval {PureCN} | R Documentation

## Function to calculate coverage from BAM file

### Description

Takes a BAM file and an interval file as input and returns coverage for each interval. Coverage should be then GC-normalized using the `correctCoverageBias` function before determining purity and ploidy with `runAbsoluteCN`. Uses the `scanBam` function and applies low quality, duplicate reads as well as secondary alignment filters.

### Usage

    calculateBamCoverageByInterval(bam.file, interval.file, output.file = NULL,
      index.file = bam.file, keep.duplicates = FALSE)

### Arguments

`bam.file` | Filename of a BAM file.
`interval.file` | File specifying the intervals. Interval is expected in first column in format CHR:START-END. The `gc.gene.file` can be used.
`output.file` | Optionally, write minimal coverage file. Can be read with the `readCoverageFile` function.
`index.file` | The bai index. This is expected without the .bai file suffix, see `?scanBam`.
`keep.duplicates` | Keep or remove duplicated reads.

### Value

Returns total and average coverage by intervals.

### Author(s)

Markus Riester

### See Also

`calculateGCContentByInterval correctCoverageBias runAbsoluteCN`

### Examples

    bam.file <- system.file("extdata", "ex1.bam", package="PureCN",
        mustWork = TRUE)
    interval.file <- system.file("extdata", "ex1_intervals.txt",
        package="PureCN", mustWork = TRUE)

    # Calculate raw coverage from BAM file. These need to be corrected for GC-bias
    # using the correctCoverageBias function before determining purity and ploidy.
    coverage <- calculateBamCoverageByInterval(bam.file=bam.file,
        interval.file=interval.file)

[Package PureCN version 1.8.1 Index]