对初步组装进行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相关推荐

  1. nextpolish安装_「三代组装」使用Pilon对基因组进行polish

    软件安装 官方提供了编译好的jar包,方便使用 wget https://github.com/broadinstitute/pilon/releases/download/v1.23/pilon-1 ...

  2. 「三代组装」Pacbio组装后如何用自身数据进行polish(更新版)

    之前那我由于需要对PacBio的组装结果进行polish,于是写了「三代组装」Pacbio组装后如何用自身数据进行polish.最近发现自己又有了需求,于是重新回顾了我之前写的这篇文章,但是在实践的时 ...

  3. 「三代组装」Pacbio组装后如何用自身数据进行polish

    三代数据由于其高错误率(目前应该是10%左右), 即便在组装前有一步纠错环节,但是组装得到序列依旧存在着许多错误,因此需要进行polish环节.polish分为两个层次,三代原始序列polish和二代 ...

  4. 「三代组装」使用新版Falcon进行三代测序基因组组装

    这里的新版指的是PacBio公司在2018年9月发布pb-assembly, 而这篇文章是在2018年9月30日发的. 今年早些时候在参加三代培训时,听说PacBio会在今年对Falcon进行一些改变 ...

  5. 米兔机器人履带机甲如何编程_「开箱组装」米兔积木机器人履带机甲

    原标题:「开箱组装」米兔积木机器人履带机甲 米兔积木机器人是一款可以自由编程的拼插类积木机器人,早在第一代发布的时候就引起了很多积木迷的注意,毕竟可拼插.可编程,这种积木的可玩性实在太高了. 而它的第 ...

  6. 米兔积木机器人与履带机甲零件差别_「开箱组装」米兔积木机器人履带机甲

    米兔积木机器人是一款可以自由编程的拼插类积木机器人,早在第一代发布的时候就引起了很多积木迷的注意,毕竟可拼插.可编程,这种积木的可玩性实在太高了. 而它的第二代"履带机甲"发布之后 ...

  7. 开发2年半,「纪念碑谷」的团队做了一款开放世界“观鸟”游戏

    如果你随便抓一个路人问他/她在手机上玩过最有艺术感的游戏是什么,相信大概率会得到的答案是「纪念碑谷」. 这款由英国Ustwo Games在2014年公布的作品极具个性和创意,在当时一众良莠不齐的移动游 ...

  8. Pacbio 纯三代组装复活草基因组

    对于植物等真核生物基因组来说,重复序列, 多倍体,高杂合度等特征在利用二代数据进行组装的时候都会有很大的问题: 利用二代数据组装出来的基因组,大多达不到完成图的水准,通常只是覆盖到编码蛋白的基因区域, ...

  9. 全球首个「活体机器人」生娃!100%青蛙基因,杀不死,可繁衍4世

      视学算法报道   编辑:桃子.David [新智元导读]活体机器人「生娃」?全球首次,实属罕见!2020年1月,美国科学家利用从青蛙胚胎中提取的活细胞,创造出全球首个活体机器人Xenobot.现在 ...

  10. 爆改古董卡西欧计算器!能联网、能聊天,「作弊神器」只要150块

    晓查 郭一璞 发自 凹非寺  量子位 报道 | 公众号 QbitAI 一台朴实无华的计算器,卡西欧牌. 看起来和千万学子们所使用的计算器一样,没什么特别的,还带着一些陈旧的使用痕迹. 但如果你用磁铁碰 ...

最新文章

  1. 动态阈值--大津法理解
  2. linux rpm 查看安装目录下,Linux如何修改RPM的安装目录?
  3. 配置Tomcat监听80端口配置Tomcat虚拟主机Tomcat日志
  4. linux安装Python36
  5. React Native 轻松集成分享功能( iOS 篇)
  6. lisp语言画阿基米德线_中国油画艺术表达着艺术家思想情感,而且展现了油画语言独特之美...
  7. java rhino 运行 js_Mozilla Rhino :如何从Java调用JS函数
  8. spring+mybatis+springMVC+redis缓存+mysql+bootstrap+异步提交----联系人小demo
  9. python变量类型有几种_python十种变量类型
  10. myeclipse9安装SVN插件
  11. 为什么Linux与众不同?
  12. 初步了解CUDA(零)
  13. Glide加载长图;WebView加载富文本(图片自适应屏幕大小)
  14. 惠普打印机驱动下载安装后不能使用,驱动人生解决方案
  15. 【螺钉和螺母问题】【算法分析与设计】假设我们有n个直径各不相同的螺钉以及n个相应的螺母...
  16. iOS 第三方库、插件、知名博客总结
  17. runge phenomenon(龙格现象)和过拟合
  18. 在小百合注册了一个帐号
  19. 热血江湖数据库MYSQL修改_热血江湖数据库修改终极教程-20070823更新
  20. 热红外探测器的全球与中国市场2022-2028年:技术、参与者、趋势、市场规模及占有率研究报告

热门文章

  1. Java NIO笔记05
  2. 程序员的必备装备——为健康加油
  3. LaTex 在线编辑编译,生成PDF —— Overleaf
  4. White Sheet(面积法)
  5. 四.驱动框架入门之LED(中)
  6. 使用python发邮箱遇到的问题
  7. base64-图片传输
  8. 迪文屏串口更新协议了解
  9. mysql 分组查询 语句_MySQL分组查询
  10. 俄勒冈州立大学计算机科学专业,2019上海软科世界一流学科排名计算机科学与工程专业排名俄勒冈州立大学排名第301-400...