BEGIN{
f="/data/images/proton2/run425/patalano/RBCL-ids.txt"
while(getline<f){
    faid[$1]=substr($0,24);
}
close(f);
FS="\t";getline;
}{
    s[$2]++; n++;
}END{print "Number of aligened reads: "n;
    for (i in s){print i" "s[i]" "100.0*s[i]/n" pct "faid[i];}
 }