meryl
meryl copied to clipboard
Use sequence to query meryl db
I want to generate a db of all kmers and their counts for a reference genome using meryl count
, then for thousands of small (~1-5 kbp) sequences I want to extract all kmers and find their counts in the genome kmer db.
Is there a way to provide a short sequence as an argument to meryl to query its kmers against an existing db?
It seems like it would not be efficient to run meryl count
on all of the short seqs and have to clean up the .meryl files between each query.