Histogram of match length distribution of reads.

Share:

Description

Histogram of match length distribution of reads.

Usage

1
histMatchLength(aln, log10Transf = 0, titleHist)

Arguments

aln

A GAlignments object of the BAM mapping file.

log10Transf

A boolean. Either 0 (default) or 1 (log10).

titleHist

a character. The main title for the histogram. Default - none.

Value

A list with 2 elements. The first element: a data.frame of the number of counts per match length distribution. The second element in the list: a ggplot2 histogram of the match length distribution.

Examples

1
2
3
4
5
6
7
#starting from a GAlignment object
data(ctrlGAlignments)
aln <- ctrlGAlignments
#no log10 scaling
matchLenDistr <- histMatchLength(aln, 0)
#to plot the histogram
matchLenDistr[[2]]