add lowsupport filter to header

ljdursi · Dec 13, 2016 · 0ac3e63 · 0ac3e63
1 parent 049ec33
commit 0ac3e63
Showing 1 changed file with 3 additions and 1 deletion.
diff --git 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]