blastx generates empty output file
I am trying to run blastx (version 2.2.28) with the following command line:
blastx -db nr -outfmt 5 -evalue 1e-3 -word_size 3 -show_gis -num_alignments 1 \
-max_hsps 20 -num_threads 5 -out transcripts_refonly.xml -query transcripts_refonly.fa
So when blastx is running, I can see the output file that grows in size, however when blastx finishes, the output file then is empty (0kb). This happens when I run it with big files like transcprit fasta file (from a RNA-seq project) of ~30Mb (21000 sequences), but if I run just a smaller file (20 sequences) it works good.
I tried also with .txt ouput format and it happens the same, and also run it into the cluster and different servers and still not working.
Does someone knows what I am doing wrong?
Thank you
1 answer
Log in to answer this question.
More posts like this
-
Issues while running blastx
written by Foad •
-
BLAST Database error: No alias or index file found for nucleotide database
written by langziv •
-
Blast command use
written by harry •
-
Puzzling NCBI BLAST+ output after transcriptome assembly
written by ella •
-
Blast run fails: memory exceeds the limit
written by langziv •
-
Blastx against nr for plants
written by m.s_tabatabaei •
-
tBlastn taking an absurd amount of time
written by Rogerio Ribeiro •
-
Blastn .XML to FASTA
written by baurumon •
-
how come that blastx with XML output format was done at reasonable time?
written by seta
-
Help For Blastx, No Output Untill The Program Finished
written by zhangmin8712 •
Sounds like you are running out of space (or hitting some limit on storage quota). Have you looked into those possibilities?
It does not seem to be a space problem... I have about 853TiB in my workspace in the cluster. When I run it in the server I had about 30GB (with no limit), do you think it will need that amount?
If you are sure that there are no storage quota then this is not making sense. And this happens reproducibly? If you are starting multiple jobs you are using unique outfile names for each job (example above seems generic but it may just be an example)?
yes, my outfile names are unique for each job