added gap primer output

shesketh · Sep 20, 2016 · ac01c37 · ac01c37
1 parent 058cee8
commit ac01c37
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/zippy/zippy.py b/zippy/zippy.py
@@ -397,6 +397,10 @@ def getPrimers(intervals, db, design, config, deep=True, rename=None, compatible
             primerTable.append([unquote(intervals[1].name)] + str(rankedPairs[0][3]).split('\t'))
             primerVariants[rankedPairs[0][2]].append(intervals[0])
             primerVariants[rankedPairs[0][3]].append(intervals[1])
+            # add gap interval
+            gapInterval = Interval(intervals[0].chrom, intervals[0].chromStart, intervals[1].chromEnd)
+            gapPrimerPair = PrimerPair([rankedPairs[0][2][0], rankedPairs[0][3][1]], name='gap_'+gapInterval.name )
+            primerTable.append([unquote(gapInterval.name)] + str(gapPrimerPair).split('\t'))
         else:
             missedIntervals = intervals
     else: