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

path editing yields dangling nodes when they are soft clips off the end of an existing node #115

Closed
ekg opened this issue Oct 7, 2015 · 1 comment

Comments

@ekg
Copy link
Member

ekg commented Oct 7, 2015

To reproduce, from test/

vg construct -v tiny/tiny.vcf.gz -r tiny/tiny.fa >t.vg
vg msga -g t.vg -s TTCTATAATATG | vg stats -s -

This shows that there are now two disjoint subgraphs (whoops):

1       1
2       55

The original graph:

screenshot from 2015-10-07 18 29 17

Post MSGA:

screenshot from 2015-10-07 18 29 36

The G node corresponds to the last base of the read. It's aligning right at the end of the node and isn't being dropped in properly.

... resolving

@ekg
Copy link
Member Author

ekg commented Oct 7, 2015

It also works in the other direction, such as a soft clip at the beginning of the alignment:

vg msga -g t.vg -s CAAATTTTCTGGAGTTCTAT

screenshot from 2015-10-07 18 35 02

So this seems to be a balanced problem.

@ekg ekg closed this as completed in ec4f2f7 Oct 8, 2015
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

1 participant