perl正则表达式匹配任意单字符_使用Perl和Python实现R match()字符串匹配替换

bf5ed01b8595092b2d5dac732e247528.png

文件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
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值