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

来源:互联网 发布:java 局部变量 编辑:程序博客网 时间:2024/06/18 12:50

题目网址在这里: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)


阅读全文
0 0
原创粉丝点击