docs: def _define_range(self, sequences):
          '''
          define_range - define the maximum range within which two hits in a db
          search can be linked. This is defined as 1.5X the average length of all
          reads in the database.
  
          Parameters
          ----------
          sequences : str
              A path to the sequences in FASTA format. This fasta file is assumed
              to be in the correct format. i.e. headers start with '>'
  
          Returns
          -------
          max_range : int
              As described above, 1.5X the size of the average length of genes
              within the database
          '''
          sequence_count = 0
          total_sequence = 0
  
          for record in SeqIO.parse(open(sequences), 'fasta'): # For every sequence in the output
              total_sequence+=1 # increment the total sequence count
              sequence_count+=len(record.seq) # add the length of that sequence to the total string count
  
          # Get the average, multiply by 1.5, and return
          max_range = (sequence_count/total_sequence)*1.5
          return max_range