vg icon indicating copy to clipboard operation
vg copied to clipboard

vg construct throwing error "what(): Found unexpected delimiter in fastq/fasta input"

Open champ1994 opened this issue 3 years ago • 0 comments

1. What were you trying to do? I was trying to construct graph using construct module from /test/1mb1kgp fasta and vcf

2. What did you want to happen? I ran the command

./vg construct -r z.fa -v z.vcf.gz -S -p > z.vg 

3. What actually happened? It threw following error

terminate called after throwing an instance of 'std::runtime_error'
  what():  Found unexpected delimiter  in fastq/fasta input
ERROR: Signal 6 occurred. VG has crashed. 

4. If you got a line like Stack trace path: /somewhere/on/your/computer/stacktrace.txt, please copy-paste the contents of that file here:

Crash report for vg v1.30.0 "Carentino"
Stack trace (most recent call last):
#13   Object "", at 0xffffffffffffffff, in 
#12   Object "/data1/VC_exp/software/vg-v1.30.0/bin/vg", at 0x55cd54e46719, in _start
#11   Object "/lib/x86_64-linux-gnu/libc-2.27.so", at 0x7f0ac5eadbf6, in __libc_start_main
#10   Object "/data1/VC_exp/software/vg-v1.30.0/bin/vg", at 0x55cd54e22bbb, in main
      Source "src/main.cpp", line 81, in main [0x55cd54e22bbb]
#9    Object "/data1/VC_exp/software/vg-v1.30.0/bin/vg", at 0x55cd554034f7, in vg::subcommand::Subcommand::operator()(int, char**) const
    | Source "src/subcommand/subcommand.cpp", line 75, in operator()
      Source "/usr/include/c++/7/bits/std_function.h", line 706, in operator() [0x55cd554034f7]
        703:     {
        704:       if (_M_empty())
        705: 	__throw_bad_function_call();
      > 706:       return _M_invoker(_M_functor, std::forward<_ArgTypes>(__args)...);
        707:     }
        708: 
        709: #if __cpp_rtti
#8    Object "/data1/VC_exp/software/vg-v1.30.0/bin/vg", at 0x55cd5540c74a, in main_view(int, char**)
      Source "src/subcommand/view_main.cpp", line 765, in main_view [0x55cd5540c74a]
#7    Object "/data1/VC_exp/software/vg-v1.30.0/bin/vg", at 0x55cd556ac662, in vg::fastq_unpaired_for_each(std::__cxx11::basic_string<char, st
d::char_traits<char>, std::allocator<char> > const&, std::function<void (vg::Alignment&)>)
      Source "src/alignment.cpp", line 272, in fastq_unpaired_for_each [0x55cd556ac662]
#6    Object "/data1/VC_exp/software/vg-v1.30.0/bin/vg", at 0x55cd556a787e, in vg::get_next_alignment_from_fastq(gzFile_s*, char*, unsigned lo
ng, vg::Alignment&)
      Source "src/alignment.cpp", line 139, in get_next_alignment_from_fastq [0x55cd556a787e]
#5    Object "/usr/lib/x86_64-linux-gnu/libstdc++.so.6.0.25", at 0x7f0ac6756d53, in __cxa_throw
#4    Object "/usr/lib/x86_64-linux-gnu/libstdc++.so.6.0.25", at 0x7f0ac6756b20, in std::terminate()
#3    Object "/usr/lib/x86_64-linux-gnu/libstdc++.so.6.0.25", at 0x7f0ac6756ae5, in 
#2    Object "/usr/lib/x86_64-linux-gnu/libstdc++.so.6.0.25", at 0x7f0ac6750956, in 
#1    Object "/lib/x86_64-linux-gnu/libc-2.27.so", at 0x7f0ac5ecc920, in abort
#0    Object "/lib/x86_64-linux-gnu/libc-2.27.so", at 0x7f0ac5ecafb7, in raise

5. What data and command can the vg dev team use to make the problem happen? Command : ./vg construct -r z.fa -v z.vcf.gz -S -p > z.vg Data : /test/1mb1kgp

6. What does running vg version say?

vg version v1.30.0 "Carentino"
Compiled with g++ (Ubuntu 7.5.0-3ubuntu1~18.04) 7.5.0 on Linux
Linked against libstd++ 20191114

champ1994 avatar Jan 12 '22 05:01 champ1994