From 0ac3e639397843e365c35df2cb082240b605a504 Mon Sep 17 00:00:00 2001
From: Jonathan Dursi <jonathan@dursi.ca>
Date: Tue, 13 Dec 2016 15:51:09 -0500
Subject: [PATCH] add lowsupport filter to header

---
 mergevcf/mergedfile.py | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/mergevcf/mergedfile.py b/mergevcf/mergedfile.py
index 4bda646..ff4b03f 100644
--- a/mergevcf/mergedfile.py
+++ b/mergevcf/mergedfile.py
@@ -123,7 +123,9 @@ def infoString(callers, infodict):
     outfile.write('##INFO=<ID=Callers,Number=.,Type=String,Description="Callers that made this call">\n')
     if output_ncallers:
         outfile.write('##INFO=<ID=NumCallers,Number=1,Type=Integer,Description="Number of callers that made this call">\n')
-    outfile.write("#CHROM  POS     ID      REF     ALT     QUAL    FILTER  INFO\n")
+    if min_num_callers > 0:
+        outfile.write('##FILTER=<ID=LOWSUPPORT,Description="Not called by enough callers in ensemble">\n')
+    outfile.write("#CHROM\tPOS\tID\tREF\tALT\tQUAL\tFILTER\tINFO\n")
 
     for variant in calldict:
         callers = variant[2]