这个问题可以写个perl脚本或者python脚本完成,分别用bioperl或者biopython包处理fasta文件会很方便,有相应课程; 不用bioperl 给你个写脚本思路:首先把ID的对应关系读到一个hash中,然后判断只接处理ID行,也就是带有>的行,其他行输出就行; 生物信息入门到精通必修基础课:、perl语言高级、python语言入门到精...
回答于 2020-03-19 18:44
这个问题可以写个perl脚本或者python脚本完成,分别用bioperl或者biopython包处理fasta文件会很方便,有相应课程; 生物信息入门到精通必修基础课:linux系统使用、biolinux搭建生物信息分析环境、实验室linux生信分析平台搭建、linux命令处理生物大数据、perl入门到精通、perl语言高级、R语言画图、R语言快速入门与提高、p...
回答于 2020-03-19 18:42
这个问题可以写个perl脚本或者python脚本完成,分别用bioperl或者biopython包就可以的,有相应课程; 生物信息入门到精通必修基础课:linux系统使用、biolinux搭建生物信息分析环境、实验室linux生信分析平台搭建、linux命令处理生物大数据、perl入门到精通、perl语言高级、R语言画图、R语言快速入门与提高、python语言入门...
回答于 2020-03-19 18:40
方法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