[BioRuby] Remote Blast
Anthony Underwood
email2ants at gmail.com
Fri Jul 4 08:13:45 EDT 2008
Hi all, is remote blast broken for bioruby?
I have some code
fasta_sequences.each do |seq_obj|
blast_report = remote_blast(seq_obj.seq, 'blastn', 'nr-nt')
..........
end
this calls the method
def remote_blast(seq_obj, program, db = 'nr-nt', options = '', server
= 'genomenet')
# create BLAST factory object
factory = Bio::Blast.remote(program, db, '-m 8' + options, server)
report = factory.query(seq_obj)
end
This fails to return anything - just times out
If I include the -V option to limit the number of hits returned
blast_report = remote_blast(seq_obj.seq, 'blastn', 'nr-nt', "-V 10")
I get no hits , or a message "RuntimeError: cannot understand
response" which I think is due to the "holding page" that occurs when
a job is running
Does anybody know a fix for this? Should I start to try and get a
ncbi remote blast working. Is there anybody who would like to help
with this?
Thanks
Anthony
More information about the BioRuby
mailing list