BEGIN{
    f="/data/results/reference/mmu/mm9/mRNA-stranded/Mus_musculus.NCBIM37.64-toMM9.headers";
    while (getline<f){
	en=$1;split($2,a,"=");split($3,b,"=");
	if (length(a[2])==0){print "ERR no gene "$0;}
	if (length(b[2])==0){print "ERR no CDS "$0;}
	split(b[2],c,"-");
	gi[en]=a[2];cdss[en]=c[1];cdse[en]=c[2];
	ok[en]=1;
    }
    close(f);
    FS="\t";
    f="/data/results/reference/mmu/Mus_musculus.NCBIM37.64-toMM9.transcript-gene_name.txt";
    while (getline<f){
	en=$1;gn[en]=$3;bt[en]=$4;
    }
    close(f);
}{
    en=$1;
    if (ok[en]==1){
	s=$2;e=$3;ncl++;iu5=0;iu3=0;icds=0;hit[en]++;rt="CDS";
	if (e<cdss[en]) {u5++;iu5=1;nu5[en]++;rt="5UTR";}
	if (s>cdse[en]) {u3++;iu3=1;nu3[en]++;rt="3UTR";}
	if ( (s>=cdss[en])&&(e<=cdse[en]) ) {cds++;icds=1;ncds[en]++;}
	print gi[en]"\t"cdss[en]"\t"cdse[en]"\t"$0;
    } else{
	rt=bt[en];
	print $0,"NA",rt,en,gn[en];
#	print "ERR unk ENS "$0;
    }
}END{
 }