当我们对NGS数据call了snp之后一般会输出成vcf格式的数据，一行代表一个突变，例如

20      2451451 .       G       T       1939.77 .

AC=1;AF=0.500;AN=2;BaseQRankSum=-10.134;DP=239;Dels=0.00;FS=2.276;HaplotypeScore=0.0000;MLEAC=1;MLEAF=0.500;MQ=60.00;MQ0=0;MQRankSum=-0.258;QD=8.12;ReadPosRankSum=0.823;SOR=0.870

GT:AD:DP:GQ:PL  0/1:150,89:239:99:1968,0,3874

#前几列记录着该突变发生在第几号染色体以及该染色体的哪个坐标，我们的参考基因组在该位点是什么碱基，我们测到的突变成了什么碱基。

最后两列是测序深度以及正负测序深度，或者ref和allele的测序深度。

只有第8列是最复杂的，可以有高达几百个数据信息，取决于我们用什么样的软件来call的snp，以及call了snp之后用什么样的软件做的注释。

接下来我们还需要探究我们找到的突变是否在其它以及公开发表的数据库里面被找到过，所以可以下载非常多的公共数据库进行比对，我所见过的有一下一些，估计完全下载有0.5T

dbsnp144 （这个是ncbi提供的最权威的啦）

cgi69

ExAC.vcf.gz（这个是broadinstitute提供的外显子联盟）

Cosmic_v73.ann.vcf.gz （这个是癌症突变信息集）

finalTCGA.vcf.gz （TCGA计划也是癌症相关的）

icgc.vcf.gz

dbNSFP2.6vcf

SCLP.ann.vcf.gz

CCLE.ann.vcf.gz

ESP6500-SIv2.vcf.gz （Variants from the Exome Sequencing Project (ESP)）

adni-sum

safs-sum.indel.vcf.gz

gonl.vcf.gz

ssm.vcf.gz

ssi.vcf.gz

uk10k.vcf.gz

1000g-ph3v5.gff.gz  （千人基因组计划）

gwasCatalog.gff.gz  \

phewascatalog.gff.gz  \

dbgap-gwas.gff.gz  \

interproDomain.gff.gz \

clinvar.gff.gz \

RegulomeDB.gff.gz \

CancerGAMAdb.gff.gz \