VCF与基因型编码 - SpiralKernel/bioinfo-wiki-migrated GitHub Wiki
VCF与基因型编码
更新:2026-10-14;状态:可用;验证范围:文档核对;操作需按所用环境验证。
VCF 以位点为记录,保存参考与替代等位基因,以及每个样本的基因型和其他字段。
先读哪些列
前八列是 CHROM、POS、ID、REF、ALT、QUAL、FILTER、INFO。包含样本时,FORMAT 定义后续样本单元格里的字段顺序。例如 FORMAT=GT:DP 的 0/1:28 表示对应基因型及深度。
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT sample1
1 100 v1 A G . PASS . GT:DP 0/1:28
此段以空格展示列,实际 VCF 使用制表符且需要元信息头。
常见误读
GT 的 0 是 REF、1 是第一个 ALT、2 是第二个 ALT,不是“携带两份 ALT”。对于二等位二倍体记录,1/1 才是两份第一 ALT;. 表示缺失。0|1 的相位顺序不自动给出父母来源。
REF 不保证常见,ALT 也不保证有害或是效应等位基因。QUAL 是位点级信息,GQ 是基因型级质量,DP 的含义需看其位置和头部定义。
交接检查
记录参考组装、倍性、多等位拆分及归一化策略;跨文件匹配优先使用组装、染色体、位置、REF、ALT 的组合,而不仅是可能重复的 ID。见SNP与基因型是什么和参考基因组与注释版本。
字段的权威定义在 VCF 4.5 规范。
记录、等位与事件
一个位置可能含多个ALT,拆分后ID可重复;一个结构事件也可能有两个BND记录。记录条数不直接等于位点数、替代等位数或事件数,见SNVIndelMNV与复杂变异、VCF中的SVTYPEEND与BND。
多等位拆分必须同步处理GT和等位维度字段,不能仅删ALT。样本子集和基因型掩蔽后,AC/AN/AF等统计要按实际来源重算或核对;外部人群频率与本队列频率分开记录。实际例子见人工VCF贯通变异分类与提取。