Python生物信息学数据管理——第5、6章课后习题

#初学小白，仅供参考
5.1一个简单的字典

dict0={'UAA':'Stop','UAG':'Stop','UGA':'Stop','AUG':'Start','GGG':'Glycin'}

5.2对起始和终止密码子进行计数

rna=''
for line in open('a.fasta'):
    if not line.startswith('>'):
        rna=rna+line.strip()
for frame in range(3):
    count_start=0
    count_end=0
    for i in range(frame,len(rna),3):
        codon=rna[i:i+3]
        if codon=='AUG':
            count_start +=1
        elif codon in ['UAA','UAG','UGA']:
            count_end +=1
        else:
            pass
    print('第{}种翻译方法有起始密码子{}个,终止密码子{}个'.format(frame+1,count_start,count_end))   

5.3搜索PubMed摘要中的关键词

fo=open('abstract.txt','r')
txt=fo.read()
i=0
if 'calmodulin' in txt: