basList = "ATCG";
  split("", basct); # {basct[bas]} is the number of times {bas} written.
  for (i = 1; i <= length(basList); i++) { basct[substr(basList,i,1)] = 0; }
  
  labList = "DEFN";
  split("", labct); # {labct[lab]} is the number of bases with label {lab} written.
  for (i = 1; i <= length(labList); i++) { labct[substr(labList,i,1)] = 0; }
  
  
  for (i = 1; i <= length(basList); i++) 
    { basct[substr(basList,i,1)] = 0; }
  for (i = 1; i <= length(basList); i++) 
    { basct[substr(basList,i,1)] = 0; }