33 lines
1.1 KiB
Text
33 lines
1.1 KiB
Text
F basecount(dna)
|
||
DefaultDict[Char, Int] d
|
||
L(c) dna
|
||
d[c]++
|
||
R sorted(d.items())
|
||
|
||
F seq_split(dna, n = 50)
|
||
R (0 .< dna.len).step(n).map(i -> @dna[i .+ @n])
|
||
|
||
F seq_pp(dna, n = 50)
|
||
L(part) seq_split(dna, n)
|
||
print(‘#5: #.’.format(L.index * n, part))
|
||
print("\n BASECOUNT:")
|
||
V tot = 0
|
||
L(base, count) basecount(dna)
|
||
print(‘ #3: #.’.format(base, count))
|
||
tot += count
|
||
V (base, count) = (‘TOT’, tot)
|
||
print(‘ #3= #.’.format(base, count))
|
||
|
||
print(‘SEQUENCE:’)
|
||
V sequence = "\
|
||
CGTAAAAAATTACAACGTCCTTTGGCTATCTCTTAAACTCCTGCTAAATG\
|
||
CTCGTGCTTTCCAATTATGTAAGCGTTCCGAGACGGGGTGGTCGATTCTG\
|
||
AGGACAAAGGTCAAGATGGAGCGCATCGAACGCAATAAGGATCATTTGAT\
|
||
GGGACGTTTCGTCGACAAAGTCTTGTTTCGAGAGTAACGGCTACCGTCTT\
|
||
CGATTCTGCTTATAACACTATGTTCTTATGAAATGGATGTTCTGAGTTGG\
|
||
TCAGTCCCAATGTGCGGGGTTTCTTTTAGTACGTCGGGAGTGGTATTATA\
|
||
TTTAATTTTTCTATATAGCGATCTGTATTTAAGCAATTCATTTAGGTTAT\
|
||
CGCCGCGATGCTCGGTTCGGACCGCCAAGCATCTGGCTCCACTGCTAGTG\
|
||
TCCTAAATTTGAATGGCAAACACAAATAAGATTTAGCAATTCGTGTAGAC\
|
||
GACCGGGGACTTGCATGATGGGAGCAGCTTTGTTAAACTACGAACGTAAT"
|
||
seq_pp(sequence)
|