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

Using blast can find results, but diamond has no results #839

Open
ggra-billionare opened this issue Nov 13, 2024 · 2 comments
Open

Using blast can find results, but diamond has no results #839

ggra-billionare opened this issue Nov 13, 2024 · 2 comments

Comments

@ggra-billionare
Copy link

ggra-billionare commented Nov 13, 2024

Hi!
I need to compare two protein sequences. NCBI's blastp has the comparison results, but Diamond can never find the results. Is it that my parameters are set incorrectly?
The two sequences that need to be compared are as follows (the two sequences are only 83 and 42 in length, respectively):

>Putative_mercuric_reductase
MSAIFVPGNAVADKVILKIEGMTUAAUPLIIRKALEGLEGVEKASVSYSKGKGEVVFDPAKVSEKDIVDQVDRIGFRAKVIEE
>JALLOF010000030.1_Nitrospinae_bacterium_AH_259_B05_G02_I21_NODE_30_length_20278_cov_88.0098
MSRREMITTLVLGLLFLGMFVVPAVAELQTVTLTIEGMVUGA

The NCBI blastp result is as follows:
image

The command I used is as follows:
diamond blastp --db ./diamond/db -q ./query.txt --out ./diamond/result --outfmt 6 stitle qseqid sseqid pident length mismatch gapopen qstart qend sstart send sseq_gapped qseq_gapped evalue bitscore --threads 120 --masking 0 --ultra-sensitive --id 20 -e 1e-2

Thank you very much for your reply!

@bbuchfink
Copy link
Owner

DIAMOND is not configured to find very short hits by default. I shared some tips how to do this here: #832

@ggra-billionare
Copy link
Author

ggra-billionare commented Nov 15, 2024

Thank you for your reply!!
I try this issue that you mentioned parameters: #607, and it can work.

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