I have been trying the last days to make my database for blast work but I really don't know how to figure out this problem.
Begining from the start of the problem, I want to make a blastp against a protein database that I have created.
Here is an example of my sequence identifiers and how they are formatted:
I try to make the blast database with the following command:
makeblastdb -in database.fa -title database -dbtype prot -out protdb -parse_seqids
It does not work and it gives me the following error:
No volumes were created because no sequences were found.
BLAST Database creation error: Defline lacks a proper ID around line 1
When I try it without '-parse_seqids', it works fine and it generated the three files with pin, psq and phr.
I put these 3 files with my database.fa fasta file in one folder called DB and then I try to run blastp with the following command:
blastp -query query.fa -db DB -out proteins_blastp_1e-30_table.txt -evalue 1e-30 -outfmt 6
blastp -query query.fa -db DB/database.fa -out proteins_blastp_1e-30_table.txt -evalue 1e-30 -outfmt 6
blastp -query query.fa -db database -out proteins_blastp_1e-30_table.txt -evalue 1e-30 -outfmt 6
But it gives me the following error:
BLAST Database error: No alias or index file found for protein database [database] in search path
How can I solve this problem or where is the origin of this problem?