DNA序列

题目描述

一个DNA序列由A/C/G/T四个字母的排列组合组成。G和C的比例（定义为GC-Ratio）是序列中G和C两个字母的总的出现次数除以总的字母数目（也就是序列长度）。在基因工程中，这个比例非常重要。因为高的GC-Ratio可能是基因的起始点。

给定一个很长的DNA序列，以及要求的最小子序列长度，研究人员经常会需要在其中找出GC-Ratio最高的子序列。

输入描述:

输入一个string型基因序列，和int型子串的长度

输出描述:

找出GC比例最高的子串,如果有多个输出第一个的子串

输入：

AACTGTGCACGACCTGA
5

输出：

GCACG

代码如下：

a = input()
n = int(input())
m = len(a)
res = 0
re = []
for i in range(m - n + 1):
    buff = a[i:i+n].count("G")+a[i:i+n].count("C")
    if res < buff:
        res = buff
        re = a[i:i+n]
print(re)