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

Recognize sample#contig paths in GFA to GBWTGraph #4190

Merged
merged 3 commits into from
Dec 10, 2023
Merged

Conversation

jltsiren
Copy link
Contributor

@jltsiren jltsiren commented Dec 10, 2023

Changelog Entry

To be copied to the draft changelog by merger:

  • GBWTGraph algorithm for parsing GFA now handles P-line names of the form sample#contig correctly.

Description

Function get_best_gbwtgraph_gfa_parsing_parameters() was missing name pattern sample#contig, which is sometimes used for reference paths stored as P-lines. This PR adds the pattern. When there are no W-lines in the GFA, paths with a name matching that pattern will be parsed as haplotype paths with no haplotype/phase number. They can be promoted to reference paths using tag RS in the header.

This resolves #4110.

@jltsiren jltsiren merged commit faf6499 into master Dec 10, 2023
1 check passed
@jltsiren jltsiren deleted the sample-contig-paths branch December 10, 2023 08:55
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

Successfully merging this pull request may close these issues.

vg gbwt -G does not parse two-part GRCh38#chr1 paths
2 participants