跳至内容
售后
用户工具
登录
站点工具
搜索
工具
显示页面
修订记录
反向链接
最近更改
媒体管理器
网站地图
登录
>
最近更改
媒体管理器
网站地图
您的足迹:
realphy构建snp进化树
编辑本页后请点击“保存”。请参阅
syntax
了解维基语法。只有在您能
改进
该页面的前提下才编辑它。如果您想尝试一些东西,请先到
playground
热身。
媒体文件
====== REALPHY ====== 使用REALPHY构建SNP进化树有两种渠道,一种通过在线网站 [[https://realphy.unibas.ch/realphy/]] 该网站可以上传基因组或GBK文件直接进行分析,具体参考示例可以参考网址Suplementary materials不分 或者可以使用本地软件进行分析 使用本地软件需要在输出文件目录下配置一个config.txt文件 <code> BOWTIE2 /PUBLIC/software/public/Alignment/bowtie2-2.0.6/bowtie2 BOWTIE2BUILDER /PUBLIC/software/public/Alignment/bowtie2-2.0.6/bowtie2-build TREEPUZZLE /TJPROJ5/META_ASS/script/yuxi/Share/tree-puzzle-5.3/bin/puzzle RAXML /TJPROJ1/META_ASS/soft/standard-RAxML-8.2.12/raxmlHPC-HYBRID-SSE3 Rscript /PUBLIC/software/public/System/R-2.15.3/bin/Rscript MaxPars /TJPROJ5/META_ASS/script/yuxi/Share/phylip-3.697/exe/dnapars PhyML /TJPROJ5/META_ASS/script/yuxi/Share/PhyML-3.1/PhyML-3.1_linux64 </code> 软件均已成功按照直接复制即可,或者以后涉及到分析需求可以直接貂绒该分析软件 **使用方法** 构建进化树 软件安装在conda下,首先激活环境 <code> source /TJPROJ1/META_ASS/soft/anaconda3/bin/activate REALPHY </code> 最基本用法,将参考基因组 reference(最好是完成的基因组) 放入 input 文件夹,可以是 fasta 和 genbank 格式,如果放置多个 fasta/gbk 文件,软件会随机选择其中一个作为参考基因组,并将其他序列作为比对序列。如果指定某个序列作为参考基因组,用 -ref 参数指定。fastq 数据作为比对数据放入 input 文件夹,如果是PE测序的数据,双端的fastq文件名要用 *_R1.fastq/*_R2.fastq 区分。fastq 数据建议先对测序数据进行一定的QC操作,尽可能减少测序错误引入的 SNP 位点。 <code> realphy input output -ref ref_genome </code> 对核心基因组构建进化树。默认的方式是对参比基因组全部区域比对,包括非编码区。如果希望仅对CDS进行比对,就可以通过添加 -genes 参数。如果添加了 -genes 参数,那么参考基因组文件格式必须是含有 CDS 的 genbank 格式数据。 <code> realphy input output -ref ref_genome -genes </code> 使用多个参考基因组。这是 REALPHY 的一个“卖点”,作者认为当参考基因组与比对基因组差异较大时,如>5%,就会造成进化树构建的偏差。所以通过选择多个参考基因组,平衡不同型别的基因组,减少远源基因组带来的比对差异。通过设置 -refN 来实现多个参考基因组,并用 -merge 参数合并结果。由于计算量较大,不建议使用太多参考基因组。 <code> realphy input output -ref1 ref1_genome -ref2 ref2_genome -ref3 ref3_genome -merge </code> 建议使用 RAxML 来构建基于 ML 的进化树。因此生成的进化树图为 “PolySeqOut/RAxML_bestTree.raxml”,将其复制到带 GUI 的计算机上,用画树软件生成树图。 <code> #figtree需要独立安装 figtree RAxML_bestTree.raxml </code> 也可以不选择参考基因组,会随机选择一条输入序列进行比对 <code> REALPHY_v113 <sequence folder> <output folder> [options] </code> <code> #参考数据路径 /TJPROJ1/META_ASS/script_Advanced_analysis/REALPHY </code> 更多使用方法可以参考官方网址 [[https://realphy.unibas.ch/docs/Realphy_manual.pdf]] 或下载以下文档 {{ :yuxi:realphy_manual.pdf |}}
保存
预览
取消
编辑摘要
当您选择开始编辑本页,即寓示你同意将你贡献的内容按下列许可协议发布:
CC Attribution-Share Alike 4.0 International
realphy构建snp进化树.txt
· 最后更改: 2023/09/27 11:46 由
yuxi
页面工具
显示页面
修订记录
反向链接
回到顶部