$genes{$gene[NAME]}{database}{$gene[DBNAME]}++;
$genes{$gene[NAME]}{hits}++;
$genes{$gene[NAME]}{terms}{$gene[KEYWORD]}++;
+ $genes{$gene[NAME]}{terms}{$gene[KEYWORD].'['.$gene[DBNAME].']'}++;
add_unique_parts($genes{$gene[NAME]},'refseq',$gene[REFSEQ]);
add_if_better($genes{$gene[NAME]},'description',$gene[DESCRIPTION]);
add_if_better($genes{$gene[NAME]},'location',$gene[LOCATION]);
print join(',',map {qq("$_")} @csv_fields),qq(\n);
for my $gene (keys %genes) {
- $genes{$gene}{rzscore} = scalar keys %{$genes{$gene}{terms}};
+ $genes{$gene}{rzscore} = scalar grep {$_ !~ /\[/} keys %{$genes{$gene}{terms}};
next if $genes{$gene}{rzscore} == 1 and exists $genes{$gene}{terms}{antigen};
$genes{$gene}{rzscore} -= 1 if exists $genes{$gene}{terms}{antigen};
print STDOUT join (',',