多序列比对-蛋白序列保守结构域比对图绘制 genedoc
蛋白质分析，经常要分析蛋白质的保守结构域，通常要展示多序列比对结果才能说明不同序列间的保守性；看到别人分析的蛋白保守结构域都有一个漂亮的多序列比对结果；今天我们就来介绍绘制这种氨基酸比对图：
1. 获得蛋白质序列
这里我直接在NCBI上搜索一下WRKY基因家族基因的蛋白质序列来进行后续的分析演示，获取方法如下图所示，搜索之后随机勾选几个基因，然后把序列下载下来，下载方法如下图所示：
2. 多序列比对
下载下来的序列首先要进行多序列比对，我们采用MEGA软件里面的clustalw来进行多序列比对。操作方法如下：打开MEGA软件，然后点击align新建一个比对，注意选择蛋白质，具体步骤如下所示：
然后，把上一步下载的fasta文件导入到MEGA，之后就可以点击W，也就是clustalW软件来做多序列比对，基因少的话，这步很快就能完成：
比对完成之后，就可以导出fasta格式的比对结果文件：
3. geneDoc软件分析基因保守结构域
geneDoc软件的下载地址：https://github.com/karlnicholas/GeneDoc，（ gd322700.zip genedoc 软件）这个软件的安装十分简单，双击一直下一步就可以完成安装，这里就不赘述了；安装好了软件打开；然后就可以导入比对好的序列，注意对于fasta的比对结果，这里采用file->import的方式导入比对文件，如下所示：
导入之后，就可以得到文章开头的图片了：黑色表示100%显示,
4. 编辑美化结果：
圆圈1处，可显示和不显示相似按钮，圆圈2处为调出设置对话框：
4.1 编辑字号大小,显示长度
选择configure 中的project选项卡,进行显示样式设置:
4.2 编辑不同相似氨基酸颜色:
选择configure 中的shade选项卡,进行彩色设置:
怎么样？这样看起来就漂亮多了吧！自己赶快来试试吧！
结果优化与导出: https://www.omicsclass.com/article/734
更多分析《基因家族分析实操课程》
gd322700.zip genedoc 软件
更多生物信息课程：
1. 文章越来越难发？是你没发现新思路，基因家族分析发2-4分文章简单快速，学习链接：基因家族分析实操课程、基因家族文献思路解读
2. 转录组数据理解不深入？图表看不懂？点击链接学习深入解读数据结果文件，学习链接：转录组（有参）结果解读；转录组（无参）结果解读
3. 转录组数据深入挖掘技能-WGCNA，提升你的文章档次，学习链接：WGCNA-加权基因共表达网络分析
4. 转录组数据怎么挖掘？学习链接：转录组标准分析后的数据挖掘、转录组文献解读
5. 微生物16S/ITS/18S分析原理及结果解读、OTU网络图绘制、cytoscape与网络图绘制课程
6. 生物信息入门到精通必修基础课：linux系统使用、perl入门到精通、perl语言高级、R语言画图
7. 医学相关数据挖掘课程，不用做实验也能发文章：TCGA-差异基因分析、GEO芯片数据挖掘、 GEO芯片数据不同平台标准化 、GSEA富集分析课程、TCGA临床数据生存分析、TCGA-转录因子分析、TCGA-ceRNA调控网络分析
8.其他，二代测序转录组数据自主分析、NCBI数据上传、二代测序数据解读、
延伸阅读
R语言快速入门与提高 | 学基因家族分析 | 基因家族分析实操线下培训班 | Excel也能生成矢量高清大图 | 投稿时-文章中图片大小分辨率如何设置 | 绘一棵超酷炫的系统发育树 | 教你绘制绘制circos图 |
- 发表于 2018-10-19 10:07
- 阅读 ( 65008 )
- 分类：软件工具
