This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to quickly find the position of k-mer on the genome

Hi, all

How to quickly find the position of k-mers on the genome? Currently, I use software such as bwa and fuzznucl, but I feel it is not fast enough.

alignment sequencing k-mer
• 3,139 views • link
updated by Christophe • Christophe • written by slbai00001 • slbai00001

Are you looking to find one specific k-mer or a map of all k-mers?

• 0 views • link
written by Malte Petersen • Malte Petersen

In fact, I first used kmc to obtain all the specific k-mers of one genome relative to the other genome, and then I wanted to obtain the position of these k-mers on the genome. The number of k-mers was about several million or even a few giga.

• 0 views • link
written by slbai00001 • slbai00001

Yes. Can I ask for your penultimate goal? Unless you're going to write a new assembler, I'm not sure what you would need that map for.

• 0 views • link
written by Malte Petersen • Malte Petersen

Did you find an efficient solution? I would like to perform the same task.

• 0 views • link
written by Christophe • Christophe

1 answer

Something similar has been asked here. In essence, you can use motif finding tools such as fuzznuc from the EMBOSS suite, BLAT with a small tile size and options -minMatch and -minIdentity, or seqkit locate, depending on your input/output format requirement.

But yes, as you already noticed, it will be slow because there are lots of k-mers. You could also format your k-mer list in fastq/fasta format and use BWA with zero mismatches allowed.

• 0 views • link
written by Malte Petersen • Malte Petersen

Log in to answer this question.

More posts like this