当前位置: 首页 » NGS(二代测序了么) » CallVariant » samtools 和 bcftools call snp

samtools 和 bcftools call snp

[文章目录] x

bcftools和samtools类似,用于处理vcf(variant call format)文件和bcf(binary call format)文件。前者为文本文件,后者为其二进制文件。官网:
https://github.com/samtools/

更多资料:

Variant calling using SAMtools

1.bam文件都要经过处理

-mpileup 是samtools中call snp的工具
-P 指platform, 现在短reads测序一般是ILLUMINA。
-f 后跟参考序列,序列文件必须提前建好index。
-E, Extended BAQ(base alignment quality) computation, 如果有的话,会提高检测出MNPs的灵敏度,当然会轻微的减下特异度。
-g Compute genotype likelihoods and output them in the binary call format(BCF).
-D Output per-sample read depth
> 是将结果保存到samtools_result.bcf文件中
还有一个参数-b 如果bam文件较多,可以放到一个文件中,每行放一个,将含有bam名字的文件的文件名放到-b 后面就可以了

2.得到bcf文件以后,第二步执行命令:

目前用的最新版1.3,跟之前的命令有些变化

报错:

参考资料:
http://www.ithao123.cn/content-5998798.html


下一篇 : 已经是最后文章

上一篇 :

2 条评论

  • flsnike 2016-08-30 回复

    博主您好,我在您以前的新浪博客看到一则关于RNA-Seq数据分析的博文http://blog.sina.com.cn/s/blog_670445240101jl08.html
    感觉对我帮助挺大,请问在那个博文中您是在哪儿学习的 课程,我也想去学学,谢谢!感激不尽!

    • sam 2017-03-22 回复 Author

      mooc

发表评论

电子邮件地址不会被公开。 必填项已用*标注

$(document).ready(function(){ $("#article-index").css('display','none');});