[ create a new paste ] login | about

Link: http://codepad.org/P42y2ZB1    [ raw code | fork ]

Perl, pasted on Jun 1:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
#open up the hmmsearch result file
my $in = Bio::SearchIO->new(-format => 'hmmer3',-file => $argument);

while (my $result = $in->next_result ){
  # get a Bio::Search::Result::HMMERResult object
  print $result->query_name, " for HMM ", $result->hmm_name, "\n";
  
	while (my $hit = $result->next_hit ){
		my $score = $hit->raw_score();
  		my $evalue = $hit->significance();
		my $hitid = $hit->name();

 	}
}


Create a new paste based on this one


Comments: