BEGIN{
    f="/data/images/proton/DKlab/run392/run402/featurecount-log.txt";
    while (getline<f){
	tr[$1]=$2;
#	print $1,$2;
    }
    close(f)
    FS=OFS="\t";
    getline;
    getline;
    print $0;
    ns=NF;
    for (i=7;i<=ns;i++){
	n=split($i,a,"/");
	sampID[i]=substr(a[n],1,length(a[n])-4);
#	print sampID[i],tr[sampID[i]]	
	tr2[i]=tr[sampID[i]];
	if (tr2[i]==0) {print "ERR : "sampID[i]" not found in "f}
    }
}{
    fl[NR]=$6;ns=NF;
    for (i=1;i<7;i++){
	r[NR,i]=$i;
    }
    for (i=7;i<=ns;i++){
	x=$i;
	nr[i]+=x;
	r[NR,i]=x;
    }
}END{
    nl=NR;
#    print "Sample","Total_mapped_reads"
#    ofn1[ctl1i]="RPKM_"ctl1"_vs_"ctl2;
    for (i=7;i<=ns;i++){
	rpm[i]=tr2[i]/1000000.0;
#	print sampID[i],nr[i] >>"
    }
    for (l=3;l<=nl;l++){
	printf("%s",r[l,1]);
	for (i=2;i<7;i++){
	    printf("\t%s",r[l,i]);
	}
	kb=fl[l]*0.001;
#	print r[l,1],kb
        for (i=7;i<=ns;i++){
	    rpkm=(r[l,i]/rpm[i])/kb;
	    if (rpkm==0){
		printf("\t0");
	    } else {
		printf("\t%f",rpkm);
	    }
	}
	printf("\n");
    }
 }