命令如下:
# 计算 π
vcftools --vcf ../SNP.vcf --site-pi --keep pop.txt --out pop
# 计算 Fst
vcftools --vcf ../SNP.vcf --weir-fst-pop pop1.txt --weir-fst-pop pop2.txt --fst-window-size 5000 --fst-window-step 500 --out pop1_vs_pop2
报错如下:
Warning: Expected at least 2 parts in FORMAT entry: ID=DP4,Number=4,Type=Integer,Description="Number of high-quality ref-fwd, ref-reverse, alt-fwd and alt-reverse bases">
Warning: Expected at least 2 parts in FORMAT entry: ID=DP4,Number=4,Type=Integer,Description="Number of high-quality ref-fwd, ref-reverse, alt-fwd and alt-reverse bases">
Warning: Expected at least 2 parts in INFO entry: ID=AC,Number=A,Type=Integer,Description="Allele count in genotypes for each ALT allele, in the same order as listed">
Warning: Expected at least 2 parts in INFO entry: ID=DP4,Number=4,Type=Integer,Description="Number of high-quality ref-forward , ref-reverse, alt-forward and alt-reverse bases">
Warning: Expected at least 2 parts in INFO entry: ID=DP4,Number=4,Type=Integer,Description="Number of high-quality ref-forward , ref-reverse, alt-forward and alt-reverse bases">
查找出现警告的原因:
找到的答案是:vcftools不知道如何处理Description标记内的逗号。 如果在Description中删除该逗号,警告将消失。 否则,通常可以忽略它。
结果
确实如回答所说,当把Description中的逗号去掉之后,就不会再出现这种警告了。