我有一本包含多个基因的字典。tuple(zip()) 给出每个位点的核苷酸。例如(A,A,A),(T,T,G)等。我正在尝试计算每个位点的核苷酸数量。这样站点 1 显示了 3 个 A,站点 2 显示了 2 个 T 和 1 G。当我运行我的代码时,它只添加到 A 而没有其他内容。
List = tuple(zip(*myDict.values()))
A = 0
T = 0
G = 0
C = 0
site = 0
for value in List:
site +=1
if 'A':
A += 1
elif 'T':
T += 1
elif 'G':
G += 1
else:
C =+ 1
print 'Site:', site
print 'A:', A
print 'T:', T
print 'G:', G
print 'C:', C