diff --git a/heterogenesis_varincorp.py b/heterogenesis_varincorp.py
index e5c0dcb..cac2c1d 100644
--- a/heterogenesis_varincorp.py
+++ b/heterogenesis_varincorp.py
@@ -146,7 +146,7 @@ def writevcffile(directory,prefix,clo,combvcfs,prochro):
             file.write('##INFO=<ID=NS,Number=1,Type=Integer,Description="Number of Samples With Data">\n')
             file.write('##FORMAT=<ID=AF,Number=A,Type=Float,Description="Alt allele frequency">\n')
             file.write('##FORMAT=<ID=TC,Number=1,Type=Integer,Description="Total copies of alt allele">\n')
-            file.write('##FORMAT=<ID=HS,Number=1,Type=Integer,Description="Haplotypes">\n')
+            file.write('##FORMAT=<ID=HS,Number=1,Type=Character,Description="Haplotypes">\n')
             file.write('##FORMAT=<ID=HC,Number=.,Type=Integer,Description="Total copies of alt allele per haplotype">\n')
             file.write('##FORMAT=<ID=CN,Number=2,Type=Integer,Description="Copy number at position">\n')
             file.write('#CHROM\tPOS\tID\tREF\tALT\tQUAL\tFILTER\tINFO\tFORMAT\t'+clo+'\n')
@@ -252,10 +252,26 @@ def invert(forward):
                 substitutions['g']='c'
                 substitutions['N']='N'
                 substitutions['n']='n'
-                substitutions['R']='R'
-                substitutions['r']='r'
-                substitutions['M']='M'
-                substitutions['m']='m'
+                substitutions['R']='Y'
+                substitutions['Y']='R'
+                substitutions['y']='r'
+                substitutions['r']='y'
+                substitutions['S']='S'
+                substitutions['s']='s'
+                substitutions['W']='W'
+                substitutions['w']='w'
+                substitutions['K']='M'
+                substitutions['M']='K'
+                substitutions['k']='m'
+                substitutions['m']='k'
+                substitutions['B']='V'
+                substitutions['V']='B'
+                substitutions['b']='v'
+                substitutions['v']='b'
+       	       	substitutions['D']='H'
+       	       	substitutions['H']='D'
+       	       	substitutions['d']='h'
+       	       	substitutions['h']='d'
                 for i in forward[::-1]:
                     reverse.append(substitutions[i])
                 return reverse