-
Notifications
You must be signed in to change notification settings - Fork 0
Open
Labels
Description
Zagros crashes with an assertion failure with the following command line:
zagros -c human/hg19/chromFa min_reproducible.bed
where the contents of min_reproducible.bed is the following single line:
chrM 0 16600 X 1026.48 + 1.56762e-15
Resultant error message is
zagros: ../rbp_common/IO.cpp:575: void extract_regions_chrom_fasta(const std::string&, const
std::string&, const std::vector<GenomicRegion, std::allocator<GenomicRegion> >&,
std::vector<std::basic_string<char, std::char_traits<char>, std::allocator<char> >,
std::allocator<std::basic_string<char, std::char_traits<char>, std::allocator<char> > > >&,
std::vector<std::basic_string<char, std::char_traits<char>, std::allocator<char> >,
std::allocator<std::basic_string<char, std::char_traits<char>, std::allocator<char> > > >&): Assertion `i-
>get_width() == sequences.back().length()' failed.
Aborted
Possibly due to exceeding the length of the chromosome? ChrM has length 16571, so this region is more than the whole chromosome...