-a option : can we concatenate custom output with precomputed blast against nr ? about genera HOT 2 OPEN

Dear Paul,
It is possible to attach two diamond outputs together using cat and feed them to genEra. The only problem with your approach is that extra_Diamond_results.bout does not contain the fifth column with the NCBI taxonomy IDs of the extra proteomes. This column is essential for the gene age inferences performed in step 3. All of this is done automatically with -a, but genEra needs to run against the NR before integrating the extra proteomes into the analysis. Maybe I can add an argument to skip the search against the NR and run step 1 using only the extra proteomes. Then you could manually attach both diamond results together with cat. Would that implementation be useful to you?
Best,
Josué.

from genera.

Dear Josué,

Yes it would be !

Thanks for your work,

Paul

from genera.