released/pIRS.old/bwasam/samflag.stat

   1 # zcat ./bwamask/mask110621_I263_FCB066DABXX_L8_HUMjrmRACDKAAPEI-3.sam.gz|head -n500000|perl -lane 'BEGIN {my %a;} next if /^@/;++$a{$F[1]}; END {print $_,"\t",$a{$_} for sort {$a<=>$b} keys %a}'|tee samflag.stat
   2
   3 01000001        65      515
   4 01000101        69      645
   5 01001001        73      1259
   6 01001101        77      7343
   7 01010001        81      832
   8 01010011        83      117732
   9 01011001        89      1260
  10 01100001        97      925
  11 01100011        99      118305
  12 01110001        113     509
  13 01110101        117     650
  14 10000001        129     515
  15 10000101        133     1259
  16 10001001        137     645
  17 10001101        141     7343
  18 10010001        145     925
  19 10010011        147     118305
  20 10011001        153     650
  21 10100001        161     832
  22 10100011        163     117732
  23 10110001        177     509
  24 10110101        181     1260