[BioRuby] help wtih bioruby wrapper for emboss tools...

Toshiaki Katayama ktym at hgc.jp
Wed Jan 9 02:04:23 UTC 2008


Hi,

On 2008/01/09, at 10:56, Toshiaki Katayama wrote:
> Real problem here is that the EBI server never finishes this work...
> 5 hours later, the process is still (reported to be) running
> thus I can't obtain the result yet.

This was due to the change of the entry name in Swiss-Prot as described at

  http://www.expasy.org/cgi-bin/idtracker?id=ALK1_HUMAN

> The ID: ALK1_HUMAN is no longer valid.
>
> ALK1_HUMAN (associated with primary accession number: P03973 from release 9.0) was renamed to SLPI_HUMAN in release 52.0.
>
> ALK1_HUMAN was associated with 1 accession number: P07757 from release 12.0.
>
> P07757 is now a secondary accession number of SLPI_HUMAN (primary accession number: P03973).

So, when I specify the USAs as 'uniprot:slpi_human' and 'uniprot:slpi_mouse'
instead of 'uniprot:alk1_human' and 'uniprot:alk1_mouse',
I could obtain the proper result.

Toshiaki


ktym at toge9064:~% bioruby
Loading config (/Users/ktym/.bioruby/shell/session/config) ... done
Loading object (/Users/ktym/.bioruby/shell/session/object) ... done
Loading history (/Users/ktym/.bioruby/shell/session/history) ... done

. . . B i o R u b y   i n   t h e   s h e l l . . .

  Version : BioRuby 1.2.1 / Ruby 1.8.6

bioruby> emboss = Bio::EBI::SOAP::Emboss.new
ignored attr: {}nillable
  ==> #<Bio::EBI::SOAP::Emboss:0x1974d38 @driver=#<SOAP::RPC::Driver:#<SOAP::RPC::Proxy:http://www.ebi.ac.uk/cgi-bin/webservices/WSEmboss>>, @log=nil, @wsdl="http://www.ebi.ac.uk/Tools/webservices/wsdl/WSEmboss.wsdl">
bioruby> hash = {
bioruby+ "tool" => "water",
bioruby+ "asequence" => "uniprot:slpi_human",
bioruby+ "bsequence" => "uniprot:slpi_mouse",
bioruby+ "email" => "ebisoap at example.org"
}
  ==> {"asequence"=>"uniprot:slpi_human", "tool"=>"water", "bsequence"=>"uniprot
:slpi_mouse", "email"=>"ebisoap at example.org"}
bioruby> poll = emboss.run(hash, [])
  ==> "emboss-20080109-01550394"
bioruby> emboss.checkStatus(poll)
  ==> "DONE"
bioruby> base = emboss.poll(poll, "tooloutput")
  ==> "IyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIwojIFByb2dyYW06ICB3YXRlcgojIFJ1bmRhdGU6ICBXZWQgSmFuIDA5IDAxOjU1OjA2IDIwMDgKIyBBbGlnbl9mb3JtYXQ6IHNyc3BhaXIKIyBSZXBvcnRfZmlsZTogL2ViaS9leHRzZXJ2L29sZC13b3JrL2VtYm9zcy0yMDA4MDEwOS0wMTU1MDM5NC5vdXRwdXQKIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIwoKIz09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PQojCiMgQWxpZ25lZF9zZXF1ZW5jZXM6IDIKIyAxOiBTTFBJX0hVTUFOCiMgMjogU0xQSV9NT1VTRQojIE1hdHJpeDogRUJMT1NVTTYyCiMgR2FwX3BlbmFsdHk6IDEwLjAKIyBFeHRlbmRfcGVuYWx0eTogMC41CiMKIyBMZW5ndGg6IDEzMgojIElkZW50aXR5OiAgICAgIDc5LzEzMiAoNTkuOCUpCiMgU2ltaWxhcml0eTogICAgOTYvMTMyICg3Mi43JSkKIyBHYXBzOiAgICAgICAgICAgMy8xMzIgKCAyLjMlKQojIFNjb3JlOiA0NjQuNQojIAojCiM9PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT0KClNMUElfSFVNQU4gICAgICAgICAxIE1LU1NHTEZQRkxWTExBTEdUTEFQV0FWRUdTR0stLVNGS0FHVkNQUEtLU0FRQ0xSWUtLICAgICA0OAogICAgICAgICAgICAgICAgICAgICB8fHwufHwufHwufHx8fHx8Lnx8fHwufHx8IHx8ICA6LnwufC58fC58fC58fHx8Oi46fApTTFBJX01PVVNFICAgICAgICAgMSBNS1NDR0xMUEZUVkxMQUxHSUxBUFdUVkVHLUdLTkRBSUtJR0FDUEFLS1BBUUNMS0xFSyAgICAgNDkKClNMUElfSFVNQU4gICAgICAgIDQ5IFBFQ1FTRFdRQ1BHS0tSQ0NQRFRDR0lLQ0xEUFZEVFBOUFRSUktQR0tDUFZUWUdRQ0xNICAgICA5OAogICAgICAgICAgICAgICAgICAgICB8Onw6Onx8Onx8fHw6fHx8LnwufHwufHw6Onx8Li4uLnwuLnx8fHw6fC4ufC4uOnw6fApTTFBJX01PVVNFICAgICAgICA1MCBQUUNSVERXRUNQR0tRUkNDUURBQ0dTS0NWTlBWUElSS1BWV1JLUEdSQ1ZLVFFBUkNNTSAgICAgOTkKClNMUElfSFVNQU4gICAgICAgIDk5IExOUFBORkNFTURHUUNLUkRMS0NDTUdNQ0dLU0NWU1BWICAgIDEzMAogICAgICAgICAgICAgICAgICAgICB8fHx8fC58Oi58fHx8Li4uLnx8fC58Onx8fC58Oi58OgpTTFBJX01PVVNFICAgICAgIDEwMCBMTlBQTlZDUVJER1FDREdLWUtDQ0VHSUNHS1ZDTFBQTSAgICAxMzEKCgojLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCiMtLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0K"
bioruby> require 'base64'
  ==> true
bioruby> puts Base64.decode64(base)
########################################
# Program:  water
# Rundate:  Wed Jan 09 01:55:06 2008
# Align_format: srspair
# Report_file: /ebi/extserv/old-work/emboss-20080109-01550394.output
########################################

#=======================================
#
# Aligned_sequences: 2
# 1: SLPI_HUMAN
# 2: SLPI_MOUSE
# Matrix: EBLOSUM62
# Gap_penalty: 10.0
# Extend_penalty: 0.5
#
# Length: 132
# Identity:      79/132 (59.8%)
# Similarity:    96/132 (72.7%)
# Gaps:           3/132 ( 2.3%)
# Score: 464.5
# 
#
#=======================================

SLPI_HUMAN         1 MKSSGLFPFLVLLALGTLAPWAVEGSGK--SFKAGVCPPKKSAQCLRYKK     48
                     |||.||.||.||||||.||||.||| ||  :.|.|.||.||.||||:.:|
SLPI_MOUSE         1 MKSCGLLPFTVLLALGILAPWTVEG-GKNDAIKIGACPAKKPAQCLKLEK     49

SLPI_HUMAN        49 PECQSDWQCPGKKRCCPDTCGIKCLDPVDTPNPTRRKPGKCPVTYGQCLM     98
                     |:|::||:||||:|||.|.||.||::||....|..||||:|..|..:|:|
SLPI_MOUSE        50 PQCRTDWECPGKQRCCQDACGSKCVNPVPIRKPVWRKPGRCVKTQARCMM     99

SLPI_HUMAN        99 LNPPNFCEMDGQCKRDLKCCMGMCGKSCVSPV    130
                     |||||.|:.||||....|||.|:|||.|:.|:
SLPI_MOUSE       100 LNPPNVCQRDGQCDGKYKCCEGICGKVCLPPM    131


#---------------------------------------
#---------------------------------------
  ==> nil






More information about the BioRuby mailing list