python寻找等位基因

B1K-04-12和Akashinriki是大麦的两个材料。请根据提供的CDS文件和GFF文件，找出两个基因组中的等位基因，并且以等位基因的物理位置为坐标，绘制七条染色体的等位基因散点图并拟合直线。

要求和注意事项：
1，简单判定等位基因的方法：两个cds文件互相做blast（第一次以B1K-04-12.cds.fasta为库文件， Akashinriki.cds.fasta 为问询序列；第二次反之）。如果B1K-04-12基因组的cds序列A与Akashinriki基因组的cds序列B，在两次blast结果里A和B都是彼此相似程度最高的序列，且A、B在同一条染色体上，那可以假定A、B是等位基因。
2，七张图输出到同一个PDF文件中。

这个问题需要用到生物信息学和数据可视化相关的知识。以下是一个简单的 Python 脚本，用于寻找 B1K-04-12 和 Akashinriki 两个基因组中的等位基因，并且绘制七条染色体的等位基因散点图并拟合直线。需要用到 Biopython 库进行序列比对和解析 GFF 文件，以及 Matplotlib 库进行数据可视化。

import os
import subprocess
from Bio import SeqIO
from Bio.Blast.Applications import NcbiblastnCommandline
import matplotlib.pyplot as plt# 定义文件路径和名称
b1k_cds_file = "B1K-04-12.cds.fasta"
aka_cds_file = "Akashinriki.cds.fasta"
b1k_gff_file = "B1K-04-12.gff"
aka_gff_file = "Akashinrik

python寻找等位基因相关推荐

python 寻找MC宝藏的具体位置
python 寻找MC宝藏具体位置这是我第一次使用 Markdown编辑器 ,之前看到好多博主写的文章中代码是写在黑方框里,我这一次也来尝试一下. 不多说,进入正题首先声明一下,这只是在有藏宝图的 ...
python寻找近义词：预训练模型 nltk+20newsbydate / gensim glove 转 word2vec
本文用python寻找英文近义词(中文:https://github.com/huyingxi/Synonyms) 使用的都是预训练模型方法一.nltk+20newsbydate (运行时下载太慢/ ...
python 知乎接口_用Python寻找知乎最美小姐姐
原标题:用Python寻找知乎最美小姐姐导读:最近知乎老是给我推送两个问答,一个是「长得好看是种什么体验?」,另一个是「女朋友长得好看是怎样的体验?」. 所以,本文将讲解如何爬取知乎这两个问题的回答 ...
Python寻找两个有序数组的中位数
Python寻找两个有序数组的中位数审题: 找出意味着这是一个查找算法题算法复杂度log级别,就是提示你是二分查找二分查找实现一般为递归 (1)递归包括递归体 (2)终止条件思路: 定理: 有 ...
python 寻找完美数_寻找完美的鼠标
python 寻找完美数 I have a small problem. I'm always looking for great computer mice. I've tried a number ...
python怎么查找文件-Python寻找路径和查找文件路径的示例
Sys.path 指定用于模块搜索路径的字符串列表也可以通过sys模块的append方法在Python环境中增加搜索路径. Sys.path.append("/usr/bin/" ...
如何使用python寻找K线中的黄昏之星
之前写过如何寻找启明星.启明星是股票在下跌过程中出现的反转信号,以一个或者两根长阴线,一个小十字星,以及一个中阳线或者大阳线并排为主要的特点.预示着下跌势头的衰减和上涨势头的确立.而黄昏之星正好相反, ...
完美假期第一步：用Python寻找最便宜的航班
大数文摘出品编译:高延.M.Y.Love.胡笳.蒋宝尚 "你喜欢旅游吗?" 这个简单的问题经常会得到一个积极的回复甚至还会额外收到一个或两个冒险的故事.通常来讲,旅行是一种体验新 ...
完美假期第一步：用Python寻找最便宜的航班！
一位资深的商业数据分析师将会带你建立一个网络爬虫项目,帮助我们找到最优惠的价格 "你喜欢旅游吗?" Python资源共享群:484031800 这个简单的问题经常会得到一个积极的回 ...

python寻找等位基因

python寻找等位基因相关推荐

最新文章

热门文章