Question: How to extract reads that match k-mer profiles from a collection of sequences?
0
gravatar for O.rka
4 days ago by
O.rka50
O.rka50 wrote:

Let's say you had 10 draft-genome assemblies from different sources with 100 contigs all together from a particular genus.

Are there any tools that allow you to use this "database" of assemblies to then grab any reads that have even a remotely similar k-mer usage to the "database"?

I know about kneaddata but that is mapping to a very specific reference sequence, I'm looking for a way to extract reads that have similar k-mer usage.

Is there a tool that I can use to do this?

sequencing • 94 views
ADD COMMENTlink modified 4 days ago by genomax57k • written 4 days ago by O.rka50

You can use mash (https://genomebiology.biomedcentral.com/articles/10.1186/s13059-016-0997-x ), sourmash ( https://sourmash.readthedocs.io/en/latest/ ) or bbsketch.sh from BBMap suite (https://sourceforge.net/projects/bbmap/ ) can help you with the classification but I am not sure if they have functionality to extract those reads.

ADD REPLYlink modified 4 days ago • written 4 days ago by genomax57k

https://github.com/will-rowe/hulk/ ?

ADD REPLYlink written 4 days ago by shenwei3564.1k
0
gravatar for genomax
4 days ago by
genomax57k
United States
genomax57k wrote:

cookiecutter (https://github.com/ad3002/Cookiecutter ) seems to do what you need. You will need to test and ascertain.

ADD COMMENTlink written 4 days ago by genomax57k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1586 users visited in the last hour