小鼠WGS数据分析准备工作
一般来说,可以选择最新版小鼠参考基因组(mm10)了,如果你实在有其它需求,也可以自行选择其它版本。...2015 ZALENDE_EiJ.mgp.v5.snps.dbSNP142.vcf.gz
这些vcf文件的理解,需要对小鼠这个实验动物背景有一点了解,实际上这个时候我们需要的vcf文件应该是来自于dbSNP数据库的...dbsnp $snp \
-O ${sample}_raw.vcf \
1>${sample}_log.HC
done
其实这样的shell脚本是很烂的, 因为这个小鼠全基因组数据太大...13, 14, 15, 16, 17, 18, 19, 2, 3, 4, 5, 6, 7, 8, 9, MT, X, Y]
也就是说我们给的vcf文件里面的染色体是没有chr这个前缀,可是我们给的参考基因组里面却有这个前缀