Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Sourmash gather returns results whereas search returns nothing #526

Closed
mw55309 opened this issue Aug 3, 2018 · 1 comment
Closed

Sourmash gather returns results whereas search returns nothing #526

mw55309 opened this issue Aug 3, 2018 · 1 comment

Comments

@mw55309
Copy link

mw55309 commented Aug 3, 2018

Sorry if this is a dumb question, maybe I am mis-interpreting the functions.

I have a genome (not reads, the genome). For the purposes of this test, I am using ecoli, so I have ecoli.fa and ecoli.sig. I am using the signature to search against the genbank database you offer for download.

I am back on version 2.0.0a6 because if I upgrade to 2.0.0a9 I have to migrate the index ;-)

sourmash gather produces results:

sourmash gather sourmash/ecoli.sig ../genbank-d2-k31.sbt.json --scaled 1000 -o ecoli.gather 

overlap     p_query p_match
---------   ------- --------
4.5 Mbp     100.0%  100.0%      CP011320.1 Escherichia coli strain SQ...

found 1 matches total;
the recovered matches hit 100.0% of the query

And sourmash search does not:

sourmash search sourmash/ecoli.sig ../genbank-d2-k31.sbt.json --scaled 1000 -o ecoli.search

0 matches:
similarity   match
----------   -----

Is this expected? If so, why? What am I mis-understanding?

Cheers
Mick

@ctb
Copy link
Contributor

ctb commented Dec 26, 2018

Sorry for delay in response - this is hopefully the bug we fixed in #484. If you get a chance to test, please do!

@ctb ctb closed this as completed Dec 26, 2018
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants