BEGIN{
    FS=OFS=" ";
    f="/data/images/proton/DKlab/orsalia/stringtie/edgeRtranscripts_wtVSko.withGID.csv";
    getline<f;h=$0;
    while(getline<f){nr++;
	id=$3;all[id]=$0;ok[id]=1;
	gid=$2;Aqval[id]=$NF;
	i=++ngid[gid];tgrp[gid,i]=id;
	if ($1!="NA"){gn[gid]=$1;gnok[gid]=1;}
    }
    close(f);
    print "ref transcripts "nr >>"/dev/stderr"
}{
    f1=$1; f=substr(f1,1,length(f1)-4);
    he[f]=1;
    getline<f1;n=0;
    while(getline<f1){
	id=$1;n++;
	partfc[f,id]=$2;
	partcpm[f,id]=$3;
	partpval[f,id]=$4;
	partqval[f,id]=$5;
	okpart[f,id]=1;
    }
    close(f1);    
    print "transcripts in "f" : "n >>"/dev/stderr"
}END{
    printf("%s",h);
    for (i in he){
	printf(" %s_fc %s_cpm %s_pval %s_qval",i,i,i,i);
    }
    printf("\n");
    for (gid in ngid){
	ng=ngid[gid];
	if (ng>1){
	    gok=0;
	    for (i=1;i<=ng;i++){
		id=tgrp[gid,i];
		if (1.0*Aqval[id]<=0.05){
		    j=++gok;Lgok=i;
#		    print "OK1"gid" "Aqval[id]" "id" "Afc[id]" "Apval[id];
		} #one sig DE iso
	    }
	    if (gok>0){
		if (gok==1){
		    gok2=0;
		    for (i=1;i<=ng;i++){
			if (i!=Lgok){
			    id=tgrp[gid,i];
			    if (1.0*Aqval[id]<=0.05){gok2=1;}
			}
		    }  
		} else{
		    gok2=1;				
		}
		if (gok2==1){
		    for (j=1;j<=ng;j++){
			id=tgrp[gid,j];
#				gid2=Agid[id];
#	if (gnok[gid2]==1){
#	    printf("%s",gn[gid2]);
#	}else{
#	    printf("%s",Agn[id]);
#	}
#	printf(",%s,%s,%s,%s",gid,Aft[id],id,Afc[id],Apval[id],Aqval[id]);
			printf("%s",all[id]);
	for (i in he){
	    if (okpart[i,id]==1){
		printf(" %s %s %s %s",partfc[i,id],partcpm[i,id],partpval[i,id],partqval[i,id]);
	    }else{
		printf(" NA NA NA NA");
	    }
	}
	printf("\n");

		    }
		}
	    }
	}
    }
}