服务器构建系统发育树
首先要用在线的MAFFT进行Alignment
MAFFT网址
- 复制粘贴所有的fasta序列到文本框中
- UPPERCASE / lowercase: 选择 Same as input
- Direction of nucleotide sequences: Help 选择 Adjust direction according to the first sequence (accurate enough for most cases)
- Output order: 选择 Same as input
- 最后点击Submit
- 我们将邮件中的结果复制粘贴,另存为一个.fasta的文件.
利用服务器建树
### 平台MAFFT比对# mafft –h
# High speed mafft in > out
# High accuracy 200条序列、2000个位点以内
# unsure mafft --auto in > out#命令:mafft --auto(strategy) 输入文件名 > 输出文件名
#eg. mafft --maxiterate 1000 --globalpair ITS\ order.fas > out.tre
#sz 输出文件mafft --maxiterate 1000 --globalpair all_compare_OTU.fa >all_compare_OTU_Aligned.fa
fasta格式转phy 参考博客文章
可在转换phylip格式时保留名称的所有字符命令:fasta2phy.pl -i 输入文件名.fasta输出文件名为: 输入文件名.phy 使用方法:
python fasta2phylip.py
建Fast树
fasttree -nt all_compare_OTU_Aligned.fa > all_compare_OTU.nwk
建Raxml树,RAxML建树命令
上载 .phy 文件
命令行:raxmlHPC-PTHREADS-SSE3 -f a -s ==输入文件名== -n tre -m ==GTRGAMMA== -x 1234 -# ==1000== -T ==40==
-n 输出文件后缀名
-m 模型
-x 1234随机数
-# bootstrap
-T CPU使用数(最大40,一般使用20左右)
输出文件为: RAxML_bipartitions.tre
nohup raxmlHPC-PTHREADS-SSE3 -f a -s alpha_pep_align_gb.phy -n tre -m GTRGAMMA -x 1234 -# 1000 -T 28
进化模型选择命令
java -jar /lustre/project/og04/shichunwei/biosoft/jmodeltest-2.1.7/jModelTest.jar -d jmodel.fasta -s 11 -i -g 4 -f -AIC -AICc -DT -BIC -a -tr 8 > jmodel.out
Bayes系统发育树
上载 .nex 文件(文件内需带Bayes命令行)
mb -h (help帮助文件)
命令: mb -i ==文件名== 回车
或mpirun -n ==数字== mb -i ==文件名== (mpirun为多线程运 行,n需为nchains的倍数,最大30)
输出文件: .tre
下载到本地
Note: ==红色字体==为自定义部分,根据自己的文件名进行修改。in red: the part named by custom, filename, number…
in blue: command
建树命令
mb -i isolated_Aligned.fasta.nex
离线跑树
nohup 命令(eg. mb -i )文件名& 回车 less nohup.out (输出文件)查看文件进程 或者 jobs
服务器构建系统发育树相关推荐
- vbn中使用的3种流程控制结构是_细菌进化树构建:从模式种序列下载到构建系统发育树一键搞定...
细菌进化树 • 构 建 细菌进化树构建:从模式种序列下载到构建系统发育树一键搞定 对于细菌新种或者新属的发现,总是那么让人期待,但是当我们批量获得16S序列后,逐一对这些尚不知分类地位的序列进行比对并 ...
- Linux多线程贝叶斯建树教程,构建系统发育树:贝叶斯法建树
写在前面:我对建树也是一知半解,这里只是想记录一下自己跟别人学习的建树方法,可能不具有普适性.但毕竟写在公众平台,大家自主参考. 用phylosuite进行.nex文件的准备,用在线建树网站CIPRE ...
- MEGA | 多序列比对及系统发育树的构建
MEGA是一个用于多序列比对和可视化.以及构建系统发育树的免费程序.自1993年发布以来,MEGA共更新9个版本 (没有第八.九版),今年发布的MEGA 11为处理更大的数据集进行了优化. 之前我们介 ...
- 构建单拷贝同源蛋白系统发育树,一条命令提序列!
由于基因重组与基因突变的存在,利用16S rRNA构建的系统发育树通常不够可靠,而使用基于全基因组的单拷贝同源蛋白系统发育树则不会存在这方面的困扰,16S rRNA序列也就1500个碱基左右,而全基因 ...
- MEGA软件——系统发育树构建方法(图文讲解) 转载
转载:http://www.plob.org/2012/12/02/4927.html 一.序列文本的准备 构树之前先将目标基因序列都分别保存为txt文本文件中(或者把所有序列保存在同一个txt文本中 ...
- 基于短核苷酸序列频率构建细菌系统发育树
谷禾健康 目前,原核生物的系统发育分类依赖于16S rRNA基因序列,这种序列在细菌中广泛存在且高度保守,但即使16S rDNA测序的结果显示某些细菌间序列相似性≥99%,但它们仍然不能被归为同一物种 ...
- 宏病毒组(五)|病毒基因组系统发育树
我们通过宏病毒组测序,分析鉴定到病毒的种类及丰度后,通常会进一步去研究样品中鉴定到的病毒之间的相关性和互作情况.病毒的进化地位等,以便能更好地了解其生物学特性,提升对检测到的病毒认知. 目前宏病毒组个 ...
- 使用GENBANK数据进行分子系统发育树的构建
一.引言 GENBANK是目前最大而权威的分子序列数据库,调用其中数据可以进行分子系统发育树的构建. 1.序列数据获取(以皿蛛系统发育树为例) 在GenBank中,每一个物种或阶元都有一个taxid, ...
- 使用MrBayes构建贝叶斯系统发育树【实践】
建树是分子系统发育与进化研究中绕不过去的一道坎,本文就本人实际经验介绍贝叶斯树的实际操作: 基于不同原理的系统发育树构建方法的比较 Bayesian Inference (BI)建树 需要使用的软件包 ...
- 如何构建古生物化石的系统发育树
构建古生物化石的系统发育树需要以下步骤: 收集数据: 首先要收集大量的古生物学数据,包括化石的形态特征和生物学特征. 数据分析: 利用各种生物学分析方法,对收集的数据进行分析. 构建分支: 根据分析结 ...
最新文章
- 微信小程序顶部tab切换以及滑动
- 核爆rpg手机版_好嗨游戏:不玩吃亏 ,20款全世界最佳移动RPG角色扮演游戏(上)...
- Cartographer学习笔记
- WinInet, WinHttp, Winsock, ws2_32的基本解释
- 飞康FreeStor推出OpenStack Cinder集成解决方案
- GitHub 热榜:文字识别神器,超轻量级中文 OCR!
- dup和dup2以及cat函数函数
- Java通过Pattern类使用正则表达式
- telnet实现本地回显
- 当一个GameObject有两个Collider组件时,Physics Material不起作用
- 洛谷1966 火柴排队
- java项目中使用kettle的JNDI连接配置
- 恒生电子offer含金量_收获8个Offer,来给大家分享备战和面试经验
- Varnish的基本应用详解
- webUploader大文件断点续传学习心得
- ROS 机器人操作系统:版本说明
- Flume 的使用场景详解
- Android 实现欢迎界面
- 详解Java 12新特性:switch表达式
- android微信qq分享,android 一键分享 QQ 微信
热门文章
- 阿里云服务器租用价格表,阿里云服务器优惠
- 什么是客户旅程分析(customer journey mapping)
- 什么是物联网技术
- 【Excel 教程系列第 15 篇】Excel 中的简单排序(升序 / 降序)、多条件排序、按颜色排序、自定义排序、以及巧用“升序“制作工资条
- Ubuntu下快速安装jdk、tomcat、mysql和Redis
- 小麦亩产一千八(kela)
- 如何禁止软件联网,防止软件自动更新
- java调用js tel,jQuery国际电话区号选择插件intlTelInput.js
- IOS 应用安全测试内容
- pandas使用dataframe的索引和数据列同时对dataframe数据进行排序、使用ascending参数指定索引和数据列的排序方向(升序或者降序)