Bioinformatics/base count: Difference between revisions
→Procedural
Line 2,236:
T: 155
TOT= 500</pre>
=== procedural ( dictionary version) ===
{{works with|Python| 3.10.5 }}
<lang Python>
"""
Python 3.10.5 (main, Jun 6 2022, 18:49:26) [GCC 12.1.0] on linux
Created on Sun 2022/08/17 11:19:31
"""
def DispCount ( D ):
S = f'\n\nBases :\n\n' + f''.join ( [ f'{i} =\t{D[i]:4d}\n' for i in D.keys() ])
return S
S = 'CGTAAAAAATTACAACGTCCTTTGGCTATCTCTTAAACTCCTGCTAAATG' \
'CTCGTGCTTTCCAATTATGTAAGCGTTCCGAGACGGGGTGGTCGATTCTG' \
'AGGACAAAGGTCAAGATGGAGCGCATCGAACGCAATAAGGATCATTTGAT' \
'GGGACGTTTCGTCGACAAAGTCTTGTTTCGAGAGTAACGGCTACCGTCTT' \
'CGATTCTGCTTATAACACTATGTTCTTATGAAATGGATGTTCTGAGTTGG' \
'TCAGTCCCAATGTGCGGGGTTTCTTTTAGTACGTCGGGAGTGGTATTATA' \
'TTTAATTTTTCTATATAGCGATCTGTATTTAAGCAATTCATTTAGGTTAT' \
'CGCCGCGATGCTCGGTTCGGACCGCCAAGCATCTGGCTCCACTGCTAGTG' \
'TCCTAAATTTGAATGGCAAACACAAATAAGATTTAGCAATTCGTGTAGAC' \
'GACCGGGGACTTGCATGATGGGAGCAGCTTTGTTAAACTACGAACGTAAT'
D = { 'A' : 0, 'C': 0 , 'G' : 0, 'T' : 0 }
D = { k : S.count ( k ) for k in sorted ( D.keys() ) }
print ( 'Sequence:\n\n')
print ( ''.join( [ f'{S [ k: k + 50 ]}\n' for k in range ( 0, len ( S ), 50 ) ]))
print ( f'{ DispCount( D ) } \n------------')
print ( f'Σ = \t {sum ( D.values() )}\n============\n')
"""
</lang>JPD 2022/08/17
Output:
Sequence:
CGTAAAAAATTACAACGTCCTTTGGCTATCTCTTAAACTCCTGCTAAATG
CTCGTGCTTTCCAATTATGTAAGCGTTCCGAGACGGGGTGGTCGATTCTG
AGGACAAAGGTCAAGATGGAGCGCATCGAACGCAATAAGGATCATTTGAT
GGGACGTTTCGTCGACAAAGTCTTGTTTCGAGAGTAACGGCTACCGTCTT
CGATTCTGCTTATAACACTATGTTCTTATGAAATGGATGTTCTGAGTTGG
TCAGTCCCAATGTGCGGGGTTTCTTTTAGTACGTCGGGAGTGGTATTATA
TTTAATTTTTCTATATAGCGATCTGTATTTAAGCAATTCATTTAGGTTAT
CGCCGCGATGCTCGGTTCGGACCGCCAAGCATCTGGCTCCACTGCTAGTG
TCCTAAATTTGAATGGCAAACACAAATAAGATTTAGCAATTCGTGTAGAC
GACCGGGGACTTGCATGATGGGAGCAGCTTTGTTAAACTACGAACGTAAT
Bases :
A = 129
C = 97
G = 119
T = 155
------------
Σ = 500
============
"""
===Functional===
| |||