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

Diamond with database of small kmers #811

Open
raeshrode opened this issue May 28, 2024 · 0 comments
Open

Diamond with database of small kmers #811

raeshrode opened this issue May 28, 2024 · 0 comments

Comments

@raeshrode
Copy link

Hello,

I am trying to run diamond blastp. I built my own diamond database with 9-12mer amino acids strings and I'm comparing those to .faa files with ~200aa. I have tried running diamond blastp with several different parameters and each time it results in 0 pairwise alignments. Below are some of the variations I've tried:

diamond blastp -d database.dmnd -q prac.faa -o match_prac.tsv

diamond blastp -d database.dmnd --matrix PAM30 -q prac.faa -o match_prac.tsv

diamond blastp -d database.dmnd --matrix PAM30 --ultra-sensitive -q prac.faa -o match_prac.tsv

diamond blastp -d database.dmnd --matrix PAM30 --ultra-sensitive --evalue 100 -q prac.faa -o match_prac.tsv

What do you suggest I do to improve my alignments? Is there a reason a database with short reads is unable to match?

Thanks,
Rachel

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

1 participant