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

"No hits found" in blast search breaks hash sorting #20

Closed
2 tasks done
bretonics opened this issue Jan 17, 2017 · 1 comment
Closed
2 tasks done

"No hits found" in blast search breaks hash sorting #20

bretonics opened this issue Jan 17, 2017 · 1 comment

Comments

@bretonics
Copy link
Owner

bretonics commented Jan 17, 2017

When a crispr sequence has no matching sequences in search sequence, then nothing is stored in hash. This causes an "undefined value" error when sorting unique identities for each match.

  • The integrity of the 'targets' hash must be checked to include all crispr occurrences (keys)
  • If no matches found for crispr, then store appropriate info in attributes
bretonics pushed a commit that referenced this issue Jan 19, 2017
Changing method calling BLAST+ using BioPerl's Bio::Tools::Run::StandAloneBlastPlus
@bretonics
Copy link
Owner Author

bretonics commented Jan 19, 2017

When using BLAST+, it will be better to use BioPerl's module's StandAloneBlastPlus to call blast and SearchIO to parse the results.

This in theory should allow identification of query sequences producing "no hits", setting values accordingly, and writes XML files as default.

bretonics pushed a commit that referenced this issue Jan 27, 2017
bretonics pushed a commit that referenced this issue Jan 27, 2017
bretonics pushed a commit that referenced this issue Jan 27, 2017
bretonics pushed a commit that referenced this issue Jan 30, 2017
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

1 participant