{
    if ($NF=="(forward-only)"){id=$1;du=0;nr++;}else{
	if (du==1){print $0;}
    }
 if ($1=="banded") {
     sc=substr($4,1,length($4)-1);
     pctid=1.0*substr($5,1,length($5)-1);
     ov=$10;split($13,a,":");
     quepos=substr(a[1],2);
     split(quepos,q,"-");
     refpos=substr(a[2],1,length(a[2])-1);
     split(refpos,r,"-");
     if ((q[1]==1)&&(r[1]==1)){
	 mms=ov*5-sc;
	 if ((q[2]>27)&&(mms<20)){
#	     print id,pctid,sc,ov,quepos,refpos,$13
#	     print $0;du=1;
	     aloc[q[2]]++;
	     aloc2[mms,q[2]]++;	     
	 }
     }
 }
}END{
    print nr;
    for (i in aloc){
	print i" "aloc[i]
    }
    for (i in aloc2){
	print i" "aloc2[i]
    }
 }