文件1,需要匹配的字符串对应列表(制表符分割的纯文本文件):
1 NC_010443.5
2 NC_010444.4
3 NC_010445.4
4 NC_010446.5
5 NC_010447.5
6 NC_010448.4
7 NC_010449.5
8 NC_010450.4
9 NC_010451.4
10 NC_010452.4
11 NC_010453.5
12 NC_010454.4
13 NC_010455.5
14 NC_010456.5
15 NC_010457.5
16 NC_010458.4
17 NC_010459.5
18 NC_010460.4
X NC_010461.5
Y NC_010462.3
文件2,vcf文件:
NC_010443.5 112 . C T 496 PASS DP=261;VDB=0.00422737;SGB=25.2245;RPB=0.537737;MQB=1;MQSB=1;BQB=0.570213;MQ0F=0;ICB=0.886034;HOB=0.0502959;AC=6;AN=52;DP4=189,12,25,1;MQ=60 GT:PL 0/0:0,12,129 0/0:0,30,224 0/0:0,18,141 0/0:0,3,37 0/0:0,24,221 0/0:0,27,163 0/0:0,24,188 0/1:114,0,19 0/1:114,0,180 0/0:0,24,162 1/1:129,15,0 0/1:86,0,136 0/0:0,21,195 0/0:0,18,188 0/0:0,15,123 0/1:142,0,89 0/0:0,15,149 0/0:0,24,162 0/0:0,15,129 0/0:0,18,165 0/0:0,60,209 0/0:0,36,195 0/0:0,36,190 0/0:0,33,187 0/0:0,69,208 0/0:0,30,187
NC_010443.5 370 . G A 999 PASS DP=480;VDB=0.765525;SGB=-9.28185;RPB=0.816783;MQB=1;MQSB=1;BQB=0.884522;MQ0F=0;ICB=0.801775;HOB=0.0650888;AC=32;AN=52;DP4=113,62,139,120;MQ=60 GT:PL 1/1:255,33,0 0/1:223,0,168 1/1:255,30,0 0/0:0,30,255 1/1:255,45,0 1/1:255,51,0 1/1:255,33,0 0/0:0,63,255 0/1:219,0,195 0/1:244,0,183 0/1:154,0,250 0/1:178,0,140 1/1:255,45,0 0/1:209,0,132 0/1:81,0,139 0/1:255,0,151 1/1:255,66,0 0/1:208,0,200 1/1:255,54,0 0/1:178,0,100 0/1:138,0,197 0/0:0,51,255 1/1:255,60,0 0/1:202,0,255 0/1:220,0,210 0/1:183,0,231
NC_010443.5 827 . T C 854 PASS DP=533;VDB=0.00165659;SGB=31.1247;RPB=0.465534;MQB=1;MQSB=1;BQB=0.993696;MQ0F=0;ICB=0.0478148;HOB=0.0184911;AC=5;AN=52;DP4=279,188,28,15;MQ=60 GT:PL 0/0:0,42,255 0/0:0,72,255 0/0:0,60,255 0/0:0,51,255 0/0:0,57,255 0/0:0,78,255 0/0:0,57,255 0/1:204,0,216 0/1:99,0,205 0/0:0,48,255 0/1:162,0,106 0/1:227,0,122 0/0:0,33,255 0/0:0,57,255 0/0:0,66,255 0/1:241,0,255 0/0:0,45,255 0/0:0,69,255 0/0:0,48,249 0/0:0,39,255 0/0:0,69,255 0/0:0,108,255 0/0:0,87,255 0/0:0,72,255 0/0:0,60,255 0/0:0,51,255
NC_010443.5 3014 . G T 999 PASS DP=454;VDB=0.958589;SGB=-9.3863;RPB=0.937155;MQB=1;MQSB=1;BQB=0.0688294;MQ0F=0;ICB=0.801775;HOB=0.0650888;AC=32;AN=52;DP4=95,70,133,124;MQ=60 GT:PL 1/1:255,51,0 0/1:120,0,228 1/1:255,51,0 0/0:0,24,203 1/1:255,66,0 1/1:255,45,0 1/1:255,33,0 0/0:0,39,255 0/1:204,0,217 0/1:255,0,230 0/1:255,0,64 0/1:110,0,128 1/1:255,39,0 0/1:220,0,242 0/1:149,0,245 0/1:207,0,145 1/1:255,42,0 0/1:152,0,152 1/1:255,51,0 0/1:166,0,130 0/1:141,0,191 0/0:0,69,255 1/1:255,48,0 0/1:171,0,173 0/1:221,0,162 0/1:192,0,177
NC_010443.5 3392 . G A 999 PASS DP=536;VDB=0.706175;SGB=-1.74937;RPB=0.609963;MQB=1;MQSB=1;BQB=0.292344;MQ0F=0;ICB=0.0478148;HOB=0.0184911;AC=47;AN=52;DP4=24,19,238,203;MQ=60 GT:PL 1/1:255,45,0 1/1:255,54,0 1/1:255,51,0 1/1:243,27,0 1/1:255,54,0 1/1:255,54,0 1/1:255,33,0 0/1:90,0,255 0/1:255,0,230 1/1:255,54,0 0/1:192,0,100 0/1:7,0,238 1/1:255,42,0 1/1:255,54,0 1/1:255,42,0 0/1:97,0,177 1/1:255,54,0 1/1:255,48,0 1/1:255,57,0 1/1:255,27,0 1/1:255,90,0 1/1:255,72,0 1/1:255,75,0 1/1:255,105,0 1/1:255,84,0 1/1:255,108,0
NC_010443.5 3430 . C T 999 PASS DP=512;VDB=0.314882;SGB=-8.61715;RPB=0.608373;MQB=1;MQSB=1;BQB=0.932661;MQ0F=0;ICB=0.0732873;HOB=0.0266272;AC=46;AN=52;DP4=28,26,212,207;MQ=60 GT:PL 1/1:255,42,0 1/1:255,48,0 1/1:255,42,0 1/1:255,30,0 1/1:255,57,0 1/1:255,63,0 1/1:255,30,0 0/1:125,0,255 0/1:250,0,239 1/1:255,54,0 0/1:157,0,131 0/1:45,0,155 1/1:255,42,0 1/1:255,48,0 1/1:255,42,0 0/1:127,0,110 1/1:255,66,0 0/1:176,0,255 1/1:255,51,0 1/1:237,30,0 1/1:255,75,0 1/1:255,72,0 1/1:255,45,0 1/1:255,111,0 1/1:255,90,0 1/1:255,102,0
NC_010443.5 4036 . A G 962 PASS DP=508;VDB=0.3649;SGB=62.5434;RPB=0.806806;MQB=1;MQSB=1;BQB=0.46193;MQ0F=0;ICB=0.0478148;HOB=0.0184911;AC=5;AN=52;DP4=226,197,27,24;MQ=60 GT:PL 0/0:0,57,255 0/0:0,72,255 0/0:0,39,255 0/0:0,51,255 0/0:0,48,255 0/0:0,33,255 0/0:0,51,255 0/1:198,0,120 0/1:255,0,177 0/0:0,57,255 0/1:220,0,242 0/1:114,0,212 0/0:0,63,255 0/0:0,39,255 0/0:0,54,255 0/1:255,0,104 0/0:0,54,255 0/0:0,48,255 0/0:0,30,255 0/0:0,36,255 0/0:0,60,255 0/0:0,105,255 0/0:0,90,255 0/0:0,72,255 0/0:0,57,255 0/0:0,39,255
NC_010443.5 4056 . C T 967 PASS DP=507;VDB=0.171323;SGB=64.0362;RPB=0.82443;MQB=1;MQSB=1;BQB=0.938503;MQ0F=0;ICB=0.0478148;HOB=0.0184911;AC=5;AN=52;DP4=226,204,32,20;MQ=60 GT:PL 0/0:0,60,255 0/0:0,60,255 0/0:0,24,255 0/0:0,51,255 0/0:0,51,255 0/0:0,39,255 0/0:0,54,255 0/1:171,0,98 0/1:244,0,196 0/0:0,60,255 0/1:231,0,214 0/1:146,0,182 0/0:0,63,255 0/0:0,45,255 0/0:0,45,255 0/1:255,0,149 0/0:0,45,255 0/0:0,42,255 0/0:0,39,255 0/0:0,33,255 0/0:0,60,255 0/0:0,108,255 0/0:0,99,255 0/0:0,75,255 0/0:0,72,255 0/0:0,51,255
NC_010443.5 4285 . C A 999 PASS DP=514;VDB=0.280171;SGB=59.5578;RPB=0.894116;MQB=1;MQSB=1;BQB=0.403017;MQ0F=0;ICB=0.0478148;HOB=0.0184911;AC=5;AN=52;DP4=207,231,17,32;MQ=60 GT:PL 0/0:0,45,255 0/0:0,57,255 0/0:0,36,255 0/0:0,36,255 0/0:0,45,255 0/0:0,48,255 0/0:0,39,255 0/1:217,0,59 0/1:255,0,255 0/0:0,48,255 0/1:235,0,101 0/1:176,0,205 0/0:0,75,255 0/0:0,81,255 0/0:0,66,255 0/1:227,0,161 0/0:0,51,255 0/0:0,57,255 0/0:0,45,255 0/0:0,33,255 0/0:0,72,255 0/0:0,84,255 0/0:0,90,255 0/0:0,72,255 0/0:0,75,255 0/0:0,57,255
NC_010443.5 4606 . T C 999 PASS DP=494;VDB=0.0629205;SGB=-17.2915;RPB=0.463604;MQB=1;MQSB=1;BQB=0.302156;MQ0F=0;ICB=0.0478148;HOB=0.0184911;AC=47;AN=52;DP4=19,19,226,183;MQ=60 GT:PL 1/1:255,60,0 1/1:255,51,0 1/1:255,75,0 1/1:255,48,0 1/1:255,39,0 1/1:244,30,0 1/1:255,78,0 0/1:161,0,229 0/1:247,0,150 1/1:255,42,0 0/1:177,0,135 0/1:152,0,146 1/1:255,36,0 1/1:255,48,0 1/1:255,51,0 0/1:144,0,205 1/1:255,57,0 1/1:255,54,0 1/1:215,24,0 1/1:255,51,0 1/1:255,69,0 1/1:255,63,0 1/1:255,63,0 1/1:255,57,0 1/1:255,78,0 1/1:255,51,0
To do: 将文件2中第一列替换成文件1中相应的名称
使用Perl代码实现字符串匹配替换:
#!/usr/bin/perl
脚本运行结果:
# 查看帮助文档:
suhui@suhui-ThinkCentre-E74:phython_Snakmake$ perl quiz1.pl -h
Discription: This script is to replace chromosome names.
Contact: Su Hui <hsu@kangpusen.com>
Example: perl quiz1.pl -chr chrname.txt -vcf input.txt -o out.txt
Useage:
Options:
-chr <file> a file containning the original and new chromosome name list
-vcf <file> a vcf-like text file modified from a normal vcf
-o <file> result file
-h help
查看结果文件:
suhui@suhui-ThinkCentre-E74:phython_Snakmake$ perl quiz1.pl -chr input1.txt -vcf input2.txt -o out.txt
suhui@suhui-ThinkCentre-E74:phython_Snakmake$ less -S out.txt
1 112 . C T 496 PASS DP=261;VDB=0.00422737;SGB=25.2245;RPB=0.537737;MQB=1;MQSB=1;BQB=0.570213;MQ0F=0;ICB=0.886034;HOB=0.0502959;AC=6;AN=52;DP4=189,12,25,1;MQ=60 GT:PL 0/0:0,12,129 0/0:0,30,224 0/0:0,18,141 0/0:0,3,37 0/0:0,24,221 0/0:0,27,163 0/0:0,24,188 0/1:114,0,19 0/1:114,0,180 0/0:0,24,162 1/1:129,15,0 0/1:86,0,136 0/0:0,21,195 0/0:0,18,188 0/0:0,15,123 0/1:142,0,89 0/0:0,15,149 0/0:0,24,162 0/0:0,15,129 0/0:0,18,165 0/0:0,60,209 0/0:0,36,195 0/0:0,36,190 0/0:0,33,187 0/0:0,69,208 0/0:0,30,187
1 370 . G A 999 PASS DP=480;VDB=0.765525;SGB=-9.28185;RPB=0.816783;MQB=1;MQSB=1;BQB=0.884522;MQ0F=0;ICB=0.801775;HOB=0.0650888;AC=32;AN=52;DP4=113,62,139,120;MQ=60 GT:PL 1/1:255,33,0 0/1:223,0,168 1/1:255,30,0 0/0:0,30,255 1/1:255,45,0 1/1:255,51,0 1/1:255,33,0 0/0:0,63,255 0/1:219,0,195 0/1:244,0,183 0/1:154,0,250 0/1:178,0,140 1/1:255,45,0 0/1:209,0,132 0/1:81,0,139 0/1:255,0,151 1/1:255,66,0 0/1:208,0,200 1/1:255,54,0 0/1:178,0,100 0/1:138,0,197 0/0:0,51,255 1/1:255,60,0 0/1:202,0,255 0/1:220,0,210 0/1:183,0,231
1 827 . T C 854 PASS DP=533;VDB=0.00165659;SGB=31.1247;RPB=0.465534;MQB=1;MQSB=1;BQB=0.993696;MQ0F=0;ICB=0.0478148;HOB=0.0184911;AC=5;AN=52;DP4=279,188,28,15;MQ=60 GT:PL 0/0:0,42,255 0/0:0,72,255 0/0:0,60,255 0/0:0,51,255 0/0:0,57,255 0/0:0,78,255 0/0:0,57,255 0/1:204,0,216 0/1:99,0,205 0/0:0,48,255 0/1:162,0,106 0/1:227,0,122 0/0:0,33,255 0/0:0,57,255 0/0:0,66,255 0/1:241,0,255 0/0:0,45,255 0/0:0,69,255 0/0:0,48,249 0/0:0,39,255 0/0:0,69,255 0/0:0,108,255 0/0:0,87,255 0/0:0,72,255 0/0:0,60,255 0/0:0,51,255
1 3014 . G T 999 PASS DP=454;VDB=0.958589;SGB=-9.3863;RPB=0.937155;MQB=1;MQSB=1;BQB=0.0688294;MQ0F=0;ICB=0.801775;HOB=0.0650888;AC=32;AN=52;DP4=95,70,133,124;MQ=60 GT:PL 1/1:255,51,0 0/1:120,0,228 1/1:255,51,0 0/0:0,24,203 1/1:255,66,0 1/1:255,45,0 1/1:255,33,0 0/0:0,39,255 0/1:204,0,217 0/1:255,0,230 0/1:255,0,64 0/1:110,0,128 1/1:255,39,0 0/1:220,0,242 0/1:149,0,245 0/1:207,0,145 1/1:255,42,0 0/1:152,0,152 1/1:255,51,0 0/1:166,0,130 0/1:141,0,191 0/0:0,69,255 1/1:255,48,0 0/1:171,0,173 0/1:221,0,162 0/1:192,0,177
1 3392 . G A 999 PASS DP=536;VDB=0.706175;SGB=-1.74937;RPB=0.609963;MQB=1;MQSB=1;BQB=0.292344;MQ0F=0;ICB=0.0478148;HOB=0.0184911;AC=47;AN=52;DP4=24,19,238,203;MQ=60 GT:PL 1/1:255,45,0 1/1:255,54,0 1/1:255,51,0 1/1:243,27,0 1/1:255,54,0 1/1:255,54,0 1/1:255,33,0 0/1:90,0,255 0/1:255,0,230 1/1:255,54,0 0/1:192,0,100 0/1:7,0,238 1/1:255,42,0 1/1:255,54,0 1/1:255,42,0 0/1:97,0,177 1/1:255,54,0 1/1:255,48,0 1/1:255,57,0 1/1:255,27,0 1/1:255,90,0 1/1:255,72,0 1/1:255,75,0 1/1:255,105,0 1/1:255,84,0 1/1:255,108,0
1 3430 . C T 999 PASS DP=512;VDB=0.314882;SGB=-8.61715;RPB=0.608373;MQB=1;MQSB=1;BQB=0.932661;MQ0F=0;ICB=0.0732873;HOB=0.0266272;AC=46;AN=52;DP4=28,26,212,207;MQ=60 GT:PL 1/1:255,42,0 1/1:255,48,0 1/1:255,42,0 1/1:255,30,0 1/1:255,57,0 1/1:255,63,0 1/1:255,30,0 0/1:125,0,255 0/1:250,0,239 1/1:255,54,0 0/1:157,0,131 0/1:45,0,155 1/1:255,42,0 1/1:255,48,0 1/1:255,42,0 0/1:127,0,110 1/1:255,66,0 0/1:176,0,255 1/1:255,51,0 1/1:237,30,0 1/1:255,75,0 1/1:255,72,0 1/1:255,45,0 1/1:255,111,0 1/1:255,90,0 1/1:255,102,0
1 4036 . A G 962 PASS DP=508;VDB=0.3649;SGB=62.5434;RPB=0.806806;MQB=1;MQSB=1;BQB=0.46193;MQ0F=0;ICB=0.0478148;HOB=0.0184911;AC=5;AN=52;DP4=226,197,27,24;MQ=60 GT:PL 0/0:0,57,255 0/0:0,72,255 0/0:0,39,255 0/0:0,51,255 0/0:0,48,255 0/0:0,33,255 0/0:0,51,255 0/1:198,0,120 0/1:255,0,177 0/0:0,57,255 0/1:220,0,242 0/1:114,0,212 0/0:0,63,255 0/0:0,39,255 0/0:0,54,255 0/1:255,0,104 0/0:0,54,255 0/0:0,48,255 0/0:0,30,255 0/0:0,36,255 0/0:0,60,255 0/0:0,105,255 0/0:0,90,255 0/0:0,72,255 0/0:0,57,255 0/0:0,39,255
1 4056 . C T 967 PASS DP=507;VDB=0.171323;SGB=64.0362;RPB=0.82443;MQB=1;MQSB=1;BQB=0.938503;MQ0F=0;ICB=0.0478148;HOB=0.0184911;AC=5;AN=52;DP4=226,204,32,20;MQ=60 GT:PL 0/0:0,60,255 0/0:0,60,255 0/0:0,24,255 0/0:0,51,255 0/0:0,51,255 0/0:0,39,255 0/0:0,54,255 0/1:171,0,98 0/1:244,0,196 0/0:0,60,255 0/1:231,0,214 0/1:146,0,182 0/0:0,63,255 0/0:0,45,255 0/0:0,45,255 0/1:255,0,149 0/0:0,45,255 0/0:0,42,255 0/0:0,39,255 0/0:0,33,255 0/0:0,60,255 0/0:0,108,255 0/0:0,99,255 0/0:0,75,255 0/0:0,72,255 0/0:0,51,255
1 4285 . C A 999 PASS DP=514;VDB=0.280171;SGB=59.5578;RPB=0.894116;MQB=1;MQSB=1;BQB=0.403017;MQ0F=0;ICB=0.0478148;HOB=0.0184911;AC=5;AN=52;DP4=207,231,17,32;MQ=60 GT:PL 0/0:0,45,255 0/0:0,57,255 0/0:0,36,255 0/0:0,36,255 0/0:0,45,255 0/0:0,48,255 0/0:0,39,255 0/1:217,0,59 0/1:255,0,255 0/0:0,48,255 0/1:235,0,101 0/1:176,0,205 0/0:0,75,255 0/0:0,81,255 0/0:0,66,255 0/1:227,0,161 0/0:0,51,255 0/0:0,57,255 0/0:0,45,255 0/0:0,33,255 0/0:0,72,255 0/0:0,84,255 0/0:0,90,255 0/0:0,72,255 0/0:0,75,255 0/0:0,57,255
1 4606 . T C 999 PASS DP=494;VDB=0.0629205;SGB=-17.2915;RPB=0.463604;MQB=1;MQSB=1;BQB=0.302156;MQ0F=0;ICB=0.0478148;HOB=0.0184911;AC=47;AN=52;DP4=19,19,226,183;MQ=60 GT:PL 1/1:255,60,0 1/1:255,51,0 1/1:255,75,0 1/1:255,48,0 1/1:255,39,0 1/1:244,30,0 1/1:255,78,0 0/1:161,0,229 0/1:247,0,150 1/1:255,42,0 0/1:177,0,135 0/1:152,0,146 1/1:255,36,0 1/1:255,48,0 1/1:255,51,0 0/1:144,0,205 1/1:255,57,0 1/1:255,54,0 1/1:215,24,0 1/1:255,51,0 1/1:255,69,0 1/1:255,63,0 1/1:255,63,0 1/1:255,57,0 1/1:255,78,0 1/1:255,51,0
使用Python代码实现字符串匹配替换:
# -*- coding: utf-8 -*-
脚本运行结果:
suhui@suhui-ThinkCentre-E74:phython_Snakmake$ python quiz1.py input1.txt input2.txt out_py.txt
suhui@suhui-ThinkCentre-E74:phython_Snakmake$ less -S out_py.txt
1 112 . C T 496 PASS DP=261;VDB=0.00422737;SGB=25.2245;RPB=0.537737;MQB=1;MQSB=1;BQB=0.570213;MQ0F=0;ICB=0.886034;HOB=0.0502959;AC=6;AN=52;DP4=189,12,25,1;MQ=60 GT:PL 0/0:0,12,129 0/0:0,30,224 0/0:0,18,141 0/0:0,3,37 0/0:0,24,221 0/0:0,27,163 0/0:0,24,188 0/1:114,0,19 0/1:114,0,180 0/0:0,24,162 1/1:129,15,0 0/1:86,0,136 0/0:0,21,195 0/0:0,18,188 0/0:0,15,123 0/1:142,0,89 0/0:0,15,149 0/0:0,24,162 0/0:0,15,129 0/0:0,18,165 0/0:0,60,209 0/0:0,36,195 0/0:0,36,190 0/0:0,33,187 0/0:0,69,208 0/0:0,30,187
1 370 . G A 999 PASS DP=480;VDB=0.765525;SGB=-9.28185;RPB=0.816783;MQB=1;MQSB=1;BQB=0.884522;MQ0F=0;ICB=0.801775;HOB=0.0650888;AC=32;AN=52;DP4=113,62,139,120;MQ=60 GT:PL 1/1:255,33,0 0/1:223,0,168 1/1:255,30,0 0/0:0,30,255 1/1:255,45,0 1/1:255,51,0 1/1:255,33,0 0/0:0,63,255 0/1:219,0,195 0/1:244,0,183 0/1:154,0,250 0/1:178,0,140 1/1:255,45,0 0/1:209,0,132 0/1:81,0,139 0/1:255,0,151 1/1:255,66,0 0/1:208,0,200 1/1:255,54,0 0/1:178,0,100 0/1:138,0,197 0/0:0,51,255 1/1:255,60,0 0/1:202,0,255 0/1:220,0,210 0/1:183,0,231
1 827 . T C 854 PASS DP=533;VDB=0.00165659;SGB=31.1247;RPB=0.465534;MQB=1;MQSB=1;BQB=0.993696;MQ0F=0;ICB=0.0478148;HOB=0.0184911;AC=5;AN=52;DP4=279,188,28,15;MQ=60 GT:PL 0/0:0,42,255 0/0:0,72,255 0/0:0,60,255 0/0:0,51,255 0/0:0,57,255 0/0:0,78,255 0/0:0,57,255 0/1:204,0,216 0/1:99,0,205 0/0:0,48,255 0/1:162,0,106 0/1:227,0,122 0/0:0,33,255 0/0:0,57,255 0/0:0,66,255 0/1:241,0,255 0/0:0,45,255 0/0:0,69,255 0/0:0,48,249 0/0:0,39,255 0/0:0,69,255 0/0:0,108,255 0/0:0,87,255 0/0:0,72,255 0/0:0,60,255 0/0:0,51,255
1 3014 . G T 999 PASS DP=454;VDB=0.958589;SGB=-9.3863;RPB=0.937155;MQB=1;MQSB=1;BQB=0.0688294;MQ0F=0;ICB=0.801775;HOB=0.0650888;AC=32;AN=52;DP4=95,70,133,124;MQ=60 GT:PL 1/1:255,51,0 0/1:120,0,228 1/1:255,51,0 0/0:0,24,203 1/1:255,66,0 1/1:255,45,0 1/1:255,33,0 0/0:0,39,255 0/1:204,0,217 0/1:255,0,230 0/1:255,0,64 0/1:110,0,128 1/1:255,39,0 0/1:220,0,242 0/1:149,0,245 0/1:207,0,145 1/1:255,42,0 0/1:152,0,152 1/1:255,51,0 0/1:166,0,130 0/1:141,0,191 0/0:0,69,255 1/1:255,48,0 0/1:171,0,173 0/1:221,0,162 0/1:192,0,177
1 3392 . G A 999 PASS DP=536;VDB=0.706175;SGB=-1.74937;RPB=0.609963;MQB=1;MQSB=1;BQB=0.292344;MQ0F=0;ICB=0.0478148;HOB=0.0184911;AC=47;AN=52;DP4=24,19,238,203;MQ=60 GT:PL 1/1:255,45,0 1/1:255,54,0 1/1:255,51,0 1/1:243,27,0 1/1:255,54,0 1/1:255,54,0 1/1:255,33,0 0/1:90,0,255 0/1:255,0,230 1/1:255,54,0 0/1:192,0,100 0/1:7,0,238 1/1:255,42,0 1/1:255,54,0 1/1:255,42,0 0/1:97,0,177 1/1:255,54,0 1/1:255,48,0 1/1:255,57,0 1/1:255,27,0 1/1:255,90,0 1/1:255,72,0 1/1:255,75,0 1/1:255,105,0 1/1:255,84,0 1/1:255,108,0
1 3430 . C T 999 PASS DP=512;VDB=0.314882;SGB=-8.61715;RPB=0.608373;MQB=1;MQSB=1;BQB=0.932661;MQ0F=0;ICB=0.0732873;HOB=0.0266272;AC=46;AN=52;DP4=28,26,212,207;MQ=60 GT:PL 1/1:255,42,0 1/1:255,48,0 1/1:255,42,0 1/1:255,30,0 1/1:255,57,0 1/1:255,63,0 1/1:255,30,0 0/1:125,0,255 0/1:250,0,239 1/1:255,54,0 0/1:157,0,131 0/1:45,0,155 1/1:255,42,0 1/1:255,48,0 1/1:255,42,0 0/1:127,0,110 1/1:255,66,0 0/1:176,0,255 1/1:255,51,0 1/1:237,30,0 1/1:255,75,0 1/1:255,72,0 1/1:255,45,0 1/1:255,111,0 1/1:255,90,0 1/1:255,102,0
1 4036 . A G 962 PASS DP=508;VDB=0.3649;SGB=62.5434;RPB=0.806806;MQB=1;MQSB=1;BQB=0.46193;MQ0F=0;ICB=0.0478148;HOB=0.0184911;AC=5;AN=52;DP4=226,197,27,24;MQ=60 GT:PL 0/0:0,57,255 0/0:0,72,255 0/0:0,39,255 0/0:0,51,255 0/0:0,48,255 0/0:0,33,255 0/0:0,51,255 0/1:198,0,120 0/1:255,0,177 0/0:0,57,255 0/1:220,0,242 0/1:114,0,212 0/0:0,63,255 0/0:0,39,255 0/0:0,54,255 0/1:255,0,104 0/0:0,54,255 0/0:0,48,255 0/0:0,30,255 0/0:0,36,255 0/0:0,60,255 0/0:0,105,255 0/0:0,90,255 0/0:0,72,255 0/0:0,57,255 0/0:0,39,255
1 4056 . C T 967 PASS DP=507;VDB=0.171323;SGB=64.0362;RPB=0.82443;MQB=1;MQSB=1;BQB=0.938503;MQ0F=0;ICB=0.0478148;HOB=0.0184911;AC=5;AN=52;DP4=226,204,32,20;MQ=60 GT:PL 0/0:0,60,255 0/0:0,60,255 0/0:0,24,255 0/0:0,51,255 0/0:0,51,255 0/0:0,39,255 0/0:0,54,255 0/1:171,0,98 0/1:244,0,196 0/0:0,60,255 0/1:231,0,214 0/1:146,0,182 0/0:0,63,255 0/0:0,45,255 0/0:0,45,255 0/1:255,0,149 0/0:0,45,255 0/0:0,42,255 0/0:0,39,255 0/0:0,33,255 0/0:0,60,255 0/0:0,108,255 0/0:0,99,255 0/0:0,75,255 0/0:0,72,255 0/0:0,51,255
1 4285 . C A 999 PASS DP=514;VDB=0.280171;SGB=59.5578;RPB=0.894116;MQB=1;MQSB=1;BQB=0.403017;MQ0F=0;ICB=0.0478148;HOB=0.0184911;AC=5;AN=52;DP4=207,231,17,32;MQ=60 GT:PL 0/0:0,45,255 0/0:0,57,255 0/0:0,36,255 0/0:0,36,255 0/0:0,45,255 0/0:0,48,255 0/0:0,39,255 0/1:217,0,59 0/1:255,0,255 0/0:0,48,255 0/1:235,0,101 0/1:176,0,205 0/0:0,75,255 0/0:0,81,255 0/0:0,66,255 0/1:227,0,161 0/0:0,51,255 0/0:0,57,255 0/0:0,45,255 0/0:0,33,255 0/0:0,72,255 0/0:0,84,255 0/0:0,90,255 0/0:0,72,255 0/0:0,75,255 0/0:0,57,255
1 4606 . T C 999 PASS DP=494;VDB=0.0629205;SGB=-17.2915;RPB=0.463604;MQB=1;MQSB=1;BQB=0.302156;MQ0F=0;ICB=0.0478148;HOB=0.0184911;AC=47;AN=52;DP4=19,19,226,183;MQ=60 GT:PL 1/1:255,60,0 1/1:255,51,0 1/1:255,75,0 1/1:255,48,0 1/1:255,39,0 1/1:244,30,0 1/1:255,78,0 0/1:161,0,229 0/1:247,0,150 1/1:255,42,0 0/1:177,0,135 0/1:152,0,146 1/1:255,36,0 1/1:255,48,0 1/1:255,51,0 0/1:144,0,205 1/1:255,57,0 1/1:255,54,0 1/1:215,24,0 1/1:255,51,0 1/1:255,69,0 1/1:255,63,0 1/1:255,63,0 1/1:255,57,0 1/1:255,78,0 1/1:255,51,0