code: def get_distance(a, b, xaxis=True):
      
      if xaxis:
          arange = ("0", a.qstart, a.qstop, a.orientation)  # 0 is the dummy chromosome
          brange = ("0", b.qstart, b.qstop, b.orientation)
      else:
          arange = ("0", a.sstart, a.sstop, a.orientation)
          brange = ("0", b.sstart, b.sstop, b.orientation)
  
      dist, oo = range_distance(arange, brange, distmode="ee")
      dist = abs(dist)
  
      return dist
docs: Returns the distance between two blast HSPs.