Quantcast
Channel: Post Feed
Viewing all articles
Browse latest Browse all 41826

Problem with blastp when blasting against custom made database

$
0
0
Hello, 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: >TRA_G345655|locus=scaffold45|99209|101155|-|translate_table|standard 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 or with: blastp -query query.fa -db DB/database.fa -out proteins_blastp_1e-30_table.txt -evalue 1e-30 -outfmt 6 or with: 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 f ...

Viewing all articles
Browse latest Browse all 41826

Trending Articles