VCF allele frequency is not sensible

Question

[wrab425]

I have a vcf file containing the variants from 64 haploid strains. I need to be able to get the allele frequencies for further work using plink but when I run the vcftools command

vcftools --vcf merged_SNPs_only.vcf --freq --out freq_analysis

I get this, The problem is I have variants that are in more than one strain but the output suggests that the frequency of the allele is 0 or 1

What is wrong with my vcf or my command?

CHROM   POS     N_ALLELES       N_CHR   {ALLELE:FREQ}
chr1    4446    2       1       C:0     T:1
chr1    4570    2       1       T:0     C:1
chr1    4575    2       1       T:0     C:1
chr1    4595    2       1       A:0     C:1
chr1    4676    2       1       A:0     G:1
chr1    4726    2       2       T:0     C:1
chr1    4865    2       7       A:0     G:1
chr1    4870    2       7       T:0     C:1
chr1    4883    2       7       A:0     T:1
chr1    4952    2       4       G:0     T:1
chr1    5015    2       4       G:0     A:1
chr1    5098    2       7       A:0     C:1
chr1    5168    2       6       G:0     A:1
chr1    5252    2       7       A:0     T:1
chr1    5253    2       4       T:0     A:1
chr1    5265    2       7       C:0     T:1
chr1    5270    2       7       G:0     T:1
chr1    5302    2       2       T:0     G:1

[Ram]

Does your VCF file have all 64 samples in it? Also, the plink tag is not relevant here as plink is not involved in your current situation.

[wrab425]

Thanks Ram, my file does have all 64 samples in it. I put the plink in there because I get the same result when I look at allele frequencies using plink. so the problem is illustrated with a snp like chr1 5265 2 7 C:0 T:1 which has seven chromosomes out of the 64 in the set but the T allele frequency is 1 when by my understanding should be 7/64.

[Ram]

What's your vcftools version?

[wrab425]

Hi Ram the version: 0.1.15

[Ram]

Yeah, the problem lies with your input file. Please follow my instructions on the gist excerpt so we can understand this better.

[wrab425]

super thanks Ram, I have entered the command and put it into a gist but if it would help to see more of the file please ask. We seem to be about to make progress

[Ram]

It would help if we had a few lines where you expect to see a different output than what you see. For example, chr1:5265. I'm not sure if my grep was badly written (I didn't test it), but it's supposed to have produced 10 lines. Remove the -m1 and replace the -A11 with -C11, and put that in the gist if you can please? The output of that along with the line for chr1:5265 should be really helpful.

[Ram]

Actually, if plink is giving you the same output, the problem probably lies with the data being piped in to the tools. Are you working on the file directly or is any other tool involved? Can you paste an except of your VCF in a GitHub gist and paste the link to the gist in your post please? Run this to get the excerpt:

grep -m1 -A11 "^#CHROM" vcf_file >file_to_be_pasted_to_a_github_gist

[wrab425]

chr1    4446    .   C   T   128 PASS    AC=1;AF=1.00;AN=1;BaseQRankSum=-0.566;ClippingRankSum=-0.382;DP=11;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=26.01;MQRankSum=0.978;QD=11.64;ReadPosRankSum=0.000;SF=2;SOR=1.802    GT:GQ:DP:PL:AD  .:.:.:.:.   .:.:.:.:.   1:99:11:.,.,.:3,8   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.   .:.:.:.:.

and

chr1 4446 . C T 128 PASS AC=1;AF=1.00;AN=1;BaseQRankSum=-0.566;ClippingRankSum=-0.382;DP=11;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=26.01;MQRankSum=0.978;QD=11.64;ReadPosRankSum=0.000;SF=2;SOR=1.802 GT:GQ:DP:PL:AD .:.:.:.:. .:.:.:.:. 1:99:11:.,.,.:3,8 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.

view raw gistfile1.txt hosted with ❤ by GitHub

[Ram]

This agrees with the frequency, in that only one sample has a variant at this locus.

[wrab425]

Hi Ram,

We do appreciate your advice as we are stuck with this. The thing is that there are snps that are present in 30 or so of the samples and they also give a MAF of 0 so this seems wrong. If you look at the second gitgist there are some snps in that sample which are present in seven samples and of course they have a MAF of 0. I think it must be because the vcf is haploid but I am not sure hw to fix it.

Thanks,

William

[Ram]

Your gist has only one row(chr1:4446) in it. I don't see a second gist.

EDIT: OK, you haven't linked it here. The gist is

	chr1 4446 . C T 128 PASS AC=1;AF=1.00;AN=1;BaseQRankSum=-0.566;ClippingRankSum=-0.382;DP=11;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=26.01;MQRankSum=0.978;QD=11.64;ReadPosRankSum=0.000;SF=2;SOR=1.802 GT:GQ:DP:PL:AD .:.:.:.:. .:.:.:.:. 1:99:11:.,.,.:3,8 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 4570 . T C 420 PASS AC=1;AF=1.00;AN=1;BaseQRankSum=-1.646;ClippingRankSum=0.690;DP=18;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=25.16;MQRankSum=-0.903;QD=23.33;ReadPosRankSum=-1.434;SF=6;SOR=2.303 GT:GQ:DP:PL:AD .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. 1:99:18:.,.,.:4,14 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 4575 . T C 420 PASS AC=1;AF=1.00;AN=1;BaseQRankSum=-1.327;ClippingRankSum=-0.053;DP=18;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=25.16;MQRankSum=-0.584;QD=23.33;ReadPosRankSum=-1.540;SF=6;SOR=2.303 GT:GQ:DP:PL:AD .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. 1:99:18:.,.,.:4,14 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 4595 . A C 689 PASS AC=1;AF=1.00;AN=1;BaseQRankSum=-0.296;ClippingRankSum=-0.619;DP=29;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=25.29;MQRankSum=-1.050;QD=23.76;ReadPosRankSum=-0.188;SF=6;SOR=2.546 GT:GQ:DP:PL:AD .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. 1:99:29:.,.,.:6,23 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 4676 . A G 393 PASS AC=1;AF=1.00;AN=1;DP=10;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=27.00;QD=29.84;SF=0;SOR=4.804 GT:GQ:DP:PL:AD 1:99:10:.,.,.:0,10 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 4726 . T C 234 PASS AC=2;AF=1.00;AN=2;BaseQRankSum=0.665;ClippingRankSum=1.096;DP=79;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=27.00;MQRankSum=0.505;QD=29.00;ReadPosRankSum=-0.197;SF=0,4;SOR=4.977 GT:GQ:DP:PL:AD 1:99:11:.,.,.:0,11 .:.:.:.:. .:.:.:.:. .:.:.:.:. 1:99:68:.,.,.:31,37 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 4865 . A G 5036.86 PASS AC=7;AF=1.00;AN=7;DP=789;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=27.64;QD=31.32;SF=0,1,2,3,4,5,6;SOR=6.518 GT:GQ:DP:PL:AD 1:99:25:.,.,.:0,25 1:99:236:.,.,.:0,236 1:99:186:.,.,.:0,186 1:99:57:.,.,.:0,57 1:99:145:.,.,.:0,145 1:99:21:.,.,.:0,21 1:99:119:.,.,.:0,119 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 4870 . T C 5146.71 PASS AC=7;AF=1.00;AN=7;DP=802;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=27.61;QD=35.16;SF=0,1,2,3,4,5,6;SOR=6.593 GT:GQ:DP:PL:AD 1:99:26:.,.,.:0,26 1:99:241:.,.,.:0,241 1:99:187:.,.,.:0,187 1:99:57:.,.,.:0,57 1:99:151:.,.,.:0,151 1:99:21:.,.,.:0,21 1:99:119:.,.,.:0,119 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 4883 . A T 5184.14 PASS AC=7;AF=1.00;AN=7;DP=819;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=27.57;QD=28.75;SF=0,1,2,3,4,5,6;SOR=6.736 GT:GQ:DP:PL:AD 1:99:28:.,.,.:0,28 1:99:247:.,.,.:0,247 1:99:192:.,.,.:0,192 1:99:57:.,.,.:0,57 1:99:151:.,.,.:0,151 1:99:21:.,.,.:0,21 1:99:123:.,.,.:0,123 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 4952 . G T 4142.47 PASS AC=4;AF=1.00;AN=4;BaseQRankSum=-0.196;ClippingRankSum=1.134;DP=891;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=26.83;MQRankSum=1.688;QD=20.67;ReadPosRankSum=-0.711;SF=1,2,3,4;SOR=9.327 GT:GQ:DP:PL:AD .:.:.:.:. 1:99:426:.,.,.:109,317 1:99:213:.,.,.:68,145 1:99:61:.,.,.:9,52 1:99:182:.,.,.:85,97 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 5015 . G A 4508.5 PASS AC=4;AF=1.00;AN=4;BaseQRankSum=-0.109;ClippingRankSum=0.855;DP=1438;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=26.79;MQRankSum=3.011;QD=14.03;ReadPosRankSum=0.415;SF=1,2,3,4;SOR=2.076 GT:GQ:DP:PL:AD .:.:.:.:. 1:99:919:.,.,.:252,667 1:99:246:.,.,.:78,168 1:99:75:.,.,.:16,59 1:99:195:.,.,.:79,116 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 5098 . A C 6470.14 PASS AC=7;AF=1.00;AN=7;BaseQRankSum=1.643;ClippingRankSum=2.140;DP=1579;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=27.00;MQRankSum=1.566;QD=30.00;ReadPosRankSum=1.964;SF=0,1,2,3,4,5,6;SOR=7.701 GT:GQ:DP:PL:AD 1:99:46:.,.,.:0,46 1:99:911:.,.,.:3,908 1:99:168:.,.,.:0,168 1:99:55:.,.,.:0,55 1:99:142:.,.,.:0,142 1:99:45:.,.,.:0,45 1:99:212:.,.,.:49,163 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 5168 . G A 13764.2 PASS AC=6;AF=1.00;AN=6;BaseQRankSum=0.090;ClippingRankSum=-1.406;DP=2363;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=27.41;MQRankSum=0.073;QD=31.13;ReadPosRankSum=0.398;SF=0,1,2,3,4,5;SOR=6.229 GT:GQ:DP:PL:AD 1:99:45:.,.,.:0,45 1:99:1638:.,.,.:0,1638 1:99:287:.,.,.:2,285 1:99:83:.,.,.:0,83 1:99:262:.,.,.:0,262 1:99:46:.,.,.:0,46 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 5252 . A T 19479.4 PASS AC=7;AF=1.00;AN=7;BaseQRankSum=1.059;ClippingRankSum=-1.382;DP=3410;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=27.49;MQRankSum=0.709;QD=34.60;ReadPosRankSum=1.041;SF=0,1,2,3,4,5,6;SOR=1.240 GT:GQ:DP:PL:AD 1:99:75:.,.,.:0,75 1:99:2107:.,.,.:0,2107 1:99:386:.,.,.:1,385 1:99:114:.,.,.:0,114 1:99:348:.,.,.:0,348 1:99:68:.,.,.:0,68 1:99:310:.,.,.:151,159 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 5253 . T A 9879.5 PASS AC=4;AF=1.00;AN=4;BaseQRankSum=3.701;ClippingRankSum=0.169;DP=2971;FS=8.722;MLEAC=1;MLEAF=1.00;MQ=32.73;MQRankSum=-4.009;QD=15.76;ReadPosRankSum=-2.961;SF=1,2,3,4;SOR=0.290 GT:GQ:DP:PL:AD .:.:.:.:. 1:99:2119:.,.,.:517,1602 1:99:388:.,.,.:126,262 1:99:114:.,.,.:31,83 1:99:347:.,.,.:161,186 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 5265 . C T 22069.9 PASS AC=7;AF=1.00;AN=7;DP=3365;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=27.55;QD=30.88;SF=0,1,2,3,4,5,6;SOR=2.338 GT:GQ:DP:PL:AD 1:99:67:.,.,.:0,67 1:99:2067:.,.,.:0,2067 1:99:386:.,.,.:0,386 1:99:114:.,.,.:0,114 1:99:347:.,.,.:0,347 1:99:69:.,.,.:0,69 1:99:315:.,.,.:0,315 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 5270 . G T 21510.3 PASS AC=7;AF=1.00;AN=7;DP=3367;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=27.55;QD=29.13;SF=0,1,2,3,4,5,6;SOR=2.483 GT:GQ:DP:PL:AD 1:99:67:.,.,.:0,67 1:99:2067:.,.,.:0,2067 1:99:387:.,.,.:0,387 1:99:111:.,.,.:0,111 1:99:347:.,.,.:0,347 1:99:67:.,.,.:0,67 1:99:321:.,.,.:0,321 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 5302 . T G 1855 PASS AC=2;AF=1.00;AN=2;BaseQRankSum=2.614;ClippingRankSum=0.790;DP=165;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=27.87;MQRankSum=0.662;QD=29.58;ReadPosRankSum=0.987;SF=0,5;SOR=2.610 GT:GQ:DP:PL:AD 1:99:89:.,.,.:0,89 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. 1:99:76:.,.,.:21,55 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 5336 . C T 13009.8 PASS AC=4;AF=1.00;AN=4;BaseQRankSum=-1.934;ClippingRankSum=1.255;DP=3557;FS=4.538;MLEAC=1;MLEAF=1.00;MQ=34.36;MQRankSum=0.850;QD=14.07;ReadPosRankSum=21.397;SF=1,2,3,4;SOR=0.389 GT:GQ:DP:PL:AD .:.:.:.:. 1:99:2529:.,.,.:506,2023 1:99:472:.,.,.:128,344 1:99:146:.,.,.:33,113 1:99:408:.,.,.:169,239 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 5361 . A G 7442 PASS AC=3;AF=1.00;AN=3;BaseQRankSum=0.163;ClippingRankSum=0.252;DP=550;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=28.22;MQRankSum=-0.610;QD=36.27;ReadPosRankSum=4.532;SF=0,5,6;SOR=2.857 GT:GQ:DP:PL:AD 1:99:99:.,.,.:0,99 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. 1:99:70:.,.,.:0,70 1:99:366:.,.,.:19,347 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 5377 . G A 8048.33 PASS AC=3;AF=1.00;AN=3;BaseQRankSum=0.870;ClippingRankSum=0.064;DP=565;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=28.08;MQRankSum=-1.414;QD=30.38;ReadPosRankSum=5.245;SF=0,5,6;SOR=3.142 GT:GQ:DP:PL:AD 1:99:95:.,.,.:0,95 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. 1:99:109:.,.,.:0,109 1:99:361:.,.,.:13,348 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 5468 . C T 5674 PASS AC=1;AF=1.00;AN=1;BaseQRankSum=-0.325;ClippingRankSum=0.369;DP=229;FS=6.550;MLEAC=1;MLEAF=1.00;MQ=26.66;MQRankSum=-0.374;QD=24.78;ReadPosRankSum=-3.785;SF=6;SOR=0.142 GT:GQ:DP:PL:AD .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. 1:99:229:.,.,.:12,217 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 5519 . T G 9194.43 PASS AC=7;AF=1.00;AN=7;DP=1653;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=27.97;QD=28.50;SF=0,1,2,3,4,5,6;SOR=2.531 GT:GQ:DP:PL:AD 1:99:42:.,.,.:0,42 1:99:1024:.,.,.:0,1024 1:99:173:.,.,.:0,173 1:99:51:.,.,.:0,51 1:99:149:.,.,.:0,149 1:99:41:.,.,.:0,41 1:99:168:.,.,.:0,168 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 5534 . A G 5614.33 PASS AC=3;AF=1.00;AN=3;BaseQRankSum=1.581;ClippingRankSum=0.632;DP=1105;FS=3.924;MLEAC=1;MLEAF=1.00;MQ=31.26;MQRankSum=9.404;QD=15.85;ReadPosRankSum=0.221;SF=1,2,3;SOR=1.043 GT:GQ:DP:PL:AD .:.:.:.:. 1:99:902:.,.,.:230,672 1:99:159:.,.,.:51,108 1:99:42:.,.,.:11,31 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.
	chr1 5950 . C G 1384 PASS AC=1;AF=1.00;AN=1;DP=34;FS=0.000;MLEAC=1;MLEAF=1.00;MQ=39.22;QD=35.28;SF=0;SOR=7.112 GT:GQ:DP:PL:AD 1:99:34:.,.,.:0,34 .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:. .:.:.:.:.

view raw gistfile1.txt hosted with ❤ by GitHub

It's weird - I guess the haploid explanation fits - I've never seen VCF with a single value in GT instead of a x/y or x|y.