omicsgene
omicsgene - 生物信息
实名认证 专业认证

性别: 北京 - 北京市 注册于 2018-04-20

擅长:重测序,遗传进化,转录组,GWAS

生物信息

向TA求助
13823金币数
77620 经验值
428个粉丝
主页被访问 75471 次

3896 个回答

0 赞同

基因家族中的get_fa_by_id.pl,运行时候出现问题

脚本已经更新请重新下载课程参考资料中对应脚本。

回答于 2018-10-28 06:29

0 赞同

提取基因上游1500bp时异常

有些基因所在的染色体或scaffold太短,导致基因的上游不足1500bp提取,所以会报错,可以查看一下所有基因的位置信息。

回答于 2018-10-28 06:27

0 赞同

利用脚本[get data by id]时出现报错

课程脚本已经更新,请重新下载脚本并运行。

回答于 2018-10-28 06:25

0 赞同

系统发育树:基因结构比较散,长度差异大,如何绘制满意的系统发...

应该是序列差异太大,可以删除一些大的gap,或者删除一些差异大的基因: https://www.omicsclass.com/article/221 如果基因全长序列差异太大,建议截取结构域所在序列构建进化树; 如果要编辑MEGA输出的多序列比对结果建议用genedoc:https://www.omicsclass.com/article/502

回答于 2018-10-27 07:15

0 赞同

序列编辑

可以用geneDOC软件编辑多序列比对后的序列: https://www.omicsclass.com/article/502 

回答于 2018-10-27 07:11

0 赞同

老师好,我在做进化树的时候报错

应该是序列差异太大,可以删除一些大的gap,或者删除一些差异大的基因: https://www.omicsclass.com/article/221 如果基因全长序列差异太大,建议截取结构域所在序列构建进化树;

回答于 2018-10-24 17:42

1 赞同

热图绘制

数据左上角,必须是NAME,不然显示数字;https://www.omicsclass.com/question/434

回答于 2018-10-23 19:06

1 赞同

热图中基因名称的显示

应该是heatmapper吧,刚刚用了下,原来是输入文件左上角,必须是NAME才行,不然是数字; http://www.heatmapper.ca/expression/

回答于 2018-10-23 18:59

1 赞同

用perl脚本提取顺势作用元件的序列是显示zsh: killed ,是因为内...

是小麦的基因组吗?小麦基因组巨大,看你的内存总共4G,一下子读入内存肯定不行; 我们的脚本是把参考基因组一下子都读进内存的,如果参考基因组文件太大,就会内存不足; 建议用get_fa_by_id.pl脚本把参考基因组分隔成不同的染色体,然后分染色体提取序列;最后  合并结果;

回答于 2018-10-23 13:00

0 赞同

在转录组中找糖苷酶基因

有参转录组:应该有基因组,如果基因组是JGI或者ensembl上的,可以直接搜索基因功能的关键字,参考:https://www.omicsclass.com/article/50 或者  https://www.omicsclass.com/article/58  再有就是基因组完成,一般会对基因的功能进行注释,会有一部分注释信息在gff文件里面你可以下载下来看看。

回答于 2018-10-22 10:49