meryl icon indicating copy to clipboard operation
meryl copied to clipboard

Use sequence to query meryl db

Open Adamtaranto opened this issue 8 months ago • 1 comments

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.

Adamtaranto avatar Jun 04 '24 07:06 Adamtaranto