0
点赞
收藏
分享

微信扫一扫

【Python小试】判断一条序列GC含量高低

题目:

随便给定一条序列,如果GC含量超过65%,则认为高。

编程:

from __future__ import division #整数除法

def is_gc_rich(dna):
length = len(dna)
G_count = dna.upper().count('G')
C_count = dna.upper().count('C')
GC_content = (G_count + C_count) / length
if GC_content > 0.65:
print('GC content is high')
else:
print("GC content is normal")

测试

dna="agcTacGAcatgcacgtgcac"
is_gc_rich(dna)

解析

Python提供了__future__模块,把下一个新版本的特性导入到当前版本,于是我们就可以在当前版本中测试一些新版本的特性。

主要解决python2版本中和python3不同的一些问题,如print、整数除法、with用法、absolute_import等。上例中即引入了整数除法。


Ref:​​​​​​https://www.liaoxuefeng.com/wiki/897692888725344/923030465280480​​



作者:Bioinfarmer

 若要及时了解动态信息,请关注同名微信公众号:Bioinfarmer。

举报

相关推荐

0 条评论