dridk / PyVCF3

A Variant Call Format reader for Python.
http://pyvcf.readthedocs.org/en/latest/index.html
Other
52 stars 9 forks source link

Header with empty VCF #5

Open antonylebechec opened 1 year ago

antonylebechec commented 1 year ago

Dear Dridk,

I develop a script to manage databases in VCF and Parquet. I need specifically to deal with VCF header, to check VCF structure (INFO fields, FORMAT...).

When I create an instance with an empty VCF, I mean only VCF header, the content of the object is empty (no variant, which is ok), but I do not have any information about the header (e.g. list of INFO fields). Apparently, only INFO field present in variant INFO column are store in the header of the object.

Is there a way to keep INFO field (and other header information) even if thy are not present in variant? Especially if the VCF is empty (no variant)?

Thank!

Best,