首页 > 代码库 > 求GC含量
求GC含量
代码如下:
from collections import OrderedDict ID_name = ‘‘ temp_dict = OrderedDict() with open(r‘D:\testpy\Rosalind\rosalind_gc.txt‘,‘r‘) as f: for line in f: if line.startswith(‘>‘): ID_name = line.strip() temp_dict[ID_name] = ‘‘ #之前一直不会将ID和seq存到序列中,这里给temp_dicr[ID_name]一个空值就是为了将ID存进字典里。 else: temp_dict[ID_name] += line for ID,seq in temp_dict.items():
GC_num = seq.count(‘G‘) + seq.count(‘C‘) GC_ration = GC_num/len(seq) print (ID + " ====> " +"%.8f"%GC_ration)
求GC含量
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。