search

marbl / meryl

A genomic k-mer counter (and sequence utility) with nice features.
113 stars 13 forks source link

meryl count assertion failed with large K #44

Open Dmitry-Antipov opened 4 months ago

Dmitry-Antipov commented 4 months ago

Tried to run with relatively large k (101), got this assertion

meryl count k=101 compress hifi.fasta.gz output 101.hpc.meryl
<...>
Start counting with THREADED method.                                                                                                                                                            
meryl: meryl/merylCountArray.C:512: uint64 merylCountArray::add(kmdata): Assertion `wordEnd <= 192' failed.

k=71 seems to work

would be nice to have max allowed k somewhere in usage/docs.

Dmitry-Antipov commented 4 months ago

actually k=71 worked only in the sense that it did not crash. Results are irrelevant unique 0 (exactly one instance of the kmer is in the input) distinct 1458 (non-redundant kmer sequences in the input) present 33707750420 (...)

k=63 provided reasonable results