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

How should we interpret SR of the L line in rGFA #19

Open
SHuang-Broad opened this issue Nov 21, 2020 · 5 comments
Open

How should we interpret SR of the L line in rGFA #19

SHuang-Broad opened this issue Nov 21, 2020 · 5 comments

Comments

@SHuang-Broad
Copy link

Hello Heng,

While I'm playing with the example assemblies on the FTP site, I noticed that (all) the L-lines in the rGFA file also have the SR tag, which isn't a required tag according to the rGFA format (only the S-Lines are required to).

So how should we interpret them?

Thanks!
Steve

@lh3
Copy link
Owner

lh3 commented Nov 21, 2020

SR on L-lines plays the same role as SR on S-lines. Think about a deletion from the reference. If you don't keep SR on L-lines/links, you don't know how to traverse the reference path.

@SHuang-Broad
Copy link
Author

What about L-lines that connects two segments of different ranks?
Would the L-line take the rank from the segment of higher priority (i.e. lower rank value)?

@lh3
Copy link
Owner

lh3 commented Nov 21, 2020

A link shouldn't have a lower rank than the two segments. Think about the minigraph incremental graph construction procedure. SR on L-line is a natural outcome.

@SHuang-Broad
Copy link
Author

Sorry, I meant a link that links together two segments having different ranks,

e.g. a link that links together one segment from GRCh38 and one segment from CHM13 (T2T).

If the rGFA file was built started with GRCh38 (hence rank 0) then with CHM13 (hence rank 1),
then should I expect the link in this example to have rank 0?

@lh3
Copy link
Owner

lh3 commented Dec 17, 2020

Probably not.

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