「三代组装」使用Pilon对基因组进行polish
对初步组装进行polish
以FASTA和BAM文件作为输入,根据比对结果对输入的参考基因组进行提高,包括
- 单碱基差异
- 小的插入缺失(indels)
- 较大的插入缺失或者block替换
- 填充参考序列中的N
- 找到局部的错误组装
最后输出polish后的FASTA文件
利用第二代数据和第三代数据进行混装(Hybrid assembly),这种方法充分发挥了第二代数据质量高和第三代数据片段长的优势,组装出来的结果质量普遍比较高。混装策略的其中一种是先用第二代数据对第三代片段进行纠错,然后将第三代长片段组装成contigs,最后基于第二代的mate-paired数据生成scaffolds,如Pilon[41]软件支持该方法。2016年Jeong-Sun等[42]使用后一种方法完成了韩国人基因组的组装,组装完成之后的contigs和scaffolds N50分别达到17.9 Mb和44.8 Mb。
consensus得到的fa与consensus+pilon得到的fa, 由于N50反映的是序列的长度,不能反映polion所做的修改单碱基差异等提升效果。Busco理论上可以反映提升效果,局部错误纠正后理论的busco值应该更大了。
参考来源:
https://www.jianshu.com/p/cceeb7d1f413
卢鹏, 金静静, 李泽锋, 等. 基于第三代测序技术的基因组组装方法及其在烟草中的应用. 烟草科技, 2018, 51(2): 87-94.
「三代组装」使用Pilon对基因组进行polish相关推荐
- nextpolish安装_「三代组装」使用Pilon对基因组进行polish
软件安装 官方提供了编译好的jar包,方便使用 wget https://github.com/broadinstitute/pilon/releases/download/v1.23/pilon-1 ...
- 「三代组装」Pacbio组装后如何用自身数据进行polish(更新版)
之前那我由于需要对PacBio的组装结果进行polish,于是写了「三代组装」Pacbio组装后如何用自身数据进行polish.最近发现自己又有了需求,于是重新回顾了我之前写的这篇文章,但是在实践的时 ...
- 「三代组装」Pacbio组装后如何用自身数据进行polish
三代数据由于其高错误率(目前应该是10%左右), 即便在组装前有一步纠错环节,但是组装得到序列依旧存在着许多错误,因此需要进行polish环节.polish分为两个层次,三代原始序列polish和二代 ...
- 「三代组装」使用新版Falcon进行三代测序基因组组装
这里的新版指的是PacBio公司在2018年9月发布pb-assembly, 而这篇文章是在2018年9月30日发的. 今年早些时候在参加三代培训时,听说PacBio会在今年对Falcon进行一些改变 ...
- 米兔机器人履带机甲如何编程_「开箱组装」米兔积木机器人履带机甲
原标题:「开箱组装」米兔积木机器人履带机甲 米兔积木机器人是一款可以自由编程的拼插类积木机器人,早在第一代发布的时候就引起了很多积木迷的注意,毕竟可拼插.可编程,这种积木的可玩性实在太高了. 而它的第 ...
- 米兔积木机器人与履带机甲零件差别_「开箱组装」米兔积木机器人履带机甲
米兔积木机器人是一款可以自由编程的拼插类积木机器人,早在第一代发布的时候就引起了很多积木迷的注意,毕竟可拼插.可编程,这种积木的可玩性实在太高了. 而它的第二代"履带机甲"发布之后 ...
- 开发2年半,「纪念碑谷」的团队做了一款开放世界“观鸟”游戏
如果你随便抓一个路人问他/她在手机上玩过最有艺术感的游戏是什么,相信大概率会得到的答案是「纪念碑谷」. 这款由英国Ustwo Games在2014年公布的作品极具个性和创意,在当时一众良莠不齐的移动游 ...
- Pacbio 纯三代组装复活草基因组
对于植物等真核生物基因组来说,重复序列, 多倍体,高杂合度等特征在利用二代数据进行组装的时候都会有很大的问题: 利用二代数据组装出来的基因组,大多达不到完成图的水准,通常只是覆盖到编码蛋白的基因区域, ...
- 全球首个「活体机器人」生娃!100%青蛙基因,杀不死,可繁衍4世
视学算法报道 编辑:桃子.David [新智元导读]活体机器人「生娃」?全球首次,实属罕见!2020年1月,美国科学家利用从青蛙胚胎中提取的活细胞,创造出全球首个活体机器人Xenobot.现在 ...
- 爆改古董卡西欧计算器!能联网、能聊天,「作弊神器」只要150块
晓查 郭一璞 发自 凹非寺 量子位 报道 | 公众号 QbitAI 一台朴实无华的计算器,卡西欧牌. 看起来和千万学子们所使用的计算器一样,没什么特别的,还带着一些陈旧的使用痕迹. 但如果你用磁铁碰 ...
最新文章
- 动态阈值--大津法理解
- linux rpm 查看安装目录下,Linux如何修改RPM的安装目录?
- 配置Tomcat监听80端口配置Tomcat虚拟主机Tomcat日志
- linux安装Python36
- React Native 轻松集成分享功能( iOS 篇)
- lisp语言画阿基米德线_中国油画艺术表达着艺术家思想情感,而且展现了油画语言独特之美...
- java rhino 运行 js_Mozilla Rhino :如何从Java调用JS函数
- spring+mybatis+springMVC+redis缓存+mysql+bootstrap+异步提交----联系人小demo
- python变量类型有几种_python十种变量类型
- myeclipse9安装SVN插件
- 为什么Linux与众不同?
- 初步了解CUDA(零)
- Glide加载长图;WebView加载富文本(图片自适应屏幕大小)
- 惠普打印机驱动下载安装后不能使用,驱动人生解决方案
- 【螺钉和螺母问题】【算法分析与设计】假设我们有n个直径各不相同的螺钉以及n个相应的螺母...
- iOS 第三方库、插件、知名博客总结
- runge phenomenon(龙格现象)和过拟合
- 在小百合注册了一个帐号
- 热血江湖数据库MYSQL修改_热血江湖数据库修改终极教程-20070823更新
- 热红外探测器的全球与中国市场2022-2028年:技术、参与者、趋势、市场规模及占有率研究报告