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

Genome Alignment: Ref doesn't show gaps in translocations #74

Open
josiahseaman opened this issue Dec 18, 2018 · 0 comments
Open

Genome Alignment: Ref doesn't show gaps in translocations #74

josiahseaman opened this issue Dec 18, 2018 · 0 comments
Assignees
Labels
Milestone

Comments

@josiahseaman
Copy link
Owner

josiahseaman commented Dec 18, 2018

Followup #33.
When a chain representing a translocation is processed, the series of AlignedSpans are inserted into the existing composite alignment. These new alignedSpans do not appear to have query_tail_size > 0.
ChainParser.py ~ line 219

image

Observations

  • Nematode alignment show no gaps anywhere
  • In test_translocations with the full chr21, there are places with no ref gaps that still contain query gaps (image above)
  • In test_whole_chromosome_alignment the master_chain shows gaps in both, other chains show gaps in neither.
  • When that same chr7 N's chain is the largest chain in the file, there are query gaps, but no ref gaps.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

1 participant