Skip to content

Commit 01bf864

Browse files
committed
fixed bug in retrieving format index
1 parent 06d6ab6 commit 01bf864

1 file changed

Lines changed: 1 addition & 1 deletion

File tree

vawk

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -165,7 +165,7 @@ def vawk(debug, header, var, vawk_string, vcf_file):
165165
else:
166166
# special case for wildcard (all samples)
167167
if s_key[0] == '*':
168-
perline = 'ALLSAMPLES_' + s_key[1] + '=""; for(i=1;i<length(fmt);++i) { if (fmt[i]=="' + s_key[1] + '") { fmt_key=i; break } }; for (COL=10;COL<=NF-1;++COL) { split($COL,samp,":"); ALLSAMPLES_' + s_key[1] + '=ALLSAMPLES_' + s_key[1] + '""samp[fmt_key]"\\t" }; split($NF,samp,":"); ALLSAMPLES_' + s_key[1] + '=ALLSAMPLES_' + s_key[1] + '""samp[fmt_key]; ' + perline
168+
perline = 'ALLSAMPLES_' + s_key[1] + '=""; for(i=1;i<=length(fmt);++i) { if (fmt[i]=="' + s_key[1] + '") { fmt_key=i; break } }; for (COL=10;COL<=NF-1;++COL) { split($COL,samp,":"); ALLSAMPLES_' + s_key[1] + '=ALLSAMPLES_' + s_key[1] + '""samp[fmt_key]"\\t" }; split($NF,samp,":"); ALLSAMPLES_' + s_key[1] + '=ALLSAMPLES_' + s_key[1] + '""samp[fmt_key]; ' + perline
169169
else:
170170
perline = 'split($SAMPLE["' + s_key[0] + '"],samp,":"); for(i=1;i<=length(fmt);++i) { if (fmt[i]=="' + s_key[1] + '") { SAMPLE_' + '_'.join(s_key) + '=samp[i]; break } } ' + perline
171171

0 commit comments

Comments
 (0)