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]]++;
	}
	if (n==1){#at most 1 mm
	    mm1[b[2]]++;	    
	    mm2[b[2]]++;	    
	    mm3[b[2]]++;	    
	}
	if (n==2){#at most 2 mm's
	    mm2[b[2]]++;
	    mm3[b[2]]++;	    
	    split(a[1],b,":");
	    mm2[b[2]]++;
	    mm3[b[2]]++;	    	    
	}
	if (n==3){#at most 3 mm's
	    mm3[b[2]]++;	    
	    split(a[1],b,":");
	    mm3[b[2]]++;	    	    
	    split(a[2],b,":");
	    mm3[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
    }
    print "max 1mm"
    for (m in mm1){
	print m,mm1[m],100.0*mm1[m]/s
    }
    print "max 2mm"
    for (m in mm2){
	print m,mm2[m],100.0*mm2[m]/s
    }
    print "max 3mm"
    for (m in mm3){
	print m,mm3[m],100.0*mm3[m]/s
    }
 }