看你的数据,应该是一个样品取了三个生物学重复(重复1,重复2,重复3),每一个重复里面的(1和2)应该是双端测序得到的read1和read2两个文件的统计结果,因为后面readnumber和bases数量(两个文件如 YX-001_1 和YX-001_2 )是一样的。 1. YX-001 的bases数可以 YX-001_1 和YX-001_2 相加; read数量你直接用16,054,...
回答于 2018-09-12 10:51
输入 以下命令就可以查看到hmmer的版本,查看帮助信息的时候就会打印版本信息,-h参数; hmmsearch -h biolinux当中的meme的版本信息可以用 -version 参数导出版本信息: /biosoft/meme/meme-v4.12.0/bin/meme -version
回答于 2018-09-12 09:09
如果一个基因的多条转录本都有结构域建议选择第一条转录本为基因的代表序列,如果只有一个转录本就用那个转录本就好; 例如下面: 更多可观看:《基因家族实操视频课程》 参考文献:https://bmcgenomics.biomedcentral.com/articles/10.1186/1471-2164-15-3
回答于 2018-09-11 18:59
以前没有遇到这个现象,检查下html的结果,是不是弄错了? 输入的文件弄混了
回答于 2018-09-10 21:47
可以参考这篇文章解决办法:https://www.omicsclass.com/article/700
回答于 2018-09-10 16:35
一般来说公开的物种的基因组序列,都有注释的,都可以从基因组网站下载到该物种所有基因的cds和PEP序列; 如果你自己组装的基因组,需要对基因组进行注释,就可以得到基因组的所有基因序列(CDS和PEP);接着就可以继续往下分析了;
回答于 2018-09-10 16:31
29767.m000208 这个是你的基因ID。你的截图中的ID,是作者对鉴定出来WRKY基因家族的基因重新编的号; 你也可以自己重新编号一下自己鉴定出来的某类基因家族;
回答于 2018-09-10 11:56
可以通过基因组网站搜索,然后批量下载相关基因,参照JGI:https://www.omicsclass.com/article/50
回答于 2018-09-06 10:19