neherlab / pangraph

A bioinformatic toolkit to align genome assemblies into pangenome graphs
https://neherlab.github.io/pangraph
MIT License
82 stars 7 forks source link

GFAv1 output #17

Open ekg opened 2 years ago

ekg commented 2 years ago

Is it possible to obtain the paths as directed walks through the graph in GFAv1 format?

nnoll commented 2 years ago

It is not currently a feature within the code, as we only implement GFAv1 and not GFAv1.1, however I think this would be easy to add in practice. My only hangup on exposing paths as a directed walk is there is always the reverse-complementation symmetry which means the path taken in reverse (with all pancontigs traversed in the opposite orientation) is also valid.

I will assign this to myself and add the feature time permitting.

ekg commented 2 years ago

I would suggest GFAv1. v1.1 adds a second path type for metadata on paths that seems relevant haplotype panels. This breaks most existing tools for no general benefit.

On Thu, Feb 24, 2022, 18:28 Nicholas Noll @.***> wrote:

It is not currently a feature within the code, as we only implement GFAv1 and not GFAv1.1, however I think this would be easy to add in practice. My only hangup on exposing paths as a directed walk is there is always the reverse-complementation symmetry which means the path taken in reverse (with all pancontigs traversed in the opposite orientation) is also valid.

I will assign this to myself and add the feature time permitting.

— Reply to this email directly, view it on GitHub https://github.com/neherlab/pangraph/issues/17#issuecomment-1050087725, or unsubscribe https://github.com/notifications/unsubscribe-auth/AABDQEKXGTWZAT47RE7JR53U4ZTDNANCNFSM5PH7DZJQ . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

You are receiving this because you authored the thread.Message ID: @.***>