my $in_blast_tab=$ARGV[0];
open(IN,$in_blast_tab) or die "cannot open $in_blast_tab\n";

my $HoFlg1={}; my $HoFlg2={};
my %maxBits;
while (my $line=<IN>) {
        next if ($line =~ /^#/);
        next unless ($line =~ /\S/);
        chomp $line;
        my ($Query_id,$strand,$Subj_id,$Perc_iden,$align_len,$num_mm,$gap,$q_start,$q_end,$s_start,$s_end,$e_value,$bit_score)=split("\t",$line); # extra field of strand
        next if ($bit_score <60);
        if (!exists $maxBits{$Query_id}){
                $maxBits{$Query_id}=$bit_score;
        }
        elsif (exists $maxBits{$Query_id} && $bit_score > $maxBits{$Query_id}){
                $maxBits{$Query_id}=$bit_score;
        }
        my ($Flag1, $Flag2 ) = &Flag( $Subj_id, \%proph_prots, \%euk_prots, \%vir_prots );
        $HoFlg1->{$Query_id}->{$Flag1}->{$bit_score}++;
        $HoFlg2->{$Query_id}->{$Flag2}->{$bit_score}++;
        print join("\t",$Query_id,$Subj_id,"bit",$bit_score,"F1",$Flag1,"F2",$Flag2,"maxBits{Query}",$maxBits{$Query_id})."\n";

}