althonos / althonos/pyrodigal

Support for circular sequences

Aperta
#67 1 commento 0 reazioni 0 assegnatari Vedi su GitHub
enhancement
Lingua principale
Cython
Stelle
189
Fork
11
Metriche di merge delle PR
Nessuna PR unita negli ultimi 30g

Descrizione

In theory, circular sequences should not contain partial genes, as genes spanning the breakpoint should appear seamlessly at both the beginning and the end of the sequence. However, Prodigal/`pyrodigal` currently do not: (1) assign the same gene ID to partial sequences at both ends of the sequence; and (2) more critically, it treats the sequence edges as independent. As a result, a partial gene can sometimes be identified at only one end.

To address this issue, I've been using [a script](https://github.com/apcamargo/reorient-circular-seq) that I wrote that iteratively changes the breakpoint to minimize gene truncation. However, this approach is obviously suboptimal and can occasionally fail (i.e. find no breakpoint that eliminates truncations), as the predicted genes may change with each breakpoint change.

Although addressing this limitation would require significant effort, it would `pyrodigal` stand out among gene prediction tools.

Guida per i contributori

Apri la guida per i contributori

Valutazione

Questa issue non è ancora stata valutata.

Ricevi le nuove issue nella tua casella

Un breve riepilogo di issue GitHub adatte ai principianti.