BEGIN{FS="\t"}{
    m=$8;n=0;
    if (length(m)>0){
	n=split(m,a,",");
	for (i=1;i<=n;i++){
	    split(a[i],b,":");
	    mm[b[2]]++;
	}
    }
    o[n]++;
    r++;
    s+=n;
}END{
    print "nreads "r
    print "nmm "s;
    print "avgmm "1.0*s/r;
    for (n in o){ print n" "o[n]}
    for (m in mm){
	print m,mm[m],100.0*mm[m]/s
    }
 }