plotNucCount: plotNucCount: Plots nucleotide counts from Fastqq objects.

In wokai/seqTools: Analysis of nucleotide, sequence and quality content on fastq files.

Description

The function creates plots from nucleotide counts from Fastqq objects.

Usage

plotNucCount(object, nucs=16, maxx,...)

Arguments

object

Fastqq: Object which contains collected values from nFiles FASTQ files.

nucs	integer(optional): Index of nucleotides for which data is plotted. When value is missing, k-mer counts for all contained data is plotted.

maxx	integer(optional): When given, nucleotide counts are plotted for the first maxx nucleotide positions. This option is used for displaying detailed plots from the first read nucleotide positions (which are sometimes not equally distributed).

...	(currently unused).

Details

Values for i must be in {1,...,nFiles}. The nucs index encodes for IUPAC characters as shown in the following table.

1	A	|	6	R	|	11	M	|	16	N
2	C	|	7	Y	|	12	B	|	17	.
3	G	|	8	S	|	13	D	|	18	-
4	T	|	9	W	|	14	H	|	19	=
5	U	|	10	K	|	15	V	|	20	''

When count values for 'A' are to be plotted, 'nucs' must be =1. When count values for 'GC' are to be plotted, 'nucs' must be c(2,3).

Value

None.

Author(s)

Wolfgang Kaisers

References

Cock PJA, Fields CJ, Goto N, Heuer ML, Rice PM The sanger FASTQ file format for sequences with quality scores and the Solexa/Illumina FASTQ variants. Nucleic Acids Research 2010 Vol.38 No.6 1767-1771

See Also

Fastqq-class

Examples

basedir <- system.file("extdata", package="seqTools")
setwd(basedir)
#
fq <- fastqq(c("g4_l101_n100.fq.gz", "g5_l101_n100.fq.gz"), k=4,
    probeLabel=c("g4", "g5"))
#
plotNucCount(fq)
plotNucCount(fq, 1)
plotNucCount(fq, 1:2)
#

wokai/seqTools documentation built on May 4, 2019, 9:46 a.m.