docs: def setup_allelespecific_database(fasta_file, database_folder, allele_list):
      """
      Since some genera have some rMLST genes missing, or two copies of some genes, genus-specific databases are needed.
      This will take only the alleles known to be part of each genus and write them to a genus-specific file.
      :param database_folder: Path to folder where rMLST_combined is stored.
      :param fasta_file: Path to fasta file to write allelespecific database to.
      :param allele_list: allele list generated by find_genusspecific_allele_list
      """
      index = SeqIO.index(os.path.join(database_folder, 'rMLST_combined.fasta'), 'fasta')
      seqs = list()
      for s in allele_list:
          try:
              seqs.append(index[s])
          except KeyError:
              logging.warning('Tried to add {} to allele-specific database, but could not find it.'.format(s))
      SeqIO.write(seqs, fasta_file, 'fasta')