用python计算字典中的基因突变

Question

我有以下格式的数据：

>；abc12
ATCGACAG

>；def34
ACCGACG公司

等等

我把每一个基因都存储在一个字典里，以gt；开头的行作为值。所以字典是类似于{'abc12'：'atcgag'等的东西

现在我希望能够比较每个基因，这样它就可以计算出每个位点的A、T、C或G的数量。在

我唯一能想到的就是把字典分解成每个核苷酸站点的列表，并使用zip（）和一个计数器。这是最好的方法吗？如果是的话，我如何将字典分解为每个站点的列表？在

Answer 1

s1 = "ATCGACAG"
s2 = "ACCGACG"   
alignment = [s1[i] if s1[i] == s2[i] else "-" for i in range(len(min([s1,s2],key=len)))]
print "".join(alignment)
A-CGAC-