Kmer screening function

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.