BEGIN{FS="\t";OFS=",";
    print "Chromosome,Strand,ClusterStart,ClusterEnd,ClusterID,ClusterSequence,ReadCount,ModeLocation,ModeScore,ConversionLocationCount,ConversionEventCount,NonConversionEventCount,FilterType,TranscriptLocation,TranscriptID,GeneName"
}{
    n=split($7,a,"_");p=a[4];loc="CDS";
    if (p<=$2){loc="5UTR";}
    if (p>=$3){loc="3UTR";}
    o[a[6]]++;s+=a[6];nr++;
#    if ((a[6]>10)&&(a[8]>=0.975)&&(a[5]>2)&&(a[3]>10)) {
#    if ((a[6]>10)&&(a[8]>=0.99)&&(a[5]>2)) {
    if ((a[6]>10)&&(a[8]>=0.99)) {
#    if ((a[8]>=0.995)) {		
    #                 ci   seq  rc   ml   ms   clc  cec  ncc
    print $4,$9,$5,$6,a[1],a[2],a[3],a[4],a[8],a[5],a[6],a[7],"NA",loc,$4,$1
    }
}END{print (1.0*s)/nr" "nr>>"/dev/stderr"}