cablastp

The database is Uniref-90 fasta file compressed using compress command. The folder essentially contains following files : coarse.fasta coarse.fasta.index coarse.links coarse.links.index coarse.seeds compressed compressed.index params The database is in another user's account (though the complete path to it given to the progam) Tried running cablastp using the following command cablastp-search /path/capblastp-UNIREF90 /path/stn1_ALL_prodigal.prot --blast-args -max_target_seqs 10 -out stn1_ALL_prodigal.blast -num_threads 6 -outfmt '6 qseqid sseqid pident qcovs length mismatch qstart qend sstart send evalue bitscore stitle' gs -max_target_seqs 10 -num_threads 6 -out stn1_ALL_prodigal.blast -outfmt '6 qseqid sseqid pident qcovs length mismatch qstart qend sstart send evalue bitscore stitle' Opening database in path/capblastp-UNIREF90... Opening compressed database... Done opening compressed database. Opening coarse database... Done opening coarse database. Done opening database in /path/capblastp-UNIREF90 Blasting query on coarse database... blastp -db path/capblastp-UNIREF90/blastdb-coarse -num_threads 24 -outfmt 5 -dbsize 4923603367 Error blasting coarse database: Error running 'blastp -db path/blastdb-coarse -num_threads 24 -outfmt 5 -dbsize 4923603367': 'exit status 2'. stderr: BLAST Database error: No alias or index file found for protein database [path/blastdb-coarse] in search path [path_to_input_file:/home/db/blastdb:] Is this an issue with the compressing of the database ? it seems to open the compressed database just fine. The path given in stderr after 'in search path' is different from what I gave for the database.

Here is another one for ya. Not sure if the real error is "signal: bus error (core dumped)" or "stderr: num_threads' is currently ignored when 'subject' is specified." Blasting query on coarse database... blastp -db /blast_databases/cablastp_database/capblastp-cogs_db_aa/blastdb-coarse -num_threads 9 -outfmt 5 -dbsize 47558392 Decompressing blast hits... Blasting query on fine database... blastp -subject /tmp/cablastp-fine-fasta604138705 -dbsize 47558392 -num_threads 9 -max_target_seqs 10 -out /blastoutput/stn_ALL_prodigal.prot_2.blast -outfmt 6 qseqid sseqid pident qcovs length mismatch qstart qend sstart send evalue bitscore Error blasting fine database: Error running 'blastp -subject /tmp/cablastp-fine-fasta604138705 -dbsize 47558392 -num_threads 9 -max_target_seqs 10 -out /blastoutput/stn_ALL_prodigal.prot_2.blast -outfmt 6 qseqid sseqid pident qcovs length mismatch qstart qend sstart send evalue bitscore': 'signal: bus error (core dumped)'. stderr: 'num_threads' is currently ignored when 'subject' is specified.

It should be set to the value of the switch. We could parse the arguments given after , but I'd rather stick to our simple solution: don't touch them.