kmer_screen: Kmer screening function

In alexpiper/taxreturn: An R package for retrieving and curating public DNA barcode reference data

Kmer screening function

Description

Kmer screening function

Usage

kmer_screen(seqs, model, threshold = 0.3, k = 5, quiet = FALSE)

Arguments

seqs	A DNAbin object
model	A Profile Hidden Markov model ("PHMM" object) generated by aphid::derivePHMM to screen the sequences against. An already derived model of COI can be loaded using data("model", package="taxreturn")
threshold	The maximum kmer distance allowed from the reference model before a sequence is discarded.
k	The k-mer size to be used for calculating the distance matrix, parsed to kmer::mbed. Note that high values of k may be slow to compute and use a lot of memory due to the large numbers of calculations required, particularly when the residue alphabet is also large.
quiet	Whether progress should be printed to the console.

alexpiper/taxreturn documentation built on Sept. 14, 2024, 7:56 p.m.