[Bioperl-l] hmmpfam

[Heikki Lehvaslaiho]

Joonas,

I am not familiar with hmmpfam, but is is part of the HMMERtoolkit, so it 
should be possible to use Bio::SearchIO modules for parsing the output. In 
addition, it is possible to run HMMER programs directly from BioPerl (the 
bioperl-run repositoty is needed):

See:

http://doc.bioperl.org/releases/bioperl-current/bioperl-run/Bio/Tools/Run/Hmmer.html
http://doc.bioperl.org/releases/bioperl-current/bioperl-live/Bio/SearchIO/hmmer.html

for code examples, and

SearchIO HOWTO
http://bio.perl.org/wiki/HOWTO%3ASearchIO

for an introduction to Bio::SearchIO.

    -Heikki

On Thursday 19 June 2008 00:52:47 Joonas Jämsen wrote:
> How to parse hmmpfam output to get gi/sequence of specified domain,
> score, E-value?
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l



-- 
______ _/      _/_____________________________________________________
      _/      _/
     _/  _/  _/  Heikki Lehvaslaiho    heikki at_sanbi _ac _za
    _/_/_/_/_/  Senior Scientist    skype: heikki_lehvaslaiho
   _/  _/  _/  SANBI, South African National Bioinformatics Institute
  _/  _/  _/  University of Western Cape, South Africa
     _/      Phone: +27 21 959 2096   FAX: +27 21 959 2512
___ _/_/_/_/_/________________________________________________________