modified: svdb/merge_vcf_module.py

	modified:   svdb/merge_vcf_module_cython.py

Author: jesper eisfeldt

svdb/merge_vcf_module.py

@@ -108,7 +108,7 @@ def print_header(vcf_list, vcf_dictionary, args, command_line):
     for entry in sorted(subheader):
         print(subheader[entry].strip())
     if not args.notag:
-        print("##INFO=<ID=VARID,Number=1,Type=String,Description=\"The variant ID of merged samples\">")
+        print("##INFO=<ID=FOUNDBY,Number=1,Type=Integer,Description=\"The number of files containing the variant\">")
         print("##INFO=<ID=set,Number=1,Type=String,Description=\"Source VCF for the merged record in SVDB\">")
         print("##INFO=<ID=svdb_origin,Number=1,Type=String,Description=\"pipe separated list of the VCF for the merged record in SVDB\">")
 

svdb/merge_vcf_module_cython.py

@@ -331,6 +331,7 @@ def merge(variants, samples, sample_order, sample_print_order, priority_order, a
             if not args.notag:
                 set_tag = determine_set_tag(priority_order, files)
                 line[7] += ";set={}".format(set_tag)
+                line[7] += ";FOUNDBY={}".format( len(set(filters_tag.keys())) )
 
             #add chrom information of all merged variants
             callers=[]