Problem with -gf and -cm parameteres about clinker HOT 3 OPEN

I was trying to figure out how to use the -gf feature so I looked at your files as examples.

I don't know if this is helpful, but the first genbank file you provided (rep_H153520460_plas_col.txt) doesn't actually have any instances of NCLOOJ_* in it - is that perhaps why some genes don't get put into the groups you defined?

Edit:
Ah nevermind I've figured out how this works and see why only one of your genbank files match the feature table :)) Disregard this

from clinker.

Hello, I had the same problem.
I have 24 groups and when i run it without -gf it assigns different colors to the corresponding groups. When I run the code with -gf and the .csv, only 5 groups are colored and it is assigning those functions incorrectly to genes that do not have those functions.

code:
clinker *gb -p -gf output4.csv

without -gf:

with the csv

csv file:
output4.csv

Can someone help me understand why they are not assigning all the functions to their respective genes?

Thanks

from clinker.

Unfortunately I think you are both hitting some weirdness with how the genes get grouped and I don't have a good resolution for it. Currently they are computed by finding disjoint sets of genes based on the gene-gene alignments - if there is an alignment between genes in two different defined function groups, they will be merged into a single group leading to the incorrect function/colour. You could try raising the identity/similarity cutoffs to avoid low quality alignments which lead to bad merges (e.g. the correct match may have 100% identity, but some gene from another group has a ~30% identity match). Alternatively, if you specify the function for EVERY gene you are looking at and use the --no_align flag (i.e. no alignments are performed), the colours should work but the gene-gene links will be lost.

from clinker.