diff --git a/codes/filter_snpindel.pl b/codes/filter_snpindel.pl index 7850d22..3405527 100755 --- a/codes/filter_snpindel.pl +++ b/codes/filter_snpindel.pl @@ -59,7 +59,15 @@ while () { my $gene = (split(";", $line[6]))[0]; if ($pipeline eq 'germline') { - $freq = (split(";", $line[-1]))[6]; + # GT:DP:VD:AD:AF:RD:ALD + # GT:DP:VD:ALD:RD:AD:AF:BIAS:PMEAN:PSTD:QUAL:QSTD:SBF:ODDRATIO:MQ:SN:HIAF:ADJAF:NM + # 信息提取 + my @content_name = split(/:/, $line[103]); + my @content = split(/:/, $line[-1]); + # 将数据与表头对应 + my %record_content; + @record_content{@content_name} = @content; + $freq = $record_content{'AF'}; } my @reason;