{
    o[$1]++;
}END{
    for (i in o){
	x=o[i];
	m[x]++;
	n++;
    }
    print "#reads "n"  #mappings "NR;
    for (i in m){
	if (i>1){nmm+=m[i];pmm+=100.0*m[i]/n;}
	print i" "m[i]" "100.0*m[i]/n
    }
    print "#multimapped reads "nmm" : ("pmm" pct)";
 }