k-mer size support >32?

[EvdH0]

Hi Brian,

Thanks for the rewrite! Is there a way to run meryl with a k-mer size larger than 32? Because running with the following example meryl count and meryl print

>test
ATGAAAATCAAAACTCGCTTCGCGCCAAGCCCAACAGGCTATCTGCA

with k=32 results in :

AAAACTCGCTTCGCGCCAAGCCCAACAGGCTA	1
AAAATCAAAACTCGCTTCGCGCCAAGCCCAAC	1
AAACTCGCTTCGCGCCAAGCCCAACAGGCTAT	1
AAATCAAAACTCGCTTCGCGCCAAGCCCAACA	1
AACTCGCTTCGCGCCAAGCCCAACAGGCTATC	1
AATCAAAACTCGCTTCGCGCCAAGCCCAACAG	1
ACTCGCTTCGCGCCAAGCCCAACAGGCTATCT	1
ATCAAAACTCGCTTCGCGCCAAGCCCAACAGG	1
ATGAAAATCAAAACTCGCTTCGCGCCAAGCCC	1
AGCCTGTTGGGCTTGGCGCGAAGCGAGTTTTG	1
CAGATAGCCTGTTGGGCTTGGCGCGAAGCGAG	1
CGCTTCGCGCCAAGCCCAACAGGCTATCTGCA	1
TCAAAACTCGCTTCGCGCCAAGCCCAACAGGC	1
TCGCTTCGCGCCAAGCCCAACAGGCTATCTGC	1
TTGGGCTTGGCGCGAAGCGAGTTTTGATTTTC	1
TGAAAATCAAAACTCGCTTCGCGCCAAGCCCA	1

but with k=33 results in:

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA	8
CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAC	7

which does not seem right.

The full commands are here

meryl count  k=32 test.fna output test.k32.out.merylCount &> test.k32.out.merylCount.std
meryl count  k=33 test.fna output test.k33.out.merylCount &> test.k33.out.merylCount.std

meryl -VV print test.k32.out.merylCount > test.k32.out.histogram 2> test.k32.out.histogram.std
meryl -VV print test.k33.out.merylCount > test.k33.out.histogram 2> test.k33.out.histogram.std

and attached a tarball with the stderr etc. kmertest.tar.gz

Thanks!

Eric