蛋白序列中氨基酸频率

【Homo sapiens】 GI386828

insulin_Homo_sapiens = "MALWMRLLPLLALLALWG\

PDPAAAFVNQHLCGSHLVEALYLVCGERGFFYTPKTRREAEDLQ\

VGQVELGGGPGAGSLQPLALEGSLQKRGIVEQCCTSICSLYQLENYCN"

Telomeres = "MSRRKRKGAREGYDFMQIRDLSRSVGQRAKL\

IGVILDLGFPRKSKGTDLCCCMRIIDETNHQIAMAVNIF\

EENVERLPRVAAVGDVIVLCCVEVKSFKGEVNATFDKRF\

SSFGLYKGKDGDDLDPYHVSSYFHHIREDESLIVKLRKWL\

MNFQPHEDSCNFPMLREIKEETSVNLACKILHFCEAAKD\

EWIIFAWDGTNTPPNVICSKLEEEINSPLPLQLEPLPLSREV\

LCTLPVVGSILRMTFDADLVKNHLHLLNVDKWVKFMNM\

RLKVVDGLWLGVFTPQSKLQYTPNEDGLIVERQRLSEEWL\

FPKPSFITEEVNQDHAIPVTLMTVLTHSEVTAKFKCVVRVVA\

ATPCQAENLLSSTGEYRMRLTLEDSTARIHAFVTAKDGEVLF\

DGYPDIDELTRKLNILLGVNEVKDAPRNPPWVCVCLKSFCV\

SKTDVWSSRTFKIFDTKIVGDT"

aadir = ['A','C','D','E','F','G','H','I','K','L','M','N','P','Q','R','S','T','V','W','Y']

for bds in aadir:

    num = Telomeres.count(bds)

    print(bds,' '*5,num)