Biostar Beta. Not for public use.
Does jellyfish 2 support kmer size > 31
2
Entering edit mode
12 months ago
Chen • 880

In Jellyfish's original paper(Marcais et al, Bioinformatics 2011), it says that "counting kmers up to 31 bases in length".

But on the other hand, in other papers such as KMC2 (Deorowicz et al, Bioinformatics 2015), Jellyfish2 is used to count kmer size = 55.

So does Jellyfish support k > 31 now? and How it is supported?(i.e. why in original version it is not supported but now it is supported)?

Thanks.

ADD COMMENTlink
4
Entering edit mode
15 months ago
John 12k
Germany

Probably because they originally used 64bits per kmer, which can store up to 32 bases in 2bit (less one for whatever reason, probably to signal reverse complementation or variable-lengths).

With an update, this man-made constraint was removed.

ADD COMMENTlink
0
Entering edit mode

Thanks again for your kind help.

ADD REPLYlink
3
Entering edit mode
12 months ago
Chen • 880

I will answer this question by myself.

I did some test on kmers greater than 31, and at least for the test dataset, the current version of jellyfish does support kmer length > 31.

ADD COMMENTlink
0
Entering edit mode

As you suggested I cleaned up the thread, please edit your answer to include that you also dived in the code to confirm this. (That part got lost by deleting comments.)
Most of us are really friendly, occasionally sarcastic. So you can always assume the best intentions ;-)

ADD REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1