docs: def find_paired_reads(fastq_directory, forward_id='_R1', reverse_id='_R2'):
      """
      Looks at a directory to try to find paired fastq files. Should be able to find anything fastq.
      :param fastq_directory: Complete path to directory containing fastq files.
      :param forward_id: Identifier for forward reads. Default R1.
      :param reverse_id: Identifier for reverse reads. Default R2.
      :return: List containing pairs of fastq files, in format [[forward_1, reverse_1], [forward_2, reverse_2]], etc.
      """
      pair_list = list()
  
      fastq_files = glob.glob(os.path.join(fastq_directory, '*.f*q*'))
      for name in sorted(fastq_files):
          if forward_id in name and os.path.isfile(name.replace(forward_id, reverse_id)):
              pair_list.append([name, name.replace(forward_id, reverse_id)])
      return pair_list