我有一个像这样的文本文件
scaffold20 6146680 . T C 44.4146 . DP=2;VDB=0.02;SGB=-0.4
scaffold20 6146696 . G A 8.13869 . DP=1;SGB=-0.379885;MQ0
scaffold20 6146760 . A G 8.13869 . DP=1;SGB=-0.379885;MQ0
scaffold20 6146785 . A G 8.13869 . DP=1;SGB=-0.379885;MQ0
scaffold20 6146864 . A C 153 . DP=7;VDB=0.637622;SGB
scaffold20 6146867 . G A 11.4845 . DP=8;VDB=0.82;SGB=-0.45
scaffold20 6146914 . G A 20.2676 . DP=5;VDB=0.06;SGB=-0.45
scaffold20 6147094 . G A 44.4146 . DP=2;VDB=0.44;SGB=-0.45
scaffold20 6147165 . C T 8.13869 . DP=1;SGB=-0.379885;MQ0F=
scaffold20 6147166 . A G 8.13869 . …Run Code Online (Sandbox Code Playgroud) 我有这样的文件(VCF)
##fileformat=VCFv4.0
##INFO=<ID=NS,Number=1,Type=Integer,Description="Number of Samples With Data">
##INFO=<ID=DP,Number=1,Type=Integer,Description="Total Depth">
##FILTER=<ID=q10,Description="Quality below 10">
##FILTER=<ID=s50,Description="Less than 50% of samples have data">
##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">
##FORMAT=<ID=GQ,Number=1,Type=Integer,Description="Genotype Quality">
##FORMAT=<ID=DP,Number=1,Type=Integer,Description="Read Depth">
##FORMAT=<ID=HQ,Number=2,Type=Integer,Description="Haplotype Quality">
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT NA00001
Chr02 259 . A . 20 . . GT:DP:A:C:G:T:PP:GQ 0/0:1:0,1:0,0:0,0:0,0:0,26,23,75,33,33,33,47,52,49:23
Chr02 260 . C . 13 . . GT:DP:A:C:G:T:PP:GQ 0/0:1:0,0:0,1:0,0:0,0:24,0,70,17,25,49,43,25,25,44:16
Chr02 261 . C . 13 . . GT:DP:A:C:G:T:PP:GQ 0/0:1:0,0:0,1:0,0:0,0:24,0,194,18,25,49,44,25,25,45:16
Chr02 262 . C A 21 . . GT:DP:A:C:G:T:PP:GQ 0/1:1:0,0:0,1:0,0:0,0:387,0,342,348,25,368,376,25,25,368:25
Chr02 263 . C …Run Code Online (Sandbox Code Playgroud)