方法1,分别用hmmer搜索,再取交集. 方法2,把三个文件cat合并在一起,得到all.hmm文件,直接用hmm搜索;但是得到的结果你也是需要筛选的; 最后,没有说一下就得到同时具有三个结构域的方法,都需要前期搜索,后期筛选才能等到。 更多分析见:学习链接:基因家族分析实操课程、基因家族文献思路解读
回答于 2020-03-18 09:16
由于vbox更新导致有些地方不兼容,建议安装老版本的vbox: 可以卸载vbox重新安装老版本的vbox: https://www.virtualbox.org/wiki/Download_Old_Builds_6_0
回答于 2020-03-18 09:07
可以用iTOL批量修改进化树名字:https://www.omicsclass.com/article/604
回答于 2020-03-17 22:36
没看到你前后代码的运行结果,不好推测原因; 你看看前一步gset是否正常生成;
回答于 2020-03-17 09:35
可参考这个:https://github.com/rajewsky-lab/mirdeep2/issues/43 Hello all! I am trying to work with miRDeep2 on a transcriptomic data (RNA-seq), however, this error pops up every-time I run the miRDeep2.pl script. No file mirdeep_runs/run_28_08_2019_t_15_49_37/tmp/precursors.fa_stack found My pa...
回答于 2020-03-16 16:56