dna$ = "CGTAAAAAATTACAACGTCCTTTGGCTATCTCTTAAACTCCTGCTAAATG" +
       "CTCGTGCTTTCCAATTATGTAAGCGTTCCGAGACGGGGTGGTCGATTCTG" +
       "AGGACAAAGGTCAAGATGGAGCGCATCGAACGCAATAAGGATCATTTGAT" +
       "GGGACGTTTCGTCGACAAAGTCTTGTTTCGAGAGTAACGGCTACCGTCTT" +
       "CGATTCTGCTTATAACACTATGTTCTTATGAAATGGATGTTCTGAGTTGG" +
       "TCAGTCCCAATGTGCGGGGTTTCTTTTAGTACGTCGGGAGTGGTATTATA" +
       "TTTAATTTTTCTATATAGCGATCTGTATTTAAGCAATTCATTTAGGTTAT" +
       "CGCCGCGATGCTCGGTTCGGACCGCCAAGCATCTGGCTCCACTGCTAGTG" +
       "TCCTAAATTTGAATGGCAAACACAAATAAGATTTAGCAATTCGTGTAGAC" +
       "GACCGGGGACTTGCATGATGGGAGCAGCTTTGTTAAACTACGAACGTAAT"

NewMap basecount.i()

If OpenConsole("")
  For i = 1 To Len(dna$)
    If (i % 50) = 1
      Print(~"\n" + RSet(Str(i - 1), 5) + " : ")
    EndIf
    t$ = Mid(dna$, i, 1)
    basecount(t$) + 1
    Print(t$)
  Next

  PrintN(~"\n\n" + Space(2) + "Base  count")
  PrintN(Space(2) + ~"----  -----")
  ForEach basecount()
    PrintN(RSet(MapKey(basecount()), 5) + " : " + RSet(Str(basecount()), 5))
    sigma + basecount()
  Next
  PrintN(~"\n" + "Total = " + RSet(Str(sigma), 5))
  Input()
EndIf