从截图看不出问题,我这只能推测:你的gff文件有可能有问题,检查一下: 1.gff文件是否包含两个物种的所有基因位置信息; 2.看你的染色体是数字,如果两个物种的染色体都是数字,应加个前缀对染色体进行区分; H1,H2...... T1,T2...... 3.用下MCScanX吧,更方便些:《MCScanX安装和使用》
回答于 2018-07-20 13:33
可到UCSC这里去下载对应的基因组和注释信息: http://hgdownload.soe.ucsc.edu/downloads.html 要浏览不同版本的基因组注释信息,可以在UCSC的这个地方: 地址:http://genome-asia.ucsc.edu/cgi-bin/hgGateway?hgsid=682957885_HpOgAL9oAwgRsbQzV0aH1BTKPNO1&redirect=manual&source=genome.ucsc.edu 操作...
回答于 2018-07-19 19:11
可以适当删除一些比对较差的GAP区域,保留比对结果中比较保守的区域,再进行系统进化树构建。至于系统进化树的好坏,可以变换下参数,再做几次,看看有没有符合预期的;
回答于 2018-07-19 10:51
基因组组装的时候,有些片段无法挂载到染色体上,因此除了染色体外还有一些片段留在参考基因组当中,有的基因组将这些片段之间加N链接成super序列,也有的直接不处理以scoffold出现,具体可查参考基因组网站的说明。
回答于 2018-07-19 10:46
可以适当删除一些比对较差的GAP区域,保留比对结果中比较保守的区域,再进行系统进化树构建。至于系统进化树的好坏,可以变换下参数,再做几次,看看有没有符合预期的;
回答于 2018-07-18 17:42
这条代码是 blast 比对: query_file:替换自己查询fasta序列文件,使用的是blastp,输入文件为蛋白质序列; database:输入blast比对数据库文件路径; xyz.blast:输出文件 其他参数默认照抄就好了 更详细的可观看:《基因家族视频课程》有详细操作说明。
回答于 2018-07-18 09:52
有很多软件可以实现,如下: MEGA、DNAStar(MegAlign)、DNAman、T-coffee( http://tcoffee.vital-it.ch/apps/tcoffee/index.html )、 ESPript(http://espript.ibcp.fr/ESPript/ESPript/); ESPript结果: MEGA比对结果:
回答于 2018-07-18 09:43
1.hg19和hg38不同的基因组版本,本身参考序列就有差异,反应在位置上也会有差异,正常的。真想看看差异,可以把不同版本基因组的对应的UTR序列截取下来看看有没有差异。
回答于 2018-07-17 14:06
你检查你的,输入文件,第一行是不是有问题,确定是不是66列;再看看文章:如何提问。现有的信息只能帮你到这里。
回答于 2018-07-12 17:39
一类基因家族的基因,往往都含有相同的保守结构域,这些保守结构域收录在pfam数据库当中。 基因家族的鉴定可采用hmmer搜索的方法: 首先你要知道你要研究的基因家族的Pfam编号,然后可到pfam数据库当中根据编号下载对应的保守结构域的隐马尔科夫模型,就可以用hmmer根据你下载的隐马尔科夫模型搜索鉴定自己物种中的基因家...
回答于 2018-07-12 17:22