[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