[Bioperl-l] parsing hmmpfam reports
    stephan.rosecker at ish.de 
    stephan.rosecker at ish.de
       
    Fri Apr 30 07:07:12 EDT 2004
    
    
  
Hi,
i'm looking for a method to get the family classification (family model)
in an hmmpfam report.Maybe some usefull stuff in
Bio::Search::Result::GenericResult ?
For example:
hmmpfam - search one or more sequences against HMM database
HMMER 2.3.2 (Oct 2003)
Copyright (C) 1992-2003 HHMI/Washington University School of Medicine
Freely distributed under the GNU General Public License (GPL)
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
HMM file:                 nst.hmms
Sequence file:            uniprot_sprot.dat
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query sequence: 128U_DROME
Accession:      P32234
Description:    GTP-binding protein 128UP.
Scores for sequence family classification (score includes all domains):
Model    Description                                    Score
E-value  N
-------- -----------                                    -----    -------
---
Drg                                                     639.4
1.6e-191   1
Parsed for domains:
Model    Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
-------- ------- ----- -----    ----- -----      -----  -------
Drg        1/1       4   367 ..     1   375 []   639.4 1.6e-191
.
.
.
What i want is the "Drg" string.
Regards,
stephan
    
    
More information about the Bioperl-l
mailing list