BEGIN{FS=OFS="\t";cut=0.1;}{ n=split($7,a,"_"); ce=a[6];nce=1.0*a[7];cf=ce/(ce+nce); if (a[9]>=cut){p2++;} if (cf>=cut){p++;} nc++; if (cf>1) { # print ce,nce,cf,$0 >>"/dev/stderr" } cf2=sprintf("%d",1000*cf);occ[cf2]++; print $1,$2,$3,$4,$5,$6,$7,cf2,$9 }END{ # for (i in occ){print i" "occ[i]>>"/dev/stderr"} print p,nc,(100.0*p)/n>>"/dev/stderr" print p2,nc,(100.0*p2)/n>>"/dev/stderr" }