有我在,别害怕!利用DNAMAN进行核酸比对的教程来啦~

有我在,别害怕!利用DNAMAN进行核酸比对的教程来啦~

小果  生信果  2023-04-15 19:00:11

DNAman是一款专业的分子生物学软件,可以用于DNA序列分析、比对、修剪、重组、克隆等操作。

首先导入比对序列,小果还是用上次做进化树用到的序列:通过NCBIhttps://www.ncbi.nlm.nih.gov/获取尼罗罗非鱼(Oreochromis niloticus)、斑点雀鳝(Lepisosteus oculatus)、青鳉(Oryzias latipes)、红鳍东方鲀(Takifugu rubripes)、斑马鱼(Danio rerio)、小鼠(Mus musculus)和人(Homo sapiens)的SOCS1基因fasta序列。

打开DNAman9软件,选择“Sequence”菜单下的“Align sequences”选项。

在“Align sequences”窗口中,点击“Mutiple Sequence Alignment”按钮

注意最下方导入的序列是DNA序列还是蛋白序列,勾选相应的即可,点击右侧File导入fasta序列,下一步

继续下一页,在这个界面一般默认即可

   * K-tuple的大小一般设置在2-6之间,具体取决于序列的长度和变异性。通常情况下,较短的K-tuple能够更好地捕捉序列中的局部相似性,但可能会在比对过程中产生误配。而较长的K-tuple则能够更好地区分序列的差异,但可能会错过一些相似性区域。

     *Gap penalty和gap extension penalty:这两个参数控制空位的罚分,即空位出现的惩罚程度。一般情况下,gap penalty的值应该比gap extension penalty高,因为在实际生物序列中,空位出现的概率较低。建议初始值分别设置为10和0。

   * delay divergent seqs%”是多序列比对中的一个参数,其作用是控制在比对过程中跳过高度分歧的序列,以减少这些序列对比对结果的影响。

完成以后等待结果,在这个窗口中,可以查看序列的比对结果,包括比对得分、相似性矩阵、序列的同源性等信息。

可以使用“File”菜单下的“Save as”选项将比对结果保存为文本文件,也可以选中结果后使用“Edit”菜单下的“Copy”选项将比对结果复制到剪贴板中,以便于在其他软件中进行分析和处理。

到这里核酸比对就结束啦,可以看到保守性比较较高的区域都用相同的色彩标记出来了,移动上方的红色框线可以查看不同区域。

如果你还在为生信分析苦恼,苦苦学习却不得要领,一定要试一下这个线上平台,什么circos图、小提琴图、曼哈顿图只要上传数据统统不在话下!

 

 

生信人R语言学习必备

立刻拥有一个Rstudio账号

开启升级模式吧

(56线程,256G内存,个人存储1T)