search

gmarcais / Jellyfish

A fast multi-threaded k-mer counter
Other
471 stars 136 forks source link

file run several times gives different results #146

Open mamora opened 5 years ago

mamora commented 5 years ago

I am running a file with several Fasta sequences with Jellyfish v. 2.2.1.

I get different results for Kmer counts for 15 5mers. Sometimes there is up to 2fold difference in numbers

some examples:

kmers | run1 | run2 | run3 TTTAG | 54884 | 22116 | 54884 GGTTT | 54977 | 41921 | 41921 TGTGT | 53702 | 41158 | 53702 TTATG | 54518 | 53494 | 54518 TATGT | 21791 | 21663 | 21663

This is my code: jellyfish count -m $i -o positive_controlinput$i.fasta -s 100000000 -t 10 I_S9_R1.txt.fasta

gmarcais commented 5 years ago

This is fixed in 2.3.0.

mamora commented 5 years ago

thanks