meryl icon indicating copy to clipboard operation
meryl copied to clipboard

meryl count assertion failed with large K

Open Dmitry-Antipov opened this issue 1 year ago • 1 comments

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 avatar May 02 '24 16:05 Dmitry-Antipov

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

Dmitry-Antipov avatar May 03 '24 18:05 Dmitry-Antipov