
基于Python的DNA全局与局部序列比对详解
5星
- 浏览量: 0
- 大小:None
- 文件类型:None
简介:
本文章详细解析了利用Python进行DNA全局及局部序列对比的方法和技术,为生物信息学研究提供有力支持。
程序能实现以下功能:
a. 允许用户自定义输入gap值以及两条需要比对的序列。
b. 计算并输出得分矩阵。
c. 输出序列比对结果。
d. 使用matplotlib绘制得分矩阵路径。
### 实现步骤
1. **用户输入**
- 输入自定义的gap值
- 输入需进行比较的第一条碱基序列(A, T, C, G),换行表示输入完成
- 输入需进行比较的第二条碱基序列(A, T, C, G),换行表示输入完成
2. **代码实现**
1. 获取用户输入的gap值、s和t。
2. 调用构建得分矩阵函数,得到得分矩阵及方向矩阵。
3. 将所得的得分矩阵与方向矩阵作为参数传递给回溯函数开始进行路径回溯。路径存储使用的是全局变量,储存的方向信息用于后续处理。
全部评论 (0)
还没有任何评论哟~


