A error during annotation: KeyError: 'MER0035756'

[aloser0621]

Hi,

After setup by the command 'DRAM-setup.py prepare_databases --output_dir ~/0_db/DRAM_db', I attempted to annotate my assemblied contig, and the command is 'DRAM.py annotate -i 3_blastx_nr_1e-5_result.fasta -o annotate'. But an error occured when getting description of hits from peptidase.

Here is the report:

1 fastas found
2022-06-15 10:24:54.076054: Annotation started
0:00:00.005417: Retrieved database locations and descriptions
0:00:00.005450: Annotating 3_blastx_nr_1e-5_result_query_potential_virus
0:00:09.874450: Turning genes from prodigal to mmseqs2 db
0:00:11.801045: Getting hits from kofam
0:20:09.002767: Getting forward best hits from peptidase
0:20:14.859413: Getting reverse best hits from peptidase
0:20:15.490429: Getting descriptions of hits from peptidase
/home/disk-mounted/liaomengen/software/miniconda3/envs/DRAM/lib/python3.10/site-packages/mag_annotator/database_handler.py:81: UserWarning: No descriptions were found for your id's. Does this MER0035756 look like an id from peptidase_description
warnings.warn("No descriptions were found for your id's. Does this %s look like an id from %s" % (list(ids)[0],
Traceback (most recent call last):
File "/home/disk-mounted/liaomengen/software/miniconda3/envs/DRAM/bin/DRAM.py", line 189, in
args.func(**args_dict)
File "/home/disk-mounted/liaomengen/software/miniconda3/envs/DRAM/lib/python3.10/site-packages/mag_annotator/annotate_bins.py", line 1040, in annotate_bins_cmd
annotate_bins(list(set(fasta_locs)), output_dir, min_contig_size, prodigal_mode, trans_table, bit_score_threshold,
File "/home/disk-mounted/liaomengen/software/miniconda3/envs/DRAM/lib/python3.10/site-packages/mag_annotator/annotate_bins.py", line 1079, in annotate_bins
all_annotations = annotate_fastas(fasta_locs, output_dir, db_handler, min_contig_size, prodigal_mode, trans_table,
File "/home/disk-mounted/liaomengen/software/miniconda3/envs/DRAM/lib/python3.10/site-packages/mag_annotator/annotate_bins.py", line 1013, in annotate_fastas
annotate_fasta(fasta_loc, fasta_name, fasta_dir, db_handler, min_contig_size, prodigal_mode, trans_table,
File "/home/disk-mounted/liaomengen/software/miniconda3/envs/DRAM/lib/python3.10/site-packages/mag_annotator/annotate_bins.py", line 921, in annotate_fasta
annotations = annotate_orfs(gene_faa, db_handler, tmp_dir, start_time, custom_db_locs, custom_hmm_locs,
File "/home/disk-mounted/liaomengen/software/miniconda3/envs/DRAM/lib/python3.10/site-packages/mag_annotator/annotate_bins.py", line 821, in annotate_orfs
annotation_list.append(do_blast_style_search(query_db, db_handler.db_locs['peptidase'], tmp_dir,
File "/home/disk-mounted/liaomengen/software/miniconda3/envs/DRAM/lib/python3.10/site-packages/mag_annotator/annotate_bins.py", line 684, in do_blast_style_search
hits = formater(hits, header_dict)
File "/home/disk-mounted/liaomengen/software/miniconda3/envs/DRAM/lib/python3.10/site-packages/mag_annotator/annotate_bins.py", line 187, in get_peptidase_description
header = header_dict[peptidase_hit]
KeyError: 'MER0035756'

Could you please help me to fix this?

Looking forward to your reply

Thanks.

Best,

Moon Liao

[rmFlynn]

Hi, usually this is solved by running DRAM-setup.py update_description_db let me know if that works for you.

[aloser0621]

Thanks! It works and completed annotations after update!