24小时热门版块排行榜    

查看: 3500  |  回复: 18
本帖产生 1 个 MicEPI ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

秋电工

新虫 (小有名气)

[求助] 请各位来看棵树

近日学习用mega4.0建系统进化树,大概会操作了,但仍然有很多疑问。坛子里关于mega的帖子看了遍,迷茫发问的朋友多,解答的大牛人也有,我的问题和很多求助的虫友相似,现系统提出,渴望知道的大侠尽早看到此贴,尽早给若干小虫指点下迷津,不甚感激!
一,附件中的小树是用mega4.0,Bootsrap Value设置1000次建成。具体步骤如下:
a、打开mega软件—Alignment 下拉—Alignment exposure/clustal—create a new alignment --然后根据实际情况(DNA or 蛋白序列)选择Yes
b、Edit—Insert sequence from file(.fastal的格式);
c、Algin by clustal W---Data—Export Alignment—Mega fomat;
d、打开上一步保存的文件—Distance—compute pairwise—compute;
e、Phylogeny—Bootstrap Test of phylogeny—Neighbor-Joining—将Boostrap(500 replicates;seed=64238修改为1000 replicates)--compute---获得进化树—设好格式---image –save as TIF file
问题来了,在坛子里看见各位大侠给出的步骤,要用Clustal X软件进行序列比对之后才能建树,请问,mega4.0中是不是嵌合了clustal么,为什么有大侠说还要下载一个clustal软件呢?上面的建树步骤C是不是算经过了clustal呢?
二,附件中小树节点处数值的含义,请各位解释下为什么HX4和AB546196.1(bootstrap值100)与HX6和AJ308316.1(bootstrap值77)两类群聚一起后又97呢?是否意味着HX4和AB546196.1亲缘关系很近,甚至是同一个东西,HX6和AJ308316.1关系较近,这两类群又聚一起,而且bootstrap值97,大于77,很费解呀很费解。
三,树下的标尺有什么用?去算枝长吗?枝长又代表什么,进化距离吗?如果枝长是进化距离,那bootstrap值又代表什么?
四,在坛子里逛的时候,看到有大侠说“遗传标度法和步长值法是用来表示两种不同的结果时用的 步长值常用在鉴定菌种的进化树中 而遗传距离法常用在表示不同菌株间进化关系的进化树中 ”。什么是遗传标度法,什么是步长值法?本人建树目的是想得到菌株间的进化关系
五,看文献过程中,发现有部分文献的进化树要么只有标尺,要么只有bootstrap值,而且也没有对数值,标尺的含义做出解释,究竟什么样的树才算规范?
昨晚码了一个求助贴,满怀希望发的时候,系统维护,伤啊伤!今天抖擞一下,再次情绪饱满的发一个,敬候各位的佳音哈!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

邱宇

新虫 (初入文坛)

大侠有mega 4软件包吗?求分享!!!
18楼2016-05-16 12:26:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答

xp198766

铁杆木虫 (著名写手)

小木虫职业打酱油滴~~!

【答案】应助回帖

★ ★ ★
秋电工(金币+6): 很感谢,解答了我大部分问题,但是方法中有些具体的操作我还是不会,希望能更详尽 2011-04-18 11:08:26
scelab(金币+3, EPI+1): good 2011-04-18 11:54:10
1、clustalx 也是为了alignment,如果在MEGA里面alignment了,就不用clustalx了;

2、bootstrap值一般是将你的序列保留一部分,把剩下部分随机打乱,拼成不同的序列,组成1000个你的alignment文件,做树,显示的77表明,1000次做树的过程有有77%次,也就是770次得到HX6和AJ308316.1聚在一起这个结果;

3、下面的标尺是枝长,也是进化距离,bootstrap值是可信度;

4、这个我水平有限,也比较迷惑,个人觉得你现在做成的树适合的进化研究;

5、一般是要有bootstrap值的,这个表示可信度,要是值低于50,一般别人是不认同的,有的文章没有标尺,可能是只想得到树的分枝情况,没有想要计算各个物种的距离,只要得到的一个树的拓扑结果吧,但是bootstrap值一般是要的。

对了,你的alignment,排列好了之后,一般要将首尾序列截成一样齐……

只是个人意见,可能有偏差,有错误欢迎指出。
3楼2011-04-18 10:52:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

秋电工

新虫 (小有名气)

引用回帖:
Originally posted by xp198766 at 2011-04-18 10:52:28:
1、clustalx 也是为了alignment,如果在MEGA里面alignment了,就不用clustalx了;

2、bootstrap值一般是将你的序列保留一部分,把剩下部分随机打乱,拼成不同的序列,组成1000个你的alignment文件,做树,显示 ...

“对了,你的alignment,排列好了之后,一般要将首尾序列截成一样齐……
”弱弱问问大侠,这个具体怎么操作?mega4.0里嵌合的clustal没有截序列的功能么?
4楼2011-04-18 11:13:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xp198766

铁杆木虫 (著名写手)

小木虫职业打酱油滴~~!

引用回帖:
Originally posted by 秋电工 at 2011-04-18 11:13:05:
“对了,你的alignment,排列好了之后,一般要将首尾序列截成一样齐……
”弱弱问问大侠,这个具体怎么操作?mega4.0里嵌合的clustal没有截序列的功能么?

MEGA里面的不知道有没有,没有怎么用MEGA,好像clustalx里面也没有,不过bioedit里面有……可以先在BIOEDIT里面截,再导入到MEGA里面,再转换成MEGA格式,做树
5楼2011-04-18 11:26:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯,如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[基金申请] 申请科研项目纯靠实力能成功吗? 200+5 David1982 2026-10-13 8/400 2026-10-14 21:38 by lhw_GPA4.0
[硕博家园] 中科院海洋所招收海洋材料相关领域博士研究生及研究助理 +3 sujiner 2026-10-08 4/200 2026-10-14 20:40 by qicaiyun
[职场人生] 40+中年人失业后,还能找到好工作吗? +3 xqy2009 2026-10-13 3/150 2026-10-14 20:17 by 43358052
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
信息提示
请填处理意见