序列多重比对工具:MUSCLE
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相关推荐
- WebMGA:超快的基因组序列聚类注释在线工具
文章目录 超快的基因组序列聚类注释在线工具WebMGA 序列聚类cd-hit 基因预测ORF prediction 蛋白COG注释 作者简介 参考文献 猜你喜欢 写在后面 超快的基因组序列聚类注释在线 ...
- Navicat Premium 12 for Mac(多重连接数据库管理工具)汉化破解版
本站第一时间为各位小伙伴带来最好的数据库管理工具Navicat Premium Mac全新汉化激活版!Navicat Premium for mac是一个可多重连接的数据库管理工具,Navicat的功 ...
- clustalw序列比对_几个多序列比对软件:Muscle,ClustalW和T-coffee的简单比较
ClustalW是现在用的最广和最经典的多序列比对软件, Muscle是速度最快的多序列比对之一,T-coffee是越来越受到欢迎的软件.它们各自有什么优点和缺点呢? ClustalW不仅可以用来做多 ...
- deepin启动盘制作工具_YUMI——多重引导制作工具
一.总述 YUMI,全称Your Universal Multiboot Installer,即你的通用多重引导安装工具,官网介绍其为MulitibootISOs的继任者,两者的主要区别在于,Mult ...
- maya渲染序列文件存在检测工具
一.说明: 渲染文件都是具有连续性的,有时渲染会丢失一些帧,但帧数多时,检查起来很麻烦耗时,所以,本工具可以允许大家快速的检测特定目录中的渲染序列的存在性. 二.使用方法: 1.在 prefix 前缀 ...
- 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 ...
- muscle多序列比对简单应用
没有深入研究,只想对多序列做个比对,然后分组,具体命令如下: 1.多序列比对(数据量不大,曾经用了一个挺大的,直接吐核,放弃了,网上有说多线程操作,没有用,照例吐) muscle -in seqs.f ...
- 【生信MOOC】生物序列比对工具——多序列比对
[生信MOOC]生物序列比对工具2--多序列比对 文章的文字/图片/代码部分/全部来源网络或学术论文,文章会持续修缮更新,仅供大家学习使用. 目录 [生信MOOC]生物序列比对工具2--多序列比对 1 ...
- clustalw序列比对_序列比对(二)
Homology is the central concept for all of biology.--David Wake. Science, 1994 前言 正如前面引用的这句话,同源性是生物学 ...
最新文章
- s1考试 图书管理系统 结构体版
- asp.net中jQuery $post用法
- cc2530期末试卷_无线传感网期末考试 zigbee期末考试试题
- JAVA基础——编程练习(二)
- TQuery组件的Open方法与ExecSQL的区别
- 用CSS画一些多边形状
- 二叉树——基本知识+python实现
- 1056. 组合数的和(15)-PAT乙级真题
- 数字电平(一):TTL和CMOS
- 使用 iPerf 测试路由器局域网性能
- html实现图片裁剪,【前端】图片裁剪(二)Jcrop实现裁剪
- 基于GIS和Python的百度地图街景爬取
- 双系统蓝牙键盘的共享配对解决办法的简要步骤:win + arch~IRK、LTK、ERand、EDIV
- 网站服务器迁移域名怎么配置文件,记录网站服务器搬家如何迁移 Let’s Encrypt 证书 | 老左笔记...
- 深圳博瑞得一定不要去
- 手撕 44 道 JavaScript 送命题
- c/c++面试官经常提的问题
- UniApp接入阿里云金融级实人认证服务
- 发光二极管和光敏二极管
- PyCharm专业版安装教程
热门文章
- Nature指数2021亚太区排名:7所中国高校挺进前10!看看有没有你的母校?
- unknown error 1130,unknown error 1045
- org.apache.ibatis.binding.BindingException: Parameter 'username' not found. Available parameters
- PNAS | 基因调控之深度学习揭示免疫细胞分化的调节机制
- 开源化学信息学库 :ScaffoldGraph
- oracle开发常用,oracle开发常用关键字
- layui上传图片列表展示_layui怎么做图片在列表中显示 | 学步园
- aes key长度_原创 | 浅谈Shiro反序列化获取Key的几种方式
- 从复旦博士生到985高校副教授,贫穷使人深知摆脱自卑的艰难
- MPB:南土所冯有智组-基于微量热曲线的微生物群落代谢特征分析