This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Extracting k-mer counts from multiple genome sequence files

Good day everyone, I am new here.

So, I have downloaded 8 completed genome fasta files for 8 strains of Bacillus subtilis spp.

My aim is to do classification based on their k-mer abundance profiles.

I am wondering, is there any tools that I can use to generate and extract the k-mer counts for each of the 8 genome files in a single output?

k-mer
• 2,943 views • link
written by yas • yas

3 answers

Two very popular tools specifically designed for k-mer counting:

  1. KMC3 [tool] [paper]

  2. Jellyfish [tool] [paper]

Other k-mer counting tools are listed in the benchmark study of k-mer counting methods.

• 0 views • link
written by Andrzej Zielezinski Andrzej Zielezinski

Thank you!

• 0 views • link
written by yas • yas

Mash is an excellent tool for this kind of thing. It's far more sophisticated than simple k-mer abundance counting..

• 0 views • link
written by 5heikki 5heikki

Thank you so much for the answer. But I need to generate the k-mer abundance profiles.

• 0 views • link
written by yas • yas

You could use kmercountmulti.sh from the BBTools Suite if you are specifically interested in k-mers.

• 0 views • link
written by Matthias Zepper Matthias Zepper

Thanks for the suggestion.

• 0 views • link
written by yas • yas

Log in to answer this question.

More posts like this