34 lines
1.1 KiB
Text
34 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)
|