码迷,mamicode.com
首页 > 其他好文 > 详细

求GC含量

时间:2017-01-12 13:13:15      阅读:191      评论:0      收藏:0      [点我收藏+]

标签:ati   ons   class   print   code   ems   ros   ict   pen   

代码如下:

from collections import OrderedDict
ID_name = ‘‘
temp_dict = OrderedDict()
with open(rD:\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含量

标签:ati   ons   class   print   code   ems   ros   ict   pen   

原文地址:http://www.cnblogs.com/nklzj/p/6275729.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!