33 lines
1.1 KiB
Plaintext
33 lines
1.1 KiB
Plaintext
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
|