mergedPhred: mergedPhred functions: Retrieving and plotting of phred...

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

Description

The Fastqq objects contain position-wise counted phred values. The mergedPhred function adds the counted values for all FASTQ files together into a single matrix. The matrix then again contains position-wise counted phred values. The mergedPhredQuantiles and plotMergedPhredQuant are analogues to the phredQuantiles and plotPhredQuant functions.

Usage

mergedPhred(object)
mergedPhredQuantiles(object, quantiles)
plotMergedPhredQuant(object, main, ...)

Arguments

object	Fastqq: Object which contains collected values from nFiles FASTQ files.
quantiles	numeric: Vector of quantiles. All values must be in [0,1].
main	character: String wich is used as figure caption. Passed internally to plot function.
...	Optional arguments which are passed to the plot function in plotMergedPhredQuant.

Details

The function adds the phred values from all contained FASTQ data.

Value

mergedPhred returns a matrix with 94 rows and (maxSeqLen + 1) columns. mergedPhredQuantiles returns a data.frame with one row for each given quantile and max(seqLen(.)) columns. plotMergedPhredQuant returns nothing.

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\ Ewing B, Green P Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Research 1998 Vol. 8 No. 3 186-194

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"))
#
ph <- mergedPhred(fq)
ph[25:35, 1:15]
pq <- mergedPhredQuantiles(fq, c(0.25, 0.5, 0.75))
plotMergedPhredQuant(fq)
#

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