【直播】我的基因组53:几个找变异的软件的效果比较

随便找一个SNP-calling的综述就可以找到一大堆软件的评价,我这里强烈推荐A survey of tools for variant analysis of next-generation genome sequencing data : http://bib.oxfordjournals.org/content/15/2/256.short

如果你在找变异这一个知识点还很模糊,甚至可以花点时间来翻译一下,同时欢迎大家交流自己的学习心得(http://www.biotrainee.com/thread-109-1-1.html 论坛回帖即可)

我很早以前处理外显子数据的时候,就比较过几个软件的找变异的效果,这里就继续沿用上次的思路!

我比较了gatk,freebayes,bcftools,varscan,都是引用率比较高的而且经受了时间的考验的好软件。

【直播】我的基因组(四):计算资源的准备

它们的下载安装方法是:

  1. ## Download and install bcftools
  2. ## http://www.htslib.org/download/
  3. ## http://www.htslib.org/doc/bcftools-1.0.html
  4. cd ~/biosoft
  5. mkdir bcftools && cd bcftools
  6. wget https://github.com/samtools/bcftools/releases/download/1.3.1/bcftools-1.3.1.tar.bz2
  7. tar xvfj bcftools-1.3.1.tar.bz2
  8. cd bcftools-1.3.1
  9. make
  10. cp bcftools /home/jianmingzeng/biosoft/myBin
  11. ~/biosoft/myBin/bin/bcftools --help
  12. ## Download and install freebayes
  13. ## https://github.com/ekg/freebayes
  14. ## http://clavius.bc.edu/~erik/CSHL-advanced-sequencing/freebayes-tutorial.html
  15. cd ~/biosoft
  16. mkdir freebayes && cd freebayes
  17. ## wget -O freebayes-master.zip https://codeload.github.com/ekg/freebayes/zip/master
  18. ## unzip freebayes-master.zip
  19. wget http://clavius.bc.edu/~erik/freebayes/freebayes-5d5b8ac0.tar.gz
  20. tar xzvf freebayes-5d5b8ac0.tar.gz
  21. cd freebayes
  22. make
  23. ~/biosoft/freebayes/freebayes/bin/freebayes
  24. cd ~/biosoft
  25. ## https://sourceforge.net/projects/varscan/files/
  26. ## http://varscan.sourceforge.net/index.html
  27. mkdir VarScan && cd VarScan
  28. wget https://sourceforge.net/projects/varscan/files/VarScan.v2.3.9.jar
  29. ## GATK有一点特殊,需要注册才能下载,微信公众号后台回复GATK也可以,我传给你哈

我这里就直接秀我以前处理的3个外显子测序结果的这4个软件的差异吧!

明显可以看出,大部分的snp位点都至少有两到三个软件支持。所以,只有测序深度达到一定级别,用什么软件来做snp-calling,其实影响并不大。

如果软件有区别的那些位点,通常我们还是需要在IGV里用肉眼仔细辨别的,毕竟这是关乎于个人健康的大事呀!

参考:http://www.bio-info-trainee.com/1150.html(阅读原文即可)

文:Jimmy

图文编辑:吃瓜群众

原文发布于微信公众号 - 生信技能树(biotrainee)

原文发表时间:2017-01-27

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏码洞

梦想不大,一笔一划!掌阅iReader Smart超级智能本使用初体验

无纸化办公提了很多年了,但对于纸张的使用依然没有减少,无论是从资源的浪费还是对环境污染的考虑,我们很多人都还在不自觉地挥霍着。所以现在有电子书阅读器,有可以无限...

953
来自专栏纯洁的微笑

gps定位 | 一张照片是怎么出卖你的?

875
来自专栏人工智能快报

科学家研发出能治疗内耳肿瘤的微型机器人

2015年11月,德国研究人员宣布开发出了一种像蠕虫般的外科手术机器人,并声称可以由此彻底改变内耳肿瘤的治疗。该机器人由弗劳恩霍夫协会的人员开发,合作方包括德国...

3005
来自专栏理论坞

60个设计师必备APP(上)

尽管安卓手机上免费的app极具诱惑力,iPhone仍然是设计圈首选的智能手机。对设计师而言,好用的不仅是最好的iPad app,苹果的Marvel也是相当便利的...

532
来自专栏CDA数据分析师

Excel到底有多厉害?

来自数据冰山,CDA以获作者授权转载 许多高级程序员瞧不上VBA。因为程序员是有鄙视链的:汇编 >C >C++ >Python >Java及C#及PHP(这三者...

2067
来自专栏成猿之路

我们为什么都应该坚持写博客,意义何在?

我很早就想写一篇为什么我坚持写博客和文字的文章给大家洗洗脑,不会总结和不懂总结的程序员不会进步太快。

822
来自专栏知晓程序

父母总信健康谣言?腾讯官方出了款小程序,用真科学粉碎它们 | 亲儿子 #31

生活中如果生病了,总是希望清清楚楚地知道自己得了什么病,以及怎样才能最有效治地好它。

731
来自专栏即时通讯技术

iOS后台唤醒实战:微信收款到账语音提醒技术总结

微信为了解决小商户老板们在频繁交易中不方便核对、确认到账的功能痛点,产品MM提出了新版本需要支持收款到账语音提醒功能。本文借此总结了iOS平台上的APP后台唤醒...

781
来自专栏企鹅号快讯

Real time是如何帮助Pxiar 完成他们作品的

这里是说的real-time不是unity引擎渲染的那种real-time,是在Pixar的每个制作环节中无时无刻都在实现的real-time,包含了电影中复杂...

2097
来自专栏前端桃园

英语很好是怎么自学前端的?

我本科学的是国际贸易,乱选的专业。毕业后做了半年外贸,实在不喜欢,然后去做英文编辑了。第二份工作也很无聊,就是写英文软文,发表在国外的行业期刊上,给公司做广告。...

812

扫码关注云+社区