Header information lost when merging profile databases

[tomdstanton]

Hi Martin, Milot and team!

First, thank you for MMSeqs! I'm a long-time listener, first-time caller 😄

I have pre-defined clusters of protein sequences that I would like to store as an MMSeqs profile database for searching against, which will replace a nHMM database.

I am using a pipeline to create individual profiles per cluster:

createdb cluster.faa clusterDB
# create fake prefilter
align clusterDB clusterDB fake_prefilter alignmentDB
result2profile clusterDB clusterDB alignmentDB profile{n}DB

This works great, however when I merge the profiles, the header information is lost during merging, which is apparent as there is no mergedDB_h file. When I perform a search against the merged profile database, the convertalis command raises the "database requires header" error.

mergedb profile{1}DB mergedDB profile{2}DB profile{3}DB profile{n}DB ...

Am I merging these databases correctly? From the command help, it suggests that the first argument is the first profile database, then the name of the merged database, then all the other profile databases.

Kind regards, Tom