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

Assertion `reference_for.count(fasta_contig)' failed. #3702

Open
8banzhuan opened this issue Jul 17, 2022 · 1 comment
Open

Assertion `reference_for.count(fasta_contig)' failed. #3702

8banzhuan opened this issue Jul 17, 2022 · 1 comment

Comments

@8banzhuan
Copy link

Hi I'm having a little trouble trying to build a vg graph
I tried to build with the following command I provided GRC38 and its VCF file
./vg construct -r fa/GRCh38_latest_genomic.fna -v vcf/GRCh38_latest_clinvar.vcf > gfa/human.vg
He gives the following error
vg: src/constructor.cpp:2230: void vg::Constructor::construct_graph(const std::vector<FastaReference*>&, const std::vectorvcflib::VariantCallFile*&, const std::vector<FastaReference*>&, const std::function<void(vg::Graph&)>&): Assertion reference_for.count(fasta_contig)' failed. ERROR: Signal 6 occurred. VG has crashed. Visit https://github.com/vgteam/vg/issues/new/choose to report a bug. Stack trace path: /tmp/vg_crash_JCvAIz/stacktrace.txt Please include the stack trace file in your bug report! The content in stacktrace.txt is Crash report for vg v1.38.0 "Canossa"
Stack trace (most recent call last):
#10 Object "/home/yaozhen/vg", at 0x5b635d, in _start
#9 Object "/home/yaozhen/vg", at 0x1cb357f, in __libc_start_main
#8 Object "/home/yaozhen/vg", at 0x58a62e, in main
#7 Object "/home/yaozhen/vg", at 0xbf999b, in vg::subcommand::Subcommand::operator()(int, char**) const
#6 Object "/home/yaozhen/vg", at 0xbb6076, in main_construct(int, char**)
#5 Object "/home/yaozhen/vg", at 0xd1d161, in vg::Constructor::construct_graph(std::vector<std::__cxx11::basic_string<char, std::char_traits, std::allocator >, std::allocator<std::__cxx11::basic_string<char, std::char_traits, std::allocator > > > const&, std::vector<std::__cxx11::basic_string<char, std::char_traits, std::allocator >, std::allocator<std::__cxx11::basic_string<char, std::char_traits, std::allocator > > > const&, std::vector<std::__cxx11::basic_string<char, std::char_traits, std::allocator >, std::allocator<std::__cxx11::basic_string<char, std::char_traits, std::allocator > > > const&, std::function<void (vg::Graph&)> const&)
#4 Object "/home/yaozhen/vg", at 0xd1c541, in vg::Constructor::construct_graph(std::vector<FastaReference*, std::allocator<FastaReference*> > const&, std::vector<vcflib::VariantCallFile*, std::allocatorvcflib::VariantCallFile* > const&, std::vector<FastaReference*, std::allocator<FastaReference*> > const&, std::function<void (vg::Graph&)> const&)
#3 Object "/home/yaozhen/vg", at 0x1cc3855, in __assert_fail
#2 Object "/home/yaozhen/vg", at 0x589a23, in __assert_fail_base.cold
#1 Object "/home/yaozhen/vg", at 0x589b53, in abort
#0 Object "/home/yaozhen/vg", at 0x136b1cb, in raise
`
Excuse me, thank you

@jeizenga
Copy link
Contributor

Can you verify that all of the chromosomes have the same name in both the FASTA and the VCF? For instance, you might have the sequence called chr12 in the FASTA and 12 in the VCF.

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