개요 #

VCF 파일은 Variant Call Format 으로 GATK tool kit의 결과 파일이며, Variation 정보를 가지고 있다. 1000 Genome Project와 같은 대규모 프로젝트와 함께 발전되었다.

파일 형식 #

VCF 파일을 열어보면 해당 데이터에 대한 설명이 ##을 통해 주석으로 달려있다. CHROM, POS, ID, REF, ALT, QUAL, FILTER, INFO, FORMAT, Sample1 등과 같은 헤더 정보가 있으며, 그 아래에 실제의 variant 정보를 확인할 수 있다.

##fileformat=VCFv4.0
##fileDate=20110705
##reference=1000GenomesPilot-NCBI37
##phasing=partial
##INFO=<ID=NS,Number=1,Type=Integer,Description="Number of Samples With Data">
##INFO=<ID=DP,Number=1,Type=Integer,Description="Total Depth">
##INFO=<ID=AF,Number=.,Type=Float,Description="Allele Frequency">
##INFO=<ID=AA,Number=1,Type=String,Description="Ancestral Allele">
##INFO=<ID=DB,Number=0,Type=Flag,Description="dbSNP membership, build 129">
##INFO=<ID=H2,Number=0,Type=Flag,Description="HapMap2 membership">
##FILTER=<ID=q10,Description="Quality below 10">
##FILTER=<ID=s50,Description="Less than 50% of samples have data">
##FORMAT=<ID=GQ,Number=1,Type=Integer,Description="Genotype Quality">
##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">
##FORMAT=<ID=DP,Number=1,Type=Integer,Description="Read Depth">
##FORMAT=<ID=HQ,Number=2,Type=Integer,Description="Haplotype Quality">
#CHROM POS    ID        REF  ALT     QUAL FILTER INFO                              FORMAT      Sample1        Sample2        Sample3
2      4370   rs6057    G    A       29   .      NS=2;DP=13;AF=0.5;DB;H2           GT:GQ:DP:HQ 0|0:48:1:52,51 1|0:48:8:51,51 1/1:43:5:.,.
2      7330   .         T    A       3    q10    NS=5;DP=12;AF=0.017               GT:GQ:DP:HQ 0|0:46:3:58,50 0|1:3:5:65,3   0/0:41:3
2      110696 rs6055    A    G,T     67   PASS   NS=2;DP=10;AF=0.333,0.667;AA=T;DB GT:GQ:DP:HQ 1|2:21:6:23,27 2|1:2:0:18,2   2/2:35:4
2      130237 .         T    .       47   .      NS=2;DP=16;AA=T                   GT:GQ:DP:HQ 0|0:54:7:56,60 0|0:48:4:56,51 0/0:61:2
2      134567 microsat1 GTCT G,GTACT 50   PASS   NS=2;DP=9;AA=G                    GT:GQ:DP    0/1:35:4       0/2:17:2       1/1:40:3

Reference #

• http://en.wikipedia.org/wiki/Variant_Call_Format

Incoming Links #

Related Data Sciences (DataScience 0) #

• SeqHBase
• 파이썬/라이브러리/PyVCF

Related Articles (Article 1) #

• Jupyter/Practice
• KOBIC-Edu33

Related Bioinformaticses (Bioinformatics 2) #

• GATK
• Google Genomics
• SRA Excel
• 생물정보 데이터 형식

Suggested Pages #

• 0.262 PICARD
• 0.148 Calico
• 0.089 GBS read mapping/BWA
• 0.052 Hyungyong Kim
• 0.040 GATK/VariantCalling
• 0.027 GBS read mapping
• 0.025 Next-generation_sequencing
• 0.025 유전체
• 0.025 dbSNP
• 0.025 whole genome sequencing
• More suggestions...