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
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
|
##fileformat=VCFv4.2
##FILTER=<ID=PASS,Description="All filters passed">
##contig=<ID=1,length=150>
##ALT=<ID=*,Description="Represents allele(s) other than observed.">
##INFO=<ID=INDEL,Number=0,Type=Flag,Description="Indicates that the variant is an INDEL.">
##INFO=<ID=IDV,Number=1,Type=Integer,Description="Maximum number of raw reads supporting an indel">
##INFO=<ID=IMF,Number=1,Type=Float,Description="Maximum fraction of raw reads supporting an indel">
##INFO=<ID=DP,Number=1,Type=Integer,Description="Raw read depth">
##INFO=<ID=VDB,Number=1,Type=Float,Description="Variant Distance Bias for filtering splice-site artefacts in RNA-seq data (bigger is better)",Version="3">
##INFO=<ID=RPBZ,Number=1,Type=Float,Description="Mann-Whitney U-z test of Read Position Bias (closer to 0 is better)">
##INFO=<ID=MQBZ,Number=1,Type=Float,Description="Mann-Whitney U-z test of Mapping Quality Bias (closer to 0 is better)">
##INFO=<ID=BQBZ,Number=1,Type=Float,Description="Mann-Whitney U-z test of Base Quality Bias (closer to 0 is better)">
##INFO=<ID=MQSBZ,Number=1,Type=Float,Description="Mann-Whitney U-z test of Mapping Quality vs Strand Bias (closer to 0 is better)">
##INFO=<ID=SCBZ,Number=1,Type=Float,Description="Mann-Whitney U-z test of Soft-Clip Length Bias (closer to 0 is better)">
##INFO=<ID=SGB,Number=1,Type=Float,Description="Segregation based metric, http://samtools.github.io/bcftools/rd-SegBias.pdf">
##INFO=<ID=MQ0F,Number=1,Type=Float,Description="Fraction of MQ0 reads (smaller is better)">
##INFO=<ID=I16,Number=16,Type=Float,Description="Auxiliary tag used for calling, see description of bcf_callret1_t in bam2bcf.h">
##INFO=<ID=QS,Number=R,Type=Float,Description="Auxiliary tag used for calling">
##FORMAT=<ID=PL,Number=G,Type=Integer,Description="List of Phred-scaled genotype likelihoods">
##INFO=<ID=SCR,Number=1,Type=Integer,Description="Number of soft-clipped reads (at high-quality bases)">
##FORMAT=<ID=SCR,Number=1,Type=Integer,Description="Per-sample number of soft-clipped reads (at high-quality bases)">
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT sample
1 64 . C <*> 0 . DP=2;SCR=0;I16=1,1,0,0,65,2117,0,0,120,7200,0,0,0,0,0,0;QS=1,0;MQ0F=0 PL:SCR 0,6,65:0
1 65 . C <*> 0 . DP=3;SCR=0;I16=2,1,0,0,96,3080,0,0,180,10800,0,0,2,2,0,0;QS=1,0;MQ0F=0 PL:SCR 0,9,90:0
1 66 . T <*> 0 . DP=4;SCR=0;I16=3,1,0,0,120,3606,0,0,240,14400,0,0,5,9,0,0;QS=1,0;MQ0F=0 PL:SCR 0,12,106:0
1 67 . A <*> 0 . DP=4;SCR=0;I16=3,1,0,0,117,3431,0,0,240,14400,0,0,9,23,0,0;QS=1,0;MQ0F=0 PL:SCR 0,12,102:0
1 68 . G <*> 0 . DP=5;SCR=0;I16=4,1,0,0,164,5398,0,0,300,18000,0,0,13,45,0,0;QS=1,0;MQ0F=0 PL:SCR 0,15,133:0
1 69 . G <*> 0 . DP=8;SCR=3;I16=6,2,0,0,264,9048,0,0,396,21072,0,0,56,662,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,180:3
1 70 . T <*> 0 . DP=9;SCR=3;I16=6,3,0,0,275,8431,0,0,456,24672,0,0,64,782,0,0;QS=1,0;MQ0F=0 PL:SCR 0,27,193:3
1 71 . C <*> 0 . DP=9;SCR=3;I16=6,3,0,0,321,11473,0,0,456,24672,0,0,73,919,0,0;QS=1,0;MQ0F=0 PL:SCR 0,27,217:3
1 72 . C <*> 0 . DP=9;SCR=3;I16=6,3,0,0,318,11254,0,0,456,24672,0,0,82,1074,0,0;QS=1,0;MQ0F=0 PL:SCR 0,27,216:3
1 73 . C <*> 0 . DP=10;SCR=3;I16=7,3,0,0,348,12146,0,0,516,28272,0,0,91,1247,0,0;QS=1,0;MQ0F=0 PL:SCR 0,30,224:3
1 74 . C <*> 0 . DP=10;SCR=3;I16=7,3,0,0,350,12278,0,0,516,28272,0,0,100,1388,0,0;QS=1,0;MQ0F=0 PL:SCR 0,30,228:3
1 75 . C G,<*> 0 . DP=12;SCR=4;I16=6,2,3,1,274,9444,132,4392,480,28800,119,3601,55,501,54,1046;QS=0.697201,0.302799,0;VDB=0.0533436;SGB=-0.556411;RPBZ=1.63069;MQBZ=-3.26599;MQSBZ=0.111111;BQBZ=-0.515877;SCBZ=3.23349;MQ0F=0 PL:SCR 73,0,170,97,182,252:4
1 76 . C <*> 0 . DP=12;SCR=4;I16=9,3,0,0,410,14076,0,0,599,32401,0,0,120,1726,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,233:4
1 77 . G <*> 0 . DP=12;SCR=4;I16=9,3,0,0,385,12495,0,0,599,32401,0,0,131,1927,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,225:4
1 78 . T <*> 0 . DP=12;SCR=4;I16=9,3,0,0,365,11201,0,0,599,32401,0,0,142,2150,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,213:4
1 79 . G <*> 0 . DP=12;SCR=4;I16=9,3,0,0,431,15825,0,0,599,32401,0,0,153,2395,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,247:4
1 80 . C <*> 0 . DP=12;SCR=4;I16=9,3,0,0,416,14894,0,0,599,32401,0,0,164,2662,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,241:4
1 81 . C <*> 0 . DP=12;SCR=4;I16=9,3,0,0,420,14802,0,0,599,32401,0,0,175,2951,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,243:4
1 82 . A <*> 0 . DP=12;SCR=4;I16=9,3,0,0,393,12893,0,0,599,32401,0,0,186,3262,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,232:4
1 83 . C <*> 0 . DP=12;SCR=4;I16=9,3,0,0,438,15990,0,0,599,32401,0,0,197,3595,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,246:4
1 84 . C <*> 0 . DP=12;SCR=4;I16=9,3,0,0,438,16014,0,0,599,32401,0,0,208,3950,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,247:4
1 85 . A <*> 0 . DP=12;SCR=4;I16=9,3,0,0,402,13470,0,0,599,32401,0,0,219,4327,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,230:4
1 86 . T <*> 0 . DP=12;SCR=4;I16=9,3,0,0,396,13074,0,0,599,32401,0,0,229,4675,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,225:4
1 87 . G <*> 0 . DP=12;SCR=4;I16=9,3,0,0,457,17411,0,0,599,32401,0,0,238,4992,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,255:4
1 88 . A <*> 0 . DP=12;SCR=4;I16=9,3,0,0,395,13003,0,0,599,32401,0,0,247,5327,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,224:4
1 89 . T <*> 0 . DP=12;SCR=4;I16=9,3,0,0,398,13204,0,0,599,32401,0,0,256,5680,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,226:4
1 90 . G <*> 0 . DP=12;SCR=4;I16=9,3,0,0,463,17867,0,0,599,32401,0,0,263,5949,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,255:4
1 91 . C <*> 0 . DP=12;SCR=4;I16=9,3,0,0,451,16953,0,0,599,32401,0,0,269,6181,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,253:4
1 92 . C <*> 0 . DP=12;SCR=4;I16=9,3,0,0,419,14667,0,0,599,32401,0,0,274,6374,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,240:4
1 93 . G <*> 0 . DP=12;SCR=4;I16=9,3,0,0,398,13262,0,0,599,32401,0,0,279,6577,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,228:4
1 94 . T <*> 0 . DP=12;SCR=4;I16=9,3,0,0,377,11857,0,0,599,32401,0,0,283,6739,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,217:4
1 95 . C <*> 0 . DP=12;SCR=4;I16=9,3,0,0,419,14695,0,0,599,32401,0,0,287,6909,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,238:4
1 96 . G <*> 0 . DP=12;SCR=4;I16=9,3,0,0,402,13550,0,0,599,32401,0,0,290,7036,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,231:4
1 97 . C <*> 0 . DP=12;SCR=4;I16=9,3,0,0,454,17178,0,0,599,32401,0,0,293,7169,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,254:4
1 98 . T <*> 0 . DP=12;SCR=4;I16=9,3,0,0,420,14714,0,0,599,32401,0,0,295,7259,0,0;QS=1,0;MQ0F=0 PL:SCR 0,36,240:4
1 99 . C <*> 0 . DP=8;SCR=0;I16=6,2,0,0,299,11183,0,0,480,28800,0,0,199,4951,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,211:0
1 100 . A <*> 0 . DP=8;SCR=0;I16=6,2,0,0,267,8927,0,0,480,28800,0,0,200,5000,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,193:0
1 101 . C <*> 0 . DP=8;SCR=0;I16=6,2,0,0,295,10879,0,0,480,28800,0,0,200,5000,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,208:0
1 102 . C <*> 0 . DP=8;SCR=0;I16=6,2,0,0,293,10737,0,0,480,28800,0,0,200,5000,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,207:0
1 103 . T <*> 0 . DP=8;SCR=0;I16=6,2,0,0,278,9676,0,0,480,28800,0,0,200,5000,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,201:0
1 104 . C <*> 0 . DP=8;SCR=0;I16=6,2,0,0,296,10952,0,0,480,28800,0,0,200,5000,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,208:0
1 105 . C <*> 0 . DP=8;SCR=0;I16=6,2,0,0,295,10891,0,0,480,28800,0,0,200,5000,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,209:0
1 106 . A <*> 0 . DP=8;SCR=0;I16=6,2,0,0,276,9524,0,0,480,28800,0,0,200,5000,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,194:0
1 107 . G <*> 0 . DP=8;SCR=0;I16=6,2,0,0,302,11410,0,0,480,28800,0,0,200,5000,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,219:0
1 108 . G <*> 0 . DP=8;SCR=0;I16=6,2,0,0,295,10879,0,0,480,28800,0,0,200,5000,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,208:0
1 109 . A <*> 0 . DP=8;SCR=0;I16=6,2,0,0,265,8779,0,0,480,28800,0,0,200,5000,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,185:0
1 110 . T <*> 0 . DP=8;SCR=0;I16=6,2,0,0,268,8980,0,0,480,28800,0,0,200,5000,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,188:0
1 111 . T <*> 0 . DP=8;SCR=0;I16=6,2,0,0,266,8846,0,0,480,28800,0,0,200,5000,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,186:0
1 112 . T <*> 0 . DP=8;SCR=0;I16=6,2,0,0,263,8647,0,0,480,28800,0,0,200,5000,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,184:0
1 113 . C <*> 0 . DP=8;SCR=0;I16=6,2,0,0,300,11256,0,0,480,28800,0,0,200,5000,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,212:0
1 114 . A <*> 0 . DP=8;SCR=0;I16=6,2,0,0,272,9252,0,0,480,28800,0,0,198,4902,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,198:0
1 115 . T <*> 0 . DP=8;SCR=0;I16=6,2,0,0,265,8779,0,0,480,28800,0,0,195,4759,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,185:0
1 116 . C <*> 0 . DP=8;SCR=0;I16=6,2,0,0,299,11179,0,0,480,28800,0,0,191,4573,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,211:0
1 117 . A <*> 0 . DP=8;SCR=0;I16=6,2,0,0,271,9187,0,0,480,28800,0,0,187,4395,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,197:0
1 118 . A <*> 0 . DP=8;SCR=0;I16=6,2,0,0,271,9187,0,0,480,28800,0,0,182,4176,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,197:0
1 119 . A <*> 0 . DP=8;SCR=0;I16=6,2,0,0,267,8935,0,0,480,28800,0,0,177,3967,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,193:0
1 120 . C <*> 0 . DP=8;SCR=0;I16=6,2,0,0,295,10879,0,0,480,28800,0,0,171,3719,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,210:0
1 121 . C <*> 0 . DP=8;SCR=0;I16=6,2,0,0,293,10735,0,0,480,28800,0,0,165,3483,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,207:0
1 122 . T <*> 0 . DP=8;SCR=0;I16=6,2,0,0,282,9946,0,0,480,28800,0,0,159,3259,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,205:0
1 123 . G <*> 0 . DP=8;SCR=0;I16=6,2,0,0,308,11864,0,0,480,28800,0,0,152,2998,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,223:0
1 124 . A <*> 0 . DP=8;SCR=0;I16=6,2,0,0,269,9049,0,0,480,28800,0,0,145,2751,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,189:0
1 125 . G <*> 0 . DP=8;SCR=0;I16=6,2,0,0,303,11483,0,0,480,28800,0,0,137,2469,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,220:0
1 126 . A <*> 0 . DP=8;SCR=0;I16=6,2,0,0,266,8854,0,0,480,28800,0,0,129,2203,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,188:0
1 127 . G <*> 0 . DP=8;SCR=0;I16=6,2,0,0,289,10547,0,0,480,28800,0,0,121,1953,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,210:0
1 128 . G <*> 0 . DP=8;SCR=0;I16=6,2,0,0,295,10879,0,0,480,28800,0,0,113,1719,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,210:0
1 129 . T <*> 0 . DP=8;SCR=0;I16=6,2,0,0,252,7952,0,0,480,28800,0,0,105,1501,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,178:0
1 130 . T <*> 0 . DP=8;SCR=0;I16=6,2,0,0,264,8714,0,0,480,28800,0,0,97,1299,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,185:0
1 131 . G <*> 0 . DP=8;SCR=0;I16=6,2,0,0,299,11199,0,0,480,28800,0,0,89,1113,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,215:0
1 132 . G <*> 0 . DP=8;SCR=0;I16=6,2,0,0,292,10666,0,0,480,28800,0,0,81,943,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,212:0
1 133 . G <*> 0 . DP=8;SCR=0;I16=6,2,0,0,292,10662,0,0,480,28800,0,0,73,789,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,208:0
1 134 . A <*> 0 . DP=8;SCR=0;I16=6,2,0,0,266,8850,0,0,480,28800,0,0,65,651,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,187:0
1 135 . G <*> 0 . DP=8;SCR=0;I16=6,2,0,0,292,10698,0,0,480,28800,0,0,57,529,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,211:0
1 136 . A <*> 0 . DP=8;SCR=0;I16=6,2,0,0,249,7793,0,0,480,28800,0,0,49,423,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,178:0
1 137 . A <*> 0 . DP=8;SCR=0;I16=6,2,0,0,252,7976,0,0,480,28800,0,0,41,333,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,181:0
1 138 . T <*> 0 . DP=8;SCR=0;I16=6,2,0,0,245,7547,0,0,480,28800,0,0,33,259,0,0;QS=1,0;MQ0F=0 PL:SCR 0,24,174:0
1 139 . C <*> 0 . DP=6;SCR=0;I16=5,1,0,0,212,7516,0,0,360,21600,0,0,27,199,0,0;QS=1,0;MQ0F=0 PL:SCR 0,18,159:0
1 140 . A <*> 0 . DP=5;SCR=0;I16=4,1,0,0,162,5266,0,0,300,18000,0,0,22,150,0,0;QS=1,0;MQ0F=0 PL:SCR 0,15,132:0
1 141 . A <*> 0 . DP=4;SCR=0;I16=3,1,0,0,132,4370,0,0,240,14400,0,0,18,110,0,0;QS=1,0;MQ0F=0 PL:SCR 0,12,116:0
1 142 . G <*> 0 . DP=4;SCR=0;I16=3,1,0,0,142,5090,0,0,240,14400,0,0,14,78,0,0;QS=1,0;MQ0F=0 PL:SCR 0,12,125:0
1 143 . G <*> 0 . DP=3;SCR=0;I16=2,1,0,0,105,3683,0,0,180,10800,0,0,11,53,0,0;QS=1,0;MQ0F=0 PL:SCR 0,9,99:0
1 144 . C <*> 0 . DP=3;SCR=0;I16=2,1,0,0,102,3486,0,0,180,10800,0,0,8,34,0,0;QS=1,0;MQ0F=0 PL:SCR 0,9,96:0
1 145 . A <*> 0 . DP=2;SCR=0;I16=2,0,0,0,63,1989,0,0,120,7200,0,0,6,20,0,0;QS=1,0;MQ0F=0 PL:SCR 0,6,58:0
1 146 . G <*> 0 . DP=2;SCR=0;I16=2,0,0,0,70,2450,0,0,120,7200,0,0,4,10,0,0;QS=1,0;MQ0F=0 PL:SCR 0,6,64:0
1 147 . A <*> 0 . DP=2;SCR=0;I16=2,0,0,0,60,1800,0,0,120,7200,0,0,2,4,0,0;QS=1,0;MQ0F=0 PL:SCR 0,6,55:0
1 148 . G <*> 0 . DP=1;SCR=0;I16=1,0,0,0,35,1225,0,0,60,3600,0,0,1,1,0,0;QS=1,0;MQ0F=0 PL:SCR 0,3,35:0
1 149 . G <*> 0 . DP=1;SCR=0;I16=1,0,0,0,34,1156,0,0,60,3600,0,0,0,0,0,0;QS=1,0;MQ0F=0 PL:SCR 0,3,34:0
|
