https://www.bilibili.com/video/BV1aK4y1f7oS?p=18

 系统发育树

多行 连接成 一行

Muscle

MUSCLE是RC Edgar开发的序列多重比对(Multiple Sequence Alignment,MSA)工具

下载和相关说明地址为http://www.drive5.com/muscle/manual/

1、比对并保存比对结果为Fasta格式文件

muscle -in seqs.fa -out seqs.afa

对于大数据集可以使用

muscle -in seqs.fa -out seqs.afa -maxiters 2

MUSCLE默认使用高准确度的比对方式,若需要更快但精度较低的方法可以使用:

用于氨基酸序列比对

muscle -in seqs.fa -out seqs.afa -maxiters 1 -diags -sv -distance1 kbit20_3

用于核算序列比对

muscle -in seqs.fa -out seqs.afa -maxiters 1 -diags

2、将比对结果转换为CLUSTALW格式的文件

muscle -in seqs.fa -clw

参数

-in 输入文件必须为fasta格式,如果序列中存在gaps,gaps将会被丢弃

-out 输出文件

3、根据多重比对结果构建UPGMA树

muscle -maketree -in seqs.afa -out seqs.phy

4、根据多重比对结果构建NJ(Neighbor-Joining)树

muscle -maketree -in seqs.afa -out seqs.phy -cluster neighborjoining

5、对已有比对结果加入新的序列

已有一个msa结果,想加入一条新的序列

muscle -profile -in1 existing_msa.afa -in2 new_seq.fa -out combined.afa

6、如果序列有多条,则先对需要加入的序列进行多重比对,然后对两个多重比对结果进行比对(同下)

muscle -in new_seqs.fa -out new_seqs.afa

muscle -profile -in1 existing_aln.afa -in2 new_seqs.afa -out combined.afas

7、两个比对结果进行比对

muscle -profile -in1 one.afa -in2 two.afa -out both.afa

8、提炼已有MSA

muscle -in msa.afa -out refined_msa.afa -refine

序列多重比对工具:MUSCLE相关推荐

  1. WebMGA:超快的基因组序列聚类注释在线工具

    文章目录 超快的基因组序列聚类注释在线工具WebMGA 序列聚类cd-hit 基因预测ORF prediction 蛋白COG注释 作者简介 参考文献 猜你喜欢 写在后面 超快的基因组序列聚类注释在线 ...

  2. Navicat Premium 12 for Mac(多重连接数据库管理工具)汉化破解版

    本站第一时间为各位小伙伴带来最好的数据库管理工具Navicat Premium Mac全新汉化激活版!Navicat Premium for mac是一个可多重连接的数据库管理工具,Navicat的功 ...

  3. clustalw序列比对_几个多序列比对软件:Muscle,ClustalW和T-coffee的简单比较

    ClustalW是现在用的最广和最经典的多序列比对软件, Muscle是速度最快的多序列比对之一,T-coffee是越来越受到欢迎的软件.它们各自有什么优点和缺点呢? ClustalW不仅可以用来做多 ...

  4. deepin启动盘制作工具_YUMI——多重引导制作工具

    一.总述 YUMI,全称Your Universal Multiboot Installer,即你的通用多重引导安装工具,官网介绍其为MulitibootISOs的继任者,两者的主要区别在于,Mult ...

  5. maya渲染序列文件存在检测工具

    一.说明: 渲染文件都是具有连续性的,有时渲染会丢失一些帧,但帧数多时,检查起来很麻烦耗时,所以,本工具可以允许大家快速的检测特定目录中的渲染序列的存在性. 二.使用方法: 1.在 prefix 前缀 ...

  6. LoRDEC: a tool for correcting errors in long sequencing reads 纠正长序列读取错误的工具

    LoRDEC: a tool for correcting errors in long sequencing reads Eric Rivals 21st August 2015 Abstract ...

  7. muscle多序列比对简单应用

    没有深入研究,只想对多序列做个比对,然后分组,具体命令如下: 1.多序列比对(数据量不大,曾经用了一个挺大的,直接吐核,放弃了,网上有说多线程操作,没有用,照例吐) muscle -in seqs.f ...

  8. 【生信MOOC】生物序列比对工具——多序列比对

    [生信MOOC]生物序列比对工具2--多序列比对 文章的文字/图片/代码部分/全部来源网络或学术论文,文章会持续修缮更新,仅供大家学习使用. 目录 [生信MOOC]生物序列比对工具2--多序列比对 1 ...

  9. clustalw序列比对_序列比对(二)

    Homology is the central concept for all of biology.--David Wake. Science, 1994 前言 正如前面引用的这句话,同源性是生物学 ...

最新文章

  1. s1考试 图书管理系统 结构体版
  2. asp.net中jQuery $post用法
  3. cc2530期末试卷_无线传感网期末考试 zigbee期末考试试题
  4. JAVA基础——编程练习(二)
  5. TQuery组件的Open方法与ExecSQL的区别
  6. 用CSS画一些多边形状
  7. 二叉树——基本知识+python实现
  8. 1056. 组合数的和(15)-PAT乙级真题
  9. 数字电平(一):TTL和CMOS
  10. 使用 iPerf 测试路由器局域网性能
  11. html实现图片裁剪,【前端】图片裁剪(二)Jcrop实现裁剪
  12. 基于GIS和Python的百度地图街景爬取
  13. 双系统蓝牙键盘的共享配对解决办法的简要步骤:win + arch~IRK、LTK、ERand、EDIV
  14. 网站服务器迁移域名怎么配置文件,记录网站服务器搬家如何迁移 Let’s Encrypt 证书 | 老左笔记...
  15. 深圳博瑞得一定不要去
  16. 手撕 44 道 JavaScript 送命题
  17. c/c++面试官经常提的问题
  18. UniApp接入阿里云金融级实人认证服务
  19. 发光二极管和光敏二极管
  20. PyCharm专业版安装教程

热门文章

  1. Nature指数2021亚太区排名:7所中国高校挺进前10!看看有没有你的母校?
  2. unknown error 1130,unknown error 1045
  3. org.apache.ibatis.binding.BindingException: Parameter 'username' not found. Available parameters
  4. PNAS | 基因调控之深度学习揭示免疫细胞分化的调节机制
  5. 开源化学信息学库 :ScaffoldGraph
  6. oracle开发常用,oracle开发常用关键字
  7. layui上传图片列表展示_layui怎么做图片在列表中显示 | 学步园
  8. aes key长度_原创 | 浅谈Shiro反序列化获取Key的几种方式
  9. 从复旦博士生到985高校副教授,贫穷使人深知摆脱自卑的艰难
  10. MPB:南土所冯有智组-基于微量热曲线的微生物群落代谢特征分析