code: def find_databases(databases):
      
      # 16 ribosomal proteins in their expected order
      proteins = ['L15', 'L18', 'L6', 'S8', 'L5', 'L24', 'L14',
              'S17', 'L16', 'S3', 'L22', 'S19', 'L2', 'L4', 'L3', 'S10']
      # curated databases
      protein_databases = {
                  'L14': 'rpL14_JGI_MDM.filtered.faa',
                  'L15': 'rpL15_JGI_MDM.filtered.faa',
                  'L16': 'rpL16_JGI_MDM.filtered.faa',
                  'L18': 'rpL18_JGI_MDM.filtered.faa',
                  'L22': 'rpL22_JGI_MDM.filtered.faa',
                  'L24': 'rpL24_JGI_MDM.filtered.faa',
                  'L2': 'rpL2_JGI_MDM.filtered.faa',
                  'L3': 'rpL3_JGI_MDM.filtered.faa',
                  'L4': 'rpL4_JGI_MDM.filtered.faa',
                  'L5': 'rpL5_JGI_MDM.filtered.faa',
                  'L6': 'rpL6_JGI_MDM.filtered.faa',
                  'S10': 'rpS10_JGI_MDM.filtered.faa',
                  'S17': 'rpS17_JGI_MDM.filtered.faa',
                  'S19': 'rpS19_JGI_MDM.filtered.faa',
                  'S3': 'rpS3_JGI_MDM.filtered.faa',
                  'S8': 'rpS8_JGI_MDM.filtered.faa'}
      protein_databases = {key: '%s/%s' % (databases, database) \
              for key, database in list(protein_databases.items())}
      return proteins, protein_databases
docs: define ribosomal proteins and location of curated databases