Fixing vcf parsing bug

polio-nanopore · Oct 31, 2023 · 850740f · 850740f
1 parent 9338444
commit 850740f
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/piranha/analysis/parseHaplos.py b/piranha/analysis/parseHaplos.py
@@ -44,8 +44,8 @@ def parseVCFCalls(vcfFileList):
 
                     elif len(ref) > 1:
                         #MNP with no del, split
-                        for pos2, ref2, alt2 in zip(pos,ref,alt):
-                            haps[i].append((pos2,alt2))
+                        haps[i].append((pos,alt[0]))
+                        haps[i].append((str(int(pos)+1),alt[1]))
 
                     else:
                         #just a normal SNP