[Bioperl-l] using HMMER

fayroz fayroz_farouk at yahoo.com
Sun Aug 8 04:24:31 EDT 2010

i need your help, i am a new perl user and want to use bioperl modules to run 
HMMER program ( HMMsearch) i have" model.hmm" and a "fasta file" to see which of 
them are similar with the model
i write this code but there is a problems

#!/usr/local/bin/perl W
use Bio::AlignIO;
use Bio::SearchIO;
use Bio::SeqIO ;
use Bio::Tools::Run::Hmmer;

# run hmmsearch (similar for hmmpfam)
my $factory = Bio::Tools::Run::Hmmer->new(-hmm => 'h6_avian.hmm',-informat => 
my $seq = Bio::SeqIO->new('-file'=> "one_seq.fa", '-format'=>'Fasta');

# Pass the factory a Bio::Seq object or a file name, returns a Bio::SearchIO
my $searchio = $factory->hmmsearch($seq);

while (my $result = $searchio->next_result){
while(my $hit = $result->next_hit){
while (my $hsp = $hit->next_hsp){
print join("\t", ( $result->query_name,
)), "\n";

MSG: Unknown kind of input 'Bio::SeqIO::fasta=HASH(0x329a504)'
STACK Bio::Tools::Run::Hmmer::_setinput 
STACK Bio::Tools::Run::Hmmer::hmmsearch 
 STACK toplevel test_bioperl.pl:12
thank you



More information about the Bioperl-l mailing list