Hi!
I noticed that the output from get_operons.py is not what I would expect, at least not in all cases.
Something can go wrong when parsing gene-gene relationship. For example, the output that I get contains wrong order of genes, which messes up operon topology (creates long, overlapping operons). For example:
...
7,ID=gene-ABHI20_RS00215
7,ID=gene-ABHI20_RS00220
7,ID=gene-ABHI20_RS00225
7,ID=gene-ABHI20_RS00230
7,ID=gene-ABHI20_RS00235
8,ID=gene-ABHI20_RS00240
9,ID=gene-ABHI20_RS0025050
10,ID=gene-ABHI20_RS0931020
11,ID=gene-ABHI20_RS0025555
12,ID=gene-ABHI20_RS0026060
Please note operon 10 (single gene) with an ID completely out of range and far from the other surrounding genes.
Hi!
I noticed that the output from
get_operons.pyis not what I would expect, at least not in all cases.Something can go wrong when parsing gene-gene relationship. For example, the output that I get contains wrong order of genes, which messes up operon topology (creates long, overlapping operons). For example:
Please note operon 10 (single gene) with an ID completely out of range and far from the other surrounding genes.