[Bioperl-l] Format of BlastHit '_description' key

Paul G Cantalupo lupey+ at pitt.edu
Fri Apr 29 13:30:54 EDT 2005


Hello,

Why does the implementation of the '_description' key of
Bio::Search::Hit::BlastHit join all the descriptions into one line? For
example, the description lines (there are 4 total) of this BlastHit


>ref|NP_040915.1| DNA polymerase [Human adenovirus A]
 emb|CAA51882.1| DNA polymerase [Human adenovirus type 12]
 pir||DJAD12 DNA-directed DNA polymerase (EC 2.7.7.7) - human adenovirus
12
 sp|P06538|DPOL_ADE12 DNA polymerase
          Length = 1061

 Score = 2172 bits (5627), Expect = 0.0
 Identities = 1049/1061 (98%), Positives = 1049/1061 (98%)



are joined into one scalar (as displayed by the Perl debugger)


'_description' => 'DNA polymerase [Human adenovirus A] emb|CAA51882.1| DNA polymerase [Human adenovirus type 12] pir||DJAD12 DNA-directed DNA polymerase (EC 2.7.7.7) - human adenovirus 12 sp|P06538|DPOL_ADE12 DNA polymerase'


The problem is that I need to be able to parse the description lines for
the one that I want. Is it possible to reconstruct the original
description lines as outputted by BLAST?

Thank you,

Paul

Paul Cantalupo
Research Specialist/Systems Programmer
559 Crawford Hall
Department of Biological Sciences
University of Pittsburgh
Pittsburgh, PA 15260
Work: 412-624-4687
Fax: 412-624-4759

Ask me about Toastmasters: www.toastmasters.org
Midday Club Treasurer


More information about the Bioperl-l mailing list