Commits

Bob Harris committed f3a8a62

fixed problem that occurred when the reference strand was negative

  • Participants
  • Parent commits 4db1702

Comments (0)

Files changed (1)

File scripts/maf_extract_ranges.py

 
     for maf in bx.align.maf.Reader( sys.stdin ):
         ref_component = maf.components[ refindex ]
+        if ref_component.strand == '-':
+            maf = block.reverse_complement()
+            ref_component = maf.components[ refindex ]
         # Find overlap with reference component
         intersections = intersecter.find( ref_component.start, ref_component.end )
         # Keep output maf ordered