-
Notifications
You must be signed in to change notification settings - Fork 80
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
[MRG] add an lca gather
signature decomposition approach
#390
Conversation
lca gather
signature decomposition approachlca gather
signature decomposition approach
Codecov Report
@@ Coverage Diff @@
## master #390 +/- ##
==========================================
- Coverage 90.14% 90.07% -0.08%
==========================================
Files 31 32 +1
Lines 4586 4824 +238
Branches 36 36
==========================================
+ Hits 4134 4345 +211
- Misses 451 478 +27
Partials 1 1
Continue to review full report at Codecov.
|
@@ -15,7 +15,7 @@ | |||
from sourmash_lib.lca.command_index import load_taxonomy_assignments | |||
|
|||
def compare_csv(args): | |||
p = argparse.ArgumentParser() | |||
p = argparse.ArgumentParser(prog="sourmash lca compare_csv") |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Oh, we can make a low-hanging fruit issue for this change (on other commands)
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Punted to #415
from sourmash_lib.lca.lca_utils import debug, set_debug | ||
from sourmash_lib.search import format_bp | ||
|
||
LCAGatherResult = namedtuple('LCAGatherResult', |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I have a feeling that we might be reaching a point where we want a Result
class to collect functionality from all these NamedTuples
...
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Punted to #416
yay w00t!
|
Ref #384. Unlike 'sourmash gather' this requires taxonomy information from an LCA database.
It is also
Unlike
sourmash gather
this will not print out arbitrary genome names or anything - it's solely about taxonomic display. It also alerts the user to the number of equivalent matches and outputs it to the CSV (but only displays one match).other additions:
sourmash gather --output-unassigned
sourmash lca
subcommands where subcommand was not output to helpsourmash lca classify
db error message isn't very friendly #387 - fail more nicely when failing to parse LCA DBschecklist:
make test
Did it pass the tests?make coverage
Is the new code covered?without a major version increment. Changing file formats also requires a
major version number increment.
changes were made?