sanger-pathogens/ariba

Summary ignores one gene

GonzaloYebra opened this issue · 0 comments

First of all, thanks to the authors for yet another great tool.

I've been running Ariba to infer presence/absence of 4 genes of interest. It runs without problems but when I run summary, the output only shows 3 yes/no columns.

From the individual reports (in which the gene seems to be analysed without problems) I can't see any pattern in the metrics that would justify ignoring the gene completely.

I've tried to run the gene "ignored" on its own and Ariba again works but in this case when running summary I get the following error:

ariba.summary.Error: Cannot calculate distance matrix to make tree for phandango. Not enough columns

I'm not sure if this is relevant, but the gene in question is always labelled as "cluster", the rest being "cluster_1", "cluster_2" and "cluster_3". This is irrespective of the order in the fasta file, I guess they are ordered by length (and indeed the problematic gene is the longest with 1,158bp).

I'd appreciate any help with this. Thanks a lot!

Gonzalo