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

Description Usage Arguments Details Value Author(s) References Examples

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

1
2
3
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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
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)
#

seqTools documentation built on Nov. 8, 2020, 5:20 p.m.