[ create a new paste ] login | about

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

Perl, pasted on Jun 1:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
#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();
		print $hitid;
}
}


Create a new paste based on this one


Comments: