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

position X+Y is not inside node X #4250

Open
Kled111 opened this issue Mar 19, 2024 · 1 comment
Open

position X+Y is not inside node X #4250

Kled111 opened this issue Mar 19, 2024 · 1 comment

Comments

@Kled111
Copy link

Kled111 commented Mar 19, 2024

I have successfully done with this command -- " vg giraffe -Z hprc-v1.1-mc-chm13.gbz -d /hprc-v1.1-mc-chm13.dist -m hprc-v1.1-mc-chm13.min -p -t 64 -f xxx.fastq > map.gam"
but when I try to do variant calling with this command -- "vg augment chr22.vg map.gam -m 4 -q 5 -Q 5 -s -A chr22_aug.gam > chr22_aug.vg", it failed.
This is the Crash report for vg v1.55.0 "Bernolda"
forwardize_breakpoints error: failure, position 45209079-3 is not inside node 45209079
vg: src/augment.cpp:514: std::unordered_map<long long int, std::set<std::tuple<long long int, bool, long unsigned int> > > vg::forwardize_breakpoints(const HandleGraph*, const std::unordered_map<long long int, std::set<std::tuple<long long int, bool, long unsigned int> > >&): Assertion `false' failed.
━━━━━━━━━━━━━━━━━━━━
What should I do?

@adamnovak
Copy link
Member

What is the relationship between hprc-v1.1-mc-chm13.gbz and chr22.vg? Are you sure that map.gam contains reads that reference nodes in chr22.vg and only nodes in chr22.vg?

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