前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >建树忽略了一个分号折腾我好久...

建树忽略了一个分号折腾我好久...

作者头像
Listenlii-生物信息知识分享
发布2021-07-30 15:09:07
8860
发布2021-07-30 15:09:07
举报
文章被收录于专栏:Listenlii的生物信息笔记

UNITE数据库做完序列比对之后用FastTree建树。

参数-nt才能识别为核酸,否则识别为蛋白。

且即使加了-nt,简并碱基还是识别不出来。

注意这里的X在序列中为N。

我把上面这些warning去掉之后,导入iTOL说不是正确的树格式。

代码语言:javascript
复制
Couldn't initialize the tree. Make sure the encoding is plain text ASCII, and that the tree is in Newick, Nexus or phyloXML format

又试着在R里面用read.tree读,报错左右括号的个数不同。

代码语言:javascript
复制
numbers of left and right parentheses in Newick string not equal

这也太奇怪了...

我先检查了数中的各种符号,发现错误原因为树文件中物种注释之间的分号(;)。

点号(.) 竖线(|)下划线(_)是允许的。

去掉分号后在R中就不报错了。

但是树直接为NULL了。

又对着文件仔细看了很久,突然意识到nwk格式的树最后有一个分号,要把这个再手动填上。

分号表示树的结束,因此树中间如果有分号只会读取到这里,就会报错左右括号数量不相同了~

折腾这么久的原因是忘了UNITE数据库注释信息本身带有分号;且全部替换后忘了最后还要再加上一个分号。

本文参与 腾讯云自媒体同步曝光计划,分享自微信公众号。
原始发表:2021-07-12,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 Listenlii 微信公众号,前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
相关产品与服务
数据库
云数据库为企业提供了完善的关系型数据库、非关系型数据库、分析型数据库和数据库生态工具。您可以通过产品选择和组合搭建,轻松实现高可靠、高可用性、高性能等数据库需求。云数据库服务也可大幅减少您的运维工作量,更专注于业务发展,让企业一站式享受数据上云及分布式架构的技术红利!
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档