alignmentStats: Get statistics on alignments

In fmicompbio/QuasR: Quantify and Annotate Short Reads in R

Get statistics on alignments

Description

Get statistics on alignments from bam file or qProject object.

Usage

alignmentStats(x, collapseBySample = TRUE)

Arguments

x	the source of alignment bam files, one of: a character vector with bam files a qProject object
collapseBySample	If TRUE and x is a qProject object, sum counts for bam files with identical sample names.

Details

Internally, alignmentStats queries the bam index files similar to 'idxstats' from samtools. Please note that this does not discriminate for example between primary and secondary alignments. If you need more statistics, see for example quickBamFlagSummary from package Rsamtools.

If x is a qProject object, the auxiliary bam files will not contain any unmapped reads, and the corresponding unmapped counts are calculated by subtracting auxiliary mapped counts from the total reads. The latter correspond to the unmapped counts from the corresponding genome bam files.

Value

A matrix with one row per bam file and three columns ("seqlength", "mapped" and "unmapped").

Author(s)

Anita Lerch, Dimos Gaidatzis and Michael Stadler

See Also

qProject, quickBamFlagSummary from package Rsamtools

Examples

## Not run: 
# see qProject manual page for an example
example(qProject)

## End(Not run)

fmicompbio/QuasR documentation built on March 19, 2024, 5:33 a.m.