[Bioperl-l] Bio::SearchIO::hmmer hsp behaviour

Chris Fields cjfields at uiuc.edu
Wed Jun 28 22:16:38 UTC 2006


Arghhhh!  Made a mistake:

> my $result_count = 1;
> while ( my $result = $searchio->next_result() ) {
>   print "Result $result_count : ",$result->query_name,"\n";
>   print "Result models: ",$result->num_hits,"\n";
>   while (my $model = $result->next_hit) {
> 	print "\tModel : ",$model->name,"\n";
> 	print "\tSignif: ",$model->significance,"\n";
> 	while (my $domain = $model->next_hsp) {
> 	  print "\t\tDomain : ",$domain->name,"\n";
                              ^^^^^^^
Should be:                    $model

> 	  print "\t\tEvalue : ",$domain->evalue,"\n";
> 	}
>   }
>   $result_count++;
> }

My bad!

Chris




More information about the Bioperl-l mailing list