docs: def extract_rmlst_genes(pair, database, forward_out, reverse_out, threads=12, logfile=None):
      """
      Given a pair of reads and an rMLST database, will extract reads that contain sequence from the database.
      :param pair: List containing path to forward reads at index 0 and path to reverse reads at index 1.
      :param database: Path to rMLST database, in FASTA format.
      :param forward_out:
      :param reverse_out:
      :param threads:
      """
      out, err, cmd = bbtools.bbduk_bait(database, pair[0], forward_out, reverse_in=pair[1],
                                         reverse_out=reverse_out, threads=str(threads), returncmd=True)
      if logfile:
          write_to_logfile(logfile, out, err, cmd)