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

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


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值