WES(2)snp-calling

来源:生信菜鸟团评论2,198

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

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

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

WES(2)snp-calling-图片1WES(2)snp-calling-图片2 WES(2)snp-calling-图片3WES(2)snp-calling-图片4 WES(2)snp-calling-图片5

WES(2)snp-calling-图片6 WES(2)snp-calling-图片7 WES(2)snp-calling-图片8 WES(2)snp-calling-图片9WES(2)snp-calling-图片10 WES(2)snp-calling-图片11WES(2)snp-calling-图片12 WES(2)snp-calling-图片13WES(2)snp-calling-图片14 WES(2)snp-calling-图片15WES(2)snp-calling-图片16 WES(2)snp-calling-图片17WES(2)snp-calling-图片18 WES(2)snp-calling-图片19 WES(2)snp-calling-图片20WES(2)snp-calling-图片21

发表评论

匿名网友