存在四个列表A,T,C,G,如果seq的第i个字符等于A,则A列表的第i个数加一,T,C,G亦然,python程序
时间: 2024-05-22 09:14:28 浏览: 34
可以这样写:
A = [0] * len(seq)
T = [0] * len(seq)
C = [0] * len(seq)
G = [0] * len(seq)
for i in range(len(seq)):
if seq[i] == "A":
A[i] += 1
elif seq[i] == "T":
T[i] += 1
elif seq[i] == "C":
C[i] += 1
elif seq[i] == "G":
G[i] += 1
这个程序会遍历整个序列,根据当前字符是哪个碱基,更新对应的列表。最后,每个列表中的第i个数,表示seq中第i个字符是对应的碱基的次数。
阅读全文