[Bioperl-l] Genbank code from Blast results

Marcelo Iwata marcelo011982 at gmail.com
Tue Aug 18 18:34:17 UTC 2009


hi all..
I was doing a script that take some information of the results of blastn
files.
Everythig was ok, but i have some dificult to pic the Genbank code number
(the 'gb' below).
I tried

$obj->each_accession_number
$hit->name

And some variation of this.



------------------------------
>gnl|UG|Gma#S23062791 gmrtDrNS01_07-B_M13R_E11_087.s1 Water stressed 5h
segment 1 gmrtDrNS01
              Glycine max cDNA 3', mRNA sequence /clone_end=3'
              /gb=CX702616 /gi=58015874 /ug=Gma.18455 /len=853
          Length = 853

 Score = 1336 bits (674), Expect = 0.0
 Identities = 793/832 (95%), Gaps = 8/832 (0%)
 Strand = Plus / Minus


Query: 294858 aaattaacaatgagactccagagtatgtgaggtcctttgaatttgatagcaaattgatgt
294917
              |||||||||||| |||||| |||||||||||||||||  ||||||||||||||||||||
Sbjct: 853    aaattaacaatgtgactcccgagtatgtgaggtccttgaaatttgatagcaaattgatgc
794
----------------------------------------


But, i still don't get it.

thank you
with regards
Miwata



More information about the Bioperl-l mailing list