WES(2)snp-calling

  • A+
所属分类:Genomics

准备文件:下载必备的软件和参考基因组数据

1、软件

ps:还有samtools,freebayes和varscan软件,我以前下载过,这次就没有再弄了,但是下面会用到

2、参考基因组

3、参考 突变数据

第一步,下载数据

第二步,bwa比对

第三步,sam转为bam,并sort好

第四步,标记PCR重复,并去除

第五步,产生需要重排的坐标记录

第六步,根据重排记录文件把比对结果重新比对

第七步,把最终的bam文件转为mpileup文件

第八步,用bcftools 来call snp

第九步,用freebayes来call snp

第十步,用gatk     来call snp

第十一步,用varscan来call snp

下面的图片是按照顺序来的,我就不整理了

1472327118-5591-image00121472327116-4886-image0022 1472327116-5528-image00321472327116-9111-image0042 1472327116-9656-image0051

1472327116-1420-image006 1472327118-7718-image007 1472327118-9587-image008 1472327119-6954-image0091472327120-8000-image010 1472327120-5693-image0111472327121-5843-image012 1472327121-8484-image0131472327122-5709-image014 1472327122-2356-image0151472327123-1249-image016 1472327123-2935-image0171472327124-6595-image018 1472327124-3468-image019 1472327125-1368-image0201472327125-7829-image021

avatar

发表评论

:?: :razz: :sad: :evil: :!: :smile: :oops: :grin: :eek: :shock: :???: :cool: :lol: :mad: :twisted: :roll: :wink: :idea: :arrow: :neutral: :cry: :mrgreen: