Description Objects from the Class Slots Methods Author(s) References Examples
Contains quality related summarizing data on fastq files.
Objects can be created by calls of the form Fastqq("test.fastq")
.
filenames
:"character"
: Vector of Fastqq file names.
probeLabel
:"character"
: Vector of probe labels.
nFiles
:"integer"
: Length of fileNamess.
k
:"integer"
: Length of counted DNA k-mers.
maxSeqLen
:"integer"
Maximum sequence length found in fastq files. Determines row-number in 'seqLenCount' matrix and column-number in 'nac' and 'phred' slot.
kmer
:"matrix"
Matrix containing DNA k-mers counts.
firstKmer
:"matrix"
Matrix containing count of incipient DNA k-mers.
nReads
:"integer"
Vector containing number of reads per file.
seqLenCount
:"matrix"
Matrix containing Counts of read lengths.
gcContent
:"matrix"
Matrix containing GC content (in percent).
nN
:"integer"
Vector containing Number of N nucleotide entries per file.
nac
:"list"
Contains counted per position alphabet frequencies.
phred
:"list"
Contains per position phred count tables (one per Fastqq file).
seqLen
:"matrix"
Contains minimal and maximal sequence length (one column per file).
collectTime
:"list"
Contains start and end time of fastq reading as 'POSIXct'.
The following methods are defined for class Fastqq
:
Basic accessors:
signature(object="Fastqq")
: Returns k-value (length of DNA k-mers) as integer
.
signature(object="Fastqq")
: Returns matrix
with 4^k rows anc nFiles
columns. For each k-mer and fastq-file, the absolute count value of the k-mer in the fastq file is given.
signature(object="Fastqq")
: Returns number of Files from which data has been collected as integer
.
signature(object="Fastqq")
: Returns integer
vector of length nFiles
. For each fastq file, the absolute number of containes 'N' nucleotide entries is given.
signature(object="Fastqq")
: Returns number of reads in each fastq file as integer
.
signature(object="Fastqq")
: Returns number names of fastq files from which data has been collected as character
.
signature(object="Fastqq")
: Returns maximum sequence length which has been found in all fastq files as integer
.
signature(object="Fastqq")
: Returns matrix which tables counted read length in all fastq files.
signature(object="Fastqq",i="numeric")
: Returns integer
vector of length 100 which countains absolute read count numbers for each percentage of GC-content. i
is the index of the fastq file for wich the values are returned. The GC content values for all files together can be obtained using gcContentMatrix
.
signature(object="Fastqq",i="integer")
: Returns matrix
which contains the absolute nucleotide count values for each nucleotide and read position. i
is the index of the fastq file for wich the values are returned.
signature(object="Fastqq")
: Returns matrix
with two rows and nFiles
columns. For each file the minimum and maximum read length is given.
signature(object="Fastqq")
: Returns a matrix
with 4^k rows and nFiles
columns. Each entry gives the absolute count of the k-mer (given as row name) in each file (given as column name).
signature(object="Fastqq",i="integer")
: Returns a matrix
with 93 rows and maxSeqLen
columns. The matrix gived the absolute counts of each phred value for each sequence position. i
is the index of the fastq file for wich the values are returned.
signature(object="Fastqq", quantiles="numeric", i="integer")
: Returns a data.frame
. The data.frame has one row for each given quantile and maxSeqLen
columns. Each value gives the quantile (given by row name) of the phred values at the sequence position (given by column name). For the quantiles
argument, a numeric vector with values in [0,1] must be given. For the i
argument, a single integer value must be given which denotes the index of the fastq file from which values are returned (value must be in {1,...,nFiles}).
signature(object="Fastqq")
: Returns character
vector which contains the probeLabel
entries for given Fastqq
object.
Wolfgang Kaisers
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
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 | 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"))
fileNames(fq)
getK(fq)
nNnucs(fq)
nFiles(fq)
nReads(fq)
maxSeqLen(fq)
collectTime(fq)
collectDur(fq)
slc<-seqLenCount(fq)
nf<-nucFreq(fq,1)
nf[1:4,1:10]
seqLen(fq)
probeLabel(fq)
probeLabel(fq)<-1:nFiles(fq)
#
kc<-kmerCount(fq)
kc[1:10,]
plotKmerCount(fq)
#
ph<-phred(fq,1)
ph[25:35,1:15]
pq<-phredQuantiles(fq,c(0.25,0.5,0.75),1)
plotNucFreq(fq,1)
# Nucleotide count
plotNucCount(fq,2:3)
# GC content
gcContent(fq,1)
#
fqq<-fq[1]
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.