use strict;
use warnings;
use Acme::Tools;
use Text::Table;

my %data;
my %freq;
my $ref_filelist = $ARGV[0];
open(BLASTFILE, $ref_filelist )
   or die "Could not open Reference filelist...($!)";
   
while (<BLASTFILE>) {
	chomp;
	
    my ($accession,$organism, @vals) = (split /\t/)[1..5];
    my $organismcut = substr( $organism, 0,75 );
    my $tot =$accession . $organismcut;
    #print "$tot\n";
    
    $freq{ $tot }++;
    
    my $col = 4;
    for my $val (@vals) {
        push @{ $data{$tot} {$col} }, $val;
        $col++;
    }
}

my @headers =  qw/ Organism Freq Median_Eval Med_Contig_Length Med_Mapped_Length /;

my $tb = Text::Table->new( map {title => $_}, @headers);

for my $test (sort {$freq{ $b } <=> $freq{ $a }} keys %freq) {
    my @row = ($test, $freq{ $test });
    for my $col (sort keys %{ $data{$test} }) {
        push @row, median(@{ $data{$test}{$col} });
        
    } 
    $tb->load( [@row] );
    #print "@row\n"; 
    
}
print $tb;