【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。