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

Guidance to search small fragments (8-mer) #368

Open
alexproteomics opened this issue Oct 15, 2024 · 2 comments
Open

Guidance to search small fragments (8-mer) #368

alexproteomics opened this issue Oct 15, 2024 · 2 comments

Comments

@alexproteomics
Copy link

Hi,
I'm new to foldseek and I wonder if you can you provide guidance on how I can obtain structure alignment of a small fragment of a protein. Basically, the protein has a 8-aa forming a beta-hairpin. I want to align this against the entire alphafold swissprot. I tried the default easy-search but nothing came back in my output file. I understand that a small fragment will align to too many sequences and maybe nothing is passing the minimum thresholds. I'm ok with having false positives since I have other means to filter the table downstream, but since I'm not getting anything back from my easy-search, I guess I need to change some parameters to allow the low confidence hits to be written to my output file.
Again, I would appreciate some guidance on how to address this.
thanks
Alex

@milot-mirdita
Copy link
Member

Please try --prefilter-mode 1. This is a different prefiltering algorithm that doesn't have a minimum length for matches. The k-mer prefilter cannot match anything shorter than ~13aa in its default settings.

Other parameters that you might need to change are the E-value threshold -e inf and minimum ungapped score --min-ungapped-score 5 (15 is default, maybe you'll need to run an even lower score).

@alexproteomics
Copy link
Author

@milot-mirdita Thank you very much for your help

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