Error calculating abundances in newly created DB

[ArnauVich]

Dear developers,

First of all, thank you for creating this tool :)

I used mOTUs-extender to add 106 genomes to mOTU v3.0.1. Database creation went smoothly without errors or warning messages. However, mOTUs return an error in the relative abundance calculation steps when using the new database to profile samples (see below). I believe it has something to do with a missing file or an error in the formatting of the file with genes length, but I cannot figure it out.

Do you have any idea what could be going wrong here? I would greatly appreciate any suggestions or advice you can give me.

Error message:

File "miniconda3/envs/motu-env/lib/python3.8/site-packages/motus/map_mOTUs_to_LGs.py", line 235, in
counts_mOTUs_j = [counts_mOTUs[x] for x in genes_list] # vector that represents the read counts of the mOTUs of the LGs
KeyError: 'unassigned.COG0552-23612NEWDB_1'

[MiguelJMorcuende]

I am having the same error

[hjruscheweyh]

Dear Both

I apologize for the very delayed answer. The issue might be a problem with the newline character in some OS. In fact this issue is generated by a faulty execution of the this line.

The quick fix is to go into the db_mOTUs folder and update one file = db_mOTU_MAP_MGCs_to_mOTUs_in-line.tsv.

You open this file in your text ediotor and add a newline character right after unassigned.COG0552-23612, save the file and rerun motus.

Let me know if this worked for you.

Best,
Hans