#!/usr/bin/perl

use strict;


open(LENGTH,"<estlength");

my %long_est;

while(<LENGTH>){
    my @entries = split;
    $long_est{$entries[0]} = $entries[1];
}



open(IDENT,"<identitats3");

my %percid;

while(<IDENT>){
    my @entries2 = split;
    $percid{$entries2[0]}{$entries2[1]} = $entries2[2];
}



open(GFF,"<hum_b_aln_length");

my %longitud;

while(<GFF>){
    my @entries = split;
    $entries[1] =~s/\;//;
    $longitud{$entries[0]}{$entries[1]} = $entries[2];

    my $coverage = sprintf"%.2f",(100*$entries[2]/$long_est{$entries[1]});



    print $entries[0],"\t",$entries[1],"\t",$percid{$entries[0]}{$entries[1]},"\t",$coverage,"\n";

}




