search

prophyle / prophex

ProPhex – an exact k-mer index using Burrows-Wheeler Transform
MIT License
6 stars 1 forks source link

Find kmer locations original sequence #29

Open Adamtaranto opened 1 year ago

Adamtaranto commented 1 year ago

Hi All,

Is it possible to query the bwa index for all locations of a single kmer using prophex?

karel-brinda commented 1 year ago

No, ProphEx doesn't output this. But you can use bwa fastmap (see eg https://github.com/karel-brinda/simplitigs-supplementary/blob/43e09b2052e64eb92fdbea8ee669becb1489ed0e/experiments/multiple_pangenomes_bwa/Snakefile#L54 for how it can be done).

Adamtaranto commented 1 year ago

I was wondering if prophex had an internal method for querying a bwa index that I might be able to re-purpose.

Thanks for the link I'll check out fastmap.