【直播】我的基因组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 条评论
登录 后参与评论

相关文章

来自专栏张善友的专栏

Silverlight + Model-View-ViewModel (MVVM)

     早在2005年,John Gossman写了一篇关于Model-View-ViewModel模式的博文,这种模式被他所在的微软的项目组用来创建Expr...

2938
来自专栏pangguoming

Spring Boot集成JasperReports生成PDF文档

由于工作需要,要实现后端根据模板动态填充数据生成PDF文档,通过技术选型,使用Ireport5.6来设计模板,结合JasperReports5.6工具库来调用渲...

1.2K7
来自专栏陈仁松博客

ASP.NET Core 'Microsoft.Win32.Registry' 错误修复

今天在发布Asp.net Core应用到Azure的时候出现错误InvalidOperationException: Cannot find compilati...

4818
来自专栏大内老A

The .NET of Tomorrow

Ed Charbeneau(http://developer.telerik.com/featured/the-net-of-tomorrow/) Exciti...

30810
来自专栏杨龙飞前端

scrollto 到指定位置

2494
来自专栏闻道于事

js登录滑动验证,不滑动无法登陆

js的判断这里是根据滑块的位置进行判断,应该是用一个flag判断 <%@ page language="java" contentType="text/html...

6648
来自专栏我和未来有约会

Kit 3D 更新

Kit3D is a 3D graphics engine written for Microsoft Silverlight. Kit3D was inita...

2506
来自专栏Ceph对象存储方案

Luminous版本PG 分布调优

Luminous版本开始新增的balancer模块在PG分布优化方面效果非常明显,操作也非常简便,强烈推荐各位在集群上线之前进行这一操作,能够极大的提升整个集群...

3095
来自专栏一个爱瞎折腾的程序猿

sqlserver使用存储过程跟踪SQL

USE [master] GO /****** Object: StoredProcedure [dbo].[sp_perfworkload_trace_s...

2000
来自专栏飞扬的花生

jsencrypt参数前端加密c#解密

      写程序时一般是通过form表单或者ajax方式将参数提交到服务器进行验证,如何防止提交的请求不被抓包后串改,虽然无法说绝对安全却给非法提交提高了难度...

3859

扫码关注云+社区