Hi All!
I would like to know how can I run BLAST locally in multiple files. Here is the command I am using for a single file:
bsub -n 16 blastp -query sequences_dsrb_sol_1_10.fasta -db database/nhr -outfmt '6 std salltitles' -num_threads 16 -max_target_seqs 1 -evalue 1e-30 -out faa_two_testing_s3_proteins_sol_1_10_taxonomy.txt
And here is the command I tried for multiple files:
bsub -n 16 ls *.fasta | parallel -a - blastp -query {} -db database/nhr -outfmt '6 std salltitles' -num_threads 16 -max_target_seqs 1 -evalue 1e-30 -out {}.txt
However, I get the following error:
/usr/bin/bash: blastp: command not found
Do you know any other methods to run BLAST in parallel for multiple fasta files?
Thank you!
2 answers
Are you using a job scheduler on a compute cluster? Based on bsub in your command, you must be. If you are on a cluster using parallel on top likely does not make sense since you may as well start individual jobs in parallel (pun intended) using the scheduler.
Since the error is about blastp not being found you first need to sort out $PATH or if your cluster uses modules system then load the module for blast first.
contact your system admin and ask them how your system is set up
So sounds like you are using modules for software management on your cluster? So relevant module would need to be added in the shell or inside your script.
If you have a job scheduler available then using that for independent jobs to run blast on the list of files makes more sense than using parallel.
You can use SEDA (https://www.sing-group.org/seda/). The "BLAST" operation (https://www.sing-group.org/seda/manual/operations.html#blast) allows to do this.
Log in to answer this question.
-
BLASTP error - too many positional arguments
written by Paula •
-
Blast "Argument "query". File is not accessible:" Error
written by logbio •
-
outfmt 6 BLAST column names
written by Paula •
-
BLAST ERROR Too many positional arguments (1)
written by Paula •
-
Blast in parallel with edible outfmt
written by K.Gee •
-
Core proteome using local Blast
written by zarodkip •
-
Blast locally on multiple fasta files with multiple database
written by fec2 •
-
Speed up BLASTp vs NCBI nr database
written by biotech
-
Blast Error: Tax database file not found
written by navillusol858 •
-
Orthomcl Blast Error: No Hits Are Being Saved
written by User 6777 •
