1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
##fileformat=VCFv4.0
##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">
##FORMAT=<ID=GQ,Number=1,Type=Integer,Description="Genotype Quality">
##FORMAT=<ID=DP,Number=1,Type=Integer,Description="Read Depth">
##FILTER=<ID=q10,Description="Quality below 10">
##INFO=<ID=HM2,Number=0,Type=Flag,Description="HapMap2 membership">
##INFO=<ID=GN,Number=1,Type=String,Description="Gene Name">
##INFO=<ID=DP,Number=0,Type=Integer,Description="Depth,etc">
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT A
1 100 id1_100 GTTT G 1806 q10 DP=5;GN=gene1;HM2 GT:GQ:DP 0/1:409:35
1 110 id2_110 C T,G 1792 PASS DP=6 GT:GQ:DP 0/1:245:32
1 110 id1_110 CAAA C 1792 PASS DP=6 GT:GQ:DP 0/1:245:32
1 120 . GA G 628 q10 DP=21 GT:GQ:DP 1/1:21:21
1 130 id1_130 G T 1016 PASS DP=7;HM2 GT:GQ:DP 0/1:212:22
1 130 id2_130 GAA GG 1016 PASS DP=7;HM2 GT:GQ:DP 0/1:212:22
1 140 id1_140 GT G 727 PASS DP=8 GT:GQ:DP 0/1:150:30
1 150 id1_150 TAAAA TA,T 246 PASS DP=9 GT:GQ:DP 1/2:12:10
1 160 . TAAAA TA,T 246 PASS DP=10 GT:GQ:DP 1/2:12:10
2 100 . GTTT G 1806 q10 DP=35 GT:GQ:DP 0/1:409:35
2 110 id2_110_150 CAAA C 1792 PASS GN=gene2;HM2 GT:GQ:DP 0/1:245:32
2 120 . GA G 628 q10 DP=21 GT:GQ:DP 1/1:21:21
2 130 . GAA G 1016 PASS DP=22 GT:GQ:DP 0/1:212:22
2 140 . GT G 727 PASS DP=30 GT:GQ:DP 0/1:150:30
2 150 . TAAAA TA,T 246 PASS DP=10 GT:GQ:DP 1/2:12:10
2 160 id2_160 TAAAA TA,TC,T 246 PASS DP=11;GN=gene3 GT:GQ:DP 0/2:12:10
|
