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

Ability to match parameters of BLAST as per the web app #58

Open
vkullu opened this issue Jan 21, 2023 · 2 comments
Open

Ability to match parameters of BLAST as per the web app #58

vkullu opened this issue Jan 21, 2023 · 2 comments
Labels
enhancement New feature or request

Comments

@vkullu
Copy link

vkullu commented Jan 21, 2023

Hello,
This request is to make the params of BLAST match exactly that of the web app (or an ability to pass them as xargs to the code)

This is an example of params sent by the code today:

   
Program blastn
Word size 28
Expect value 10
Hitlist size 100
Match/Mismatch scores 1,-2
Gapcosts 0,2.5
Low Complexity Filter Yes
Filter string L;m;
Genetic Code 1

and these are the params if BLAST is done on the web directly

Program blastn
Word size 11
Expect value 0.05
Hitlist size 100
Match/Mismatch scores 2,-3
Gapcosts 5,2
Low Complexity Filter Yes
Filter string L;m;
Genetic Code 1

Thank you very much

@lauraluebbert
Copy link
Member

Hi Vinayak, thanks for reaching out. gget currently matches the standard settings for megablast (on by default). It looks like the settings from the web you posted below are the standard settings for blastn without megablast. One exception is the expect value which was originally set to 10 as default (reference), but NCBI seems to have changed it to 0.05 (the expect value can already be modified by the user). I'll work on adding more arguments to adjust all of the options manually.

Could you please post your gget version, gget blast call, and how you returned the parameters?

@lauraluebbert lauraluebbert added the enhancement New feature or request label Jan 22, 2023
@vkullu
Copy link
Author

vkullu commented Jan 23, 2023

Thank you Laura for a quick response. I was able to recreate near identical params to online version following your prompt.
Thank you for such a useful piece of software!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
enhancement New feature or request
Projects
None yet
Development

No branches or pull requests

2 participants