homopolymerHist: Create A Histogram Of The Homopolymer Stretches

In R453Plus1Toolbox: A package for importing and analyzing data from Roche's Genome Sequencer System

Description

This function creates a histogram for the different lengths of the homopolymer stretches with one bar for each nucleotide in the flow. With the height giving the number of this homopolymer stretch in the flowgram.

Usage

  homopolymerHist(x, range=c(0, length(flowgram(x))), xlab="Homopolymer length", 
    ylab="Number of homopolymers", col=c(A="black", C="red", G="blue", T="green"), ...)

Arguments

x	An object of class SFFRead.
range	Two positions between which the flows should be plotted.
xlab	The X axis label.
ylab	The Y axis label.
col	The colors in which the four nucleotids should be plotted.
...	Arguments to be passed to methods, such as graphical parameters (see ‘par’).

Author(s)

Christian Ruckert

R453Plus1Toolbox documentation built on Nov. 8, 2020, 5:59 p.m.