countFastaKmers: countFastaKmers function: Counts DNA k-mers from (compressed)...

Description Usage Arguments Details Value Note Author(s) References Examples

Description

Reads (compressed) fasta files and counts for DNA k-mers in the sequence.

Usage

1
countFastaKmers(filenames,k=4)

Arguments

filenames

character: Vector of fasta file names. Files can be gz compressed.

k

Length of counted DNA k-mers.

Details

Maximal allowed value for k is 12.

Value

matrix.

Note

The static size of the retured k-mer array is 4^k.

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-177

Examples

1
2
3
4
basedir<-system.file("extdata",package="seqTools")
filename<-file.path(basedir,"small.fa")
#writeFai(filename,"small.fa.fai")
res<-countFastaKmers(filename,k=2)

seqTools documentation built on May 2, 2019, 4:45 p.m.