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

low sensitibity using local search mode #68

Closed
yye88 opened this issue Nov 10, 2021 · 4 comments
Closed

low sensitibity using local search mode #68

yye88 opened this issue Nov 10, 2021 · 4 comments

Comments

@yye88
Copy link

yye88 commented Nov 10, 2021

Hello, thank you for previous helps. The search went great in NCBI BLAST API. Then I tried the local database, but found that the sensitivity is quite low. I am sure my local genome library has at least ten clusters containing both gene_a and gene_b, but the search returned me zero hit. I adjusted the –min_identity argument, but made no difference. Would you please tell me how to increase the sensitivity of the local search mode? Thank you.

@gamcil
Copy link
Owner

gamcil commented Nov 13, 2021

Hi, most of the time this is the result of the clustering arguments - could you try lowering the --min_hits and --unique arguments? By default, they are set to 3, so clusters containing less than 3 unique hits are discarded.

For the actual sequence search sensititivity, I'll be adding an argument for controlling the sensitivity of DIAMOND in the next update.

@yye88
Copy link
Author

yye88 commented Nov 13, 2021

Thank you very much for the suggestion and future update!

@gamcil
Copy link
Owner

gamcil commented Dec 13, 2021

Closing this issue as you can now adjust the DIAMOND sensitivity as of v1.3.12 (https://github.com/gamcil/cblaster/releases/tag/v1.3.12).

@gamcil gamcil closed this as completed Dec 13, 2021
@yye88
Copy link
Author

yye88 commented Dec 13, 2021 via email

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