Bioinformatics/base count: Difference between revisions
→procedural ( dictionary version)
Line 2,252:
def DispCount ( D ) :
S = f'\n\nBases :\n\n' + f''.join ( [ f'{i} =\t{D [ i ]:4d}\n' for i in
return S
S = 'CGTAAAAAATTACAACGTCCTTTGGCTATCTCTTAAACTCCTGCTAAATGCTCGTGCTTTCCAATTATGTAAGCGTTCCGAGACGGGGTGGTCGATTCTG' \
'AGGACAAAGGTCAAGATGGAGCGCATCGAACGCAATAAGGATCATTTGATGGGACGTTTCGTCGACAAAGTCTTGTTTCGAGAGTAACGGCTACCGTCTT' \
'CGATTCTGCTTATAACACTATGTTCTTATGAAATGGATGTTCTGAGTTGGTCAGTCCCAATGTGCGGGGTTTCTTTTAGTACGTCGGGAGTGGTATTATA' \
'TTTAATTTTTCTATATAGCGATCTGTATTTAAGCAATTCATTTAGGTTATCGCCGCGATGCTCGGTTCGGACCGCCAAGCATCTGGCTCCACTGCTAGTG' \
'TCCTAAATTTGAATGGCAAACACAAATAAGATTTAGCAATTCGTGTAGACGACCGGGGACTTGCATGATGGGAGCAGCTTTGTTAAACTACGAACGTAAT'
All = set( S )
BoI = set ( [ "A","C","G","T" ] )
other = All - BoI
D = {k : S.count ( k ) for k in All}
print ( 'Sequence:\n\n')
print ( f'
print ( '' if ( other == set () ) else f'Other\t{sum ( [ D [ k ] for k in sorted ( other ) ] ):4d}\n\n' )
print ( f'Σ = \t {sum ( [ D[k] for k in set ( [ "A","C","G","T" ] ) ] ) } \n============\n')▼
▲ print ( f'Σ = \t {sum ( [ D [ k ] for k in
pass
Line 2,295 ⟶ 2,296:
if ( __name__ == '__main__' ) :
main () if LIVE else test ()
|