【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://blog.csdn.net/zzc15806/article/details/81133045
https://www.liaoxuefeng.com/wiki/897692888725344/923030465280480
本文来自博客园,作者:生物信息与育种,转载请注明原文链接:https://www.cnblogs.com/miyuanbiotech/p/12741462.html。若要及时了解动态信息,请关注同名微信公众号:生物信息与育种。
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】博客园社区专享云产品让利特惠,阿里云新客6.5折上折
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步