{ if ($NF=="(forward-only)"){ print $0; 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]>24)&&(mms<10)){ # print id,pctid,sc,ov,quepos,refpos,$13 print $0;du=1; aloc[q[2]]++; } } } }END{ print nr; for (i in aloc){ print i" "aloc[i] } }