题目网址在这里:http://rosalind.info/problems/hamm/
我对这道题的思路,用readlines()将两条序列的序列信息分别保存在两个变量中,然后再一一比对,比较适合这道题。但如果序列长度较长,行数较多,可能并不适用。
这是我的代码,可以参考下。
#_*_ coding: utf-8 _*
file = open('E:\\bioinfo\study\data\\text6.txt', 'r')
seq = file.readlines()
seq1, seq2 = seq[0].strip(), seq[1].strip()
mutation = [i for i in range(len(seq1)) if seq1[i] != seq2[i]]
print len(mutation)