Commit activity
Last 12 weeks · 0 commits
Community health
2 of 6 standards met
42
✓README✓License○Contributing○Code of Conduct○Issue Template○PR Template
Recent PRs & issues
Quiet · 4 discussions · Last activity 11y ago
See all on GitHub →
Hi, I just realized that cablastp-extract appends characters to my protein sequnces: Input Output Is there a reason for that, or is that a bug?
While I was going over the NW code in biogo, I noticed that the lookup for backtracking is incorrect, the error is also in the cablastp backtracking code. The lookup into matrix is transposed between left and up - this is unlikely to make a difference in any normal matrix. I also corrected the approach to assessing which path to take on the backtrack - the inequality tests were a kludge. Please take a look at the current state of the NW backtracking code in the align package: http://code.google.com/p/biogo/source/browse/align/nw_letters.go#130
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.
Recent fixes
I tried blastp against downloaded compressed nr database. However got an error during the run.... something about not being able to find a file. The only file I downloaded was the Compressed NCBI NR dated Jun 30, 2013 (http://groups.csail.mit.edu/cb/cablastp/cablastp-nr20130630.tar.gz) Is there any way I can check if this is not an issue with the downloaded compressed database ? Opening database in /nv/hp10/nsarode3/blast_databases/cablastp_database/cablastp-nr20130630... Opening compressed database... Done opening compressed database. Opening coarse database... Done opening coarse database. Done opening database in /nv/hp10/nsarode3/blast_databases/cablastp_database/cablastp-nr20130630. Blasting query on coarse database... blastp -db /nv/hp10/nsarode3/blast_databases/cablastp_database/cablastp-nr20130630/blastdb-coarse -num_threads 10 -outfmt 5 -dbsize 9281362451 Error blasting coarse database: Error running 'blastp -db /nv/hp10/nsarode3/blast_databases/cablastp_database/cablastp-nr20130630/blastdb-coarse -num_threads 10 -outfmt 5 -dbsize 9281362451': 'exit status 2'. BLAST Database error: Could not find volume or alias file (cablastp-nr20130630/blastdb-coarse.00) referenced in alias file (/nv/pb4/bio-stewart/cablastp_database/cablastp-nr20130630/blastdb-coarse). Thanks, Neha
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.
Structured data for AI agents
Repository: BurntSushi/cablastp. Description: Performs BLAST on compressed proteomic data. Stars: 18, Forks: 5. Primary language: Go. Languages: Go (92.9%), Ruby (4.3%), Python (2.1%), Makefile (0.4%), Shell (0.3%). License: GPL-2.0. Open PRs: 0, open issues: 4. Last activity: 11y ago. Community health: 42%. Top contributors: BurntSushi, ndaniels.