BlastX Database: output search is empty
Dear all,
I am running a blastx search using a cluster. I uploaded the nr database as well as the taxon indices taxdb.btd and taxdb.bti into a local folder and ran a script as
blastx -db <path_to_nr_database/nr> -num_threads 30 -max_target_seqs 10 -outfmt '6 sskingdoms' -query <file>
I also exported the variable for the database
export BLASTDB="<path_to_nr_database/nr>"
but the output is empty. Is there something wrong in the syntax?
• 2,406 views
•
link
written
by
marongiu.luigi
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
BLASTX script doesn't print results and never ends
written by langziv •
-
Failing to run tblastn
written by langziv •
-
How to solve Warning: [blastx] Taxonomy name
written by marongiu.luigi
-
magicblast aganist local nr database
written by Assa Yeroslaviz
-
BlastX Database error: No alias or index file found for protein database
written by marongiu.luigi
-
error when running blastn
written by Nifaste •
-
Taxid will not function with other databases (including custom)
written by emilyc •
-
BLAST options error: tried to set option (130) and value (1), line (537)
written by oliver.bayfield •
-
Blast Error: Tax database file not found
written by navillusol858 •
-
Help For Blastx, No Output Untill The Program Finished
written by zhangmin8712 •
How about just
the same, it is still empty, as is in "N/A"
does that cmdline provide any feedback/errors/warnings ?
on the other hand if there are no hits, then with the
-outfmt 6you will indeed get an empty file (as it only prints hits info to the output in that case)no errors, but when I ran a check on my machine I got a hit, thus I am afraid something is not right when running in the cluster...
Is
$BLASTDBdirectory available on all cluster nodes? Sometimes parts of a cluster may not be able to access all storage.I am using
bsub ... -env "all, BLASTDB=<path>"andexport BLASTDB=<pathin the shell script but when I doecho $BLASTDBthe result is emtpy...If your
$BLASTDBvariable is not propagating correctly in your shell script then that is the problem. It is not clear how you submit these jobs so once you figure out why that is happening the problem should be fixed.clearly blastx does not find a hit in the taxon database. Is there a way to force blastx to find the taxon database more directly than with the environmental variable BLASTDB? Also, are taxdb.btd and taxdb.bti the correct databases for protein hits? Thanks
can you first check if you have the taxon part of the DB? I seem to remember you had to download that one separately, it was (is?) not part of the default blastdb release.
from the NCBI blast website:
The md5sum of the files in the link is different from those of the files I had, I guess it is only a different update, but I will try to see if there is a difference. Thank you