From ffca2202b8b9489a01209f80d145d6b3207e2e95 Mon Sep 17 00:00:00 2001
From: pjentsch <pjentsch@uwaterloo.ca>
Date: Fri, 20 May 2022 16:24:13 -0700
Subject: [PATCH] add FORMAT to header as per VCFv4.2 spec

---
 msa2vcf/msa2vcf.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/msa2vcf/msa2vcf.py b/msa2vcf/msa2vcf.py
index aa7df9d..31d2705 100644
--- a/msa2vcf/msa2vcf.py
+++ b/msa2vcf/msa2vcf.py
@@ -174,7 +174,7 @@ def make_vcf(snps, qname, rname, keep_n):
     vcflines.append("##contig=<ID=" + rname + ">")
     vcflines.append("##INFO=<ID=DP,Number=1,Type=Integer,Description=\"Total Depth\">")
     vcflines.append("##INFO=<ID=AF,Number=A,Type=Float,Description=\"Allele Frequency\">")
-    vcflines.append("#CHROM\tPOS\tID\tREF\tALT\tQUAL\tFILTER\tINFO\t"+qname)
+    vcflines.append("#CHROM\tPOS\tID\tREF\tALT\tQUAL\tFILTER\tINFO\tFORMAT\t"+qname)
 
     # variants
     for line in snps: