The +munge plugin appears to fail to read the last column of a tab-delimited header line if it starts with a # character.
So the following (tab-delimited, b37 positions) issued a warning about not finding the p-value column:
#CHROM POS REF ALT BETA SE P
13 19141990 C G 0.5 0.001 0.05
13 19144814 T C 0.5 0.03 0.05
13 19147641 T C 0.8 0.02 0.04
13 19149280 C T -1.1 0.001 1.0e-8
If it starts with CHROM instead it will work fine.
If I add a nonsense column at the end, it will work fine.
If I move e.g. the ALT column to the end, it will error out due to not finding the alternative allele column.
Gemini suggests it could be something to do with the HTSlib parser but I haven't been able to verify.
The
+mungeplugin appears to fail to read the last column of a tab-delimited header line if it starts with a#character.So the following (tab-delimited, b37 positions) issued a warning about not finding the p-value column:
If it starts with
CHROMinstead it will work fine.If I add a nonsense column at the end, it will work fine.
If I move e.g. the
ALTcolumn to the end, it will error out due to not finding the alternative allele column.Gemini suggests it could be something to do with the HTSlib parser but I haven't been able to verify.