Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Feature Request: Orient a GFA #4

Open
mir-cat opened this issue Aug 14, 2023 · 1 comment
Open

Feature Request: Orient a GFA #4

mir-cat opened this issue Aug 14, 2023 · 1 comment

Comments

@mir-cat
Copy link

mir-cat commented Aug 14, 2023

I have a request related to the orientation of the sequence in GFA graphs. With the GFA graphs I am working with, the strand represented in the GFA segment sequence is essentially random. One way this manifests in the graph is that links from the same node end arrive have opposite "to orientation" For example, for nodes utig4-1529 to utig4-1530 and utig4-1527 in this fork structure

Screenshot 2023-08-14 at 3 50 36 PM

The edge from 1529 to 1530 has a + "to orientation" while the edge to 1527 has a - "to orientation"

Screenshot 2023-08-14 at 4 05 08 PM

I believe that this could be done, in the acyclic case, by traversing the graph from one end to the other and taking reverse complements when relative misorientations are detected.

@tobiasmarschall
Copy link
Member

I encouraged @mir-cat to post this feature request because I thought that this might fit into the order_gfa subcommand (see https://github.com/marschall-lab/gaftools/blob/main/gaftools/cli/order_gfa.py).

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants