-
Notifications
You must be signed in to change notification settings - Fork 6
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
Feature request: Utilize Diamond's contig features #10
Comments
Thanks for using the software and for making the suggestion! Yes most of the call to diamond is made with default parameters, and configurable options are so far only I briefly remember thinking about the Again, thanks for making the suggestion. It's interesting to think about and discuss these things. Adding the range-culling feature as an option to contigtax should not be a problem as it doesn't appear to affect the output format and thus doesn't cause problems with downstream assignments. I noticed that the feature requires diamond to be run with a frameshift penalty which is mostly recommended for error prone long-read sequence output and not the assembled short-read sequences I had in mind when designing contigtax, but nevertheless it may have it's benefits. I've found that the |
Thanks for the awesome software! If I understand the code correctly, Diamond is executed using largely default parameters. I'd suggest adding in
‐‐range‐culling ‐‐top 10 -F 15
(source), but this will likely require rewrites of other areas of contigtax. These parameters will perform local Diamond alignment, retaining the top hit (within 10%) in each area of the query contig. We'd then have to not filter by bitscore in contigtax, and also should rely on the evalue parameter of Diamond instead of filtering on that. Just wanted to get the discussion started - there's probably a bunch of other design decisions that I'm missing.The text was updated successfully, but these errors were encountered: