BEGIN{FS="\t";} { m=$8; n=split(m,a,","); if (n>0){ s=$2; for (i=1;i<=n;i++){ split(a[i],b,":") oc[b[2]" "s]++;nm++; } } }END{ for (i in oc){print i" "oc[i]" "oc[i]*1.0/nm} }