meryl
meryl copied to clipboard
meryl count assertion failed with large K
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.
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