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