[BioRuby] Blast Problem

Toshiaki Katayama ktym at hgc.jp
Wed Apr 11 13:51:05 UTC 2007


Hi,

To use GenomeNet server, you need to select the database name
which the server can accept (described in the bio/appl/blast.rb
or see the HTML source of http://blast.genome.jp/ form).

------------------------------------------------------------
program = 'blastp'
database = 'swissprot'     # <- here
options = '-e 0.0001'
server = 'genomenet'

serv = Bio::Blast.remote(program, database, options, server)
report = serv.query(sequence_text)

p report
------------------------------------------------------------

Toshiaki

On 2007/04/11, at 8:47, george wrote:

> Hi,
> am trying to run a remote blast job as follows:
>
> require 'bio'
> sequence_text='MRKLYCVLLLSAFEFTYMINFGRGQNYWEHPYQNSDVYRPINEHREHPKEYEYPLHQEHT
> YQQEDSGEDENTLQHAYPIDHEGAEPAPQEQNLFSSIEIVERSNYMGNPWTEYMAKYDIE
> EVHGSGIRVDLGEDAEVAGTQYRLPSGKCPVFGKGIIIENSNTTFLTPVATGNQYLKDGG
> FAFPPTEPLMSPMTLDEMRHFYKDNKYVKNLDELTLCSRHAGNMIPDNDKNSNYKYPAVY
> DDKDKKCHILYIAAQENNGPRYCNKDESKRNSMFCFRPAKDISFQNYTYLSKNVVDNWEK
> VCPRKNLQNAKFGLWVDGNCEDIPHVNEFPAIDLFECNKLVFELSASDQPKQYEQHLTDY
> EKIKEGFKNKNASMIKSAFLPTGAFKADRYKSHGKGYNWGNYNTETQKCEIFNVKPTCLI
> NNSSYIATTALSHPIEVENNFPCSLYKDEIMKEIERESKRIKLNDNDDEGNKKIIAPRIF
> ISDDKDSLKCPCDPEMVSNSTCRFFVCKCVERRAEVTSNNEVVVKEEYKDEYADIPEHKP
> TYDKMKIIIASSAAVAVLATILMVYLYKRKGNAEKYDKMDEPQDYGKSNSRNDEMLDPEA
> SFWGEEKRASHTTPVLMEKPYY'
>
>
> program = 'blastp'
>   database = 'SwissProt'
>   options = '-e 0.0001'
>   serv = Bio::Blast.new(program, database, options)
>
>   server = 'genomenet'
>   genomenet = Bio::Blast.remote(program, database, options, server)
>
>   report = serv.query(sequence_text)
>
> But i get this:
>
>> ruby blast.rb
> [NULL_Caption] WARNING: Could not find index files for database SwissProt
> [NULL_Caption] WARNING: Could not find index files for database SwissProt
> [NULL_Caption] WARNING: Could not find index files for database SwissProt
> c:/ruby/lib/ruby/gems/1.8/gems/bio-1.0.0/lib/bio/appl/blast/format8.rb:77:in 
> `tab_parse_hsp': undefined method `strip' for nil:NilClass (NoMethodError)
>     from 
> c:/ruby/lib/ruby/gems/1.8/gems/bio-1.0.0/lib/bio/appl/blast/format8.rb:40:in 
> `tab_parse'
>     from 
> c:/ruby/lib/ruby/gems/1.8/gems/bio-1.0.0/lib/bio/appl/blast/format8.rb:38:in 
> `each'
>     from 
> c:/ruby/lib/ruby/gems/1.8/gems/bio-1.0.0/lib/bio/appl/blast/format8.rb:38:in 
> `tab_parse'
>     from 
> c:/ruby/lib/ruby/gems/1.8/gems/bio-1.0.0/lib/bio/appl/blast/report.rb:54:in 
> `auto_parse'
>     from 
> c:/ruby/lib/ruby/gems/1.8/gems/bio-1.0.0/lib/bio/appl/blast/report.rb:70:in 
> `initialize'
>     from 
> c:/ruby/lib/ruby/gems/1.8/gems/bio-1.0.0/lib/bio/appl/blast.rb:202:in `new'
>     from 
> c:/ruby/lib/ruby/gems/1.8/gems/bio-1.0.0/lib/bio/appl/blast.rb:202:in 
> `parse_result'
>     from 
> c:/ruby/lib/ruby/gems/1.8/gems/bio-1.0.0/lib/bio/appl/blast.rb:216:in 
> `exec_local'
>     from 
> c:/ruby/lib/ruby/gems/1.8/gems/bio-1.0.0/lib/bio/appl/blast.rb:182:in `send'
>     from 
> c:/ruby/lib/ruby/gems/1.8/gems/bio-1.0.0/lib/bio/appl/blast.rb:182:in 
> `query'
>     from blast.rb:23
>> Exit code: 1
>
>
> What could be the reason?
>
> Thanks
>
> _______________________________________________
> BioRuby mailing list
> BioRuby at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioruby




More information about the BioRuby mailing list