您的位置:首页 > 编程语言 > Python开发

python练习题2 从两个DNA序列中找到突变位点的数量

2017-11-13 20:37 513 查看
题目网址在这里: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)
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
相关文章推荐