diff --git a/scripts/combine_vcf_variants.py b/scripts/combine_vcf_variants.py
index 511705c..c6da7f0 100755
--- a/scripts/combine_vcf_variants.py
+++ b/scripts/combine_vcf_variants.py
@@ -173,11 +173,13 @@ def has_md_tag(bam):
             hap = ''.join(alts)
             if 'DP4' in v.info:
                 haplotypes_by_strand[hap] = int(v.info['DP4'][2])+int(v.info['DP4'][3])
+            elif 'SAF' in v.info:
+                haplotypes_by_strand[hap] = int(v.info['SAF'][0])
+                haplotypes_by_strand[hap.lower()] = int(v.info['SAR'][0])
             elif 'AD' in v.samples[0]:
                 haplotypes_by_strand[hap] = int(v.samples[0]['AD'][1])
 
 
-
     haplotypes = Counter()
     for h,count in haplotypes_by_strand.items():
         haplotypes[h.upper()] += count