zjshi / gt-pro

MIT License
23 stars 7 forks source link

GT_pro optimize terminate called after throwing an instance of 'std::length_error' what(): vector::_M_default_append #60

Open HaimAshk opened 9 months ago

HaimAshk commented 9 months ago

Hi,

I try to optimize a custom build DB (2 species) and get the following error. I would appreciate any advice!

GT_pro optimize --db db_build/Pseudomonas_Sphingomonas_db --in S2MetagenomicR1R2.fq.gz --tmp /tmp
[OK] start initial optimization
[OK] database found
[OK] optimize from db_build/Pseudomonas_Sphingomonas_db_optimized_db_kmer_index.bin and db_build/Pseudomonas_Sphingomonas_db_optimized_db_snps.bin

Error: the following returned non-zero status: 'gt_pro -d db_build/Pseudomonas_Sphingomonas_db </dev/null':

b"gt_pro\tdb_build/Pseudomonas_Sphingomonas_db\t256\tno_overwrite\n1699636911863:  [Info] Starting to load DB: db_build/Pseudomonas_Sphingomonas_db\n1699636911863:  [Info] MMAPPING db_build/Pseudomonas_Sphingomonas_db_optimized_db_snps.bin\n1699636911868:  [Info] MMAPPING db_build/Pseudomonas_Sphingomonas_db_optimized_db_kmer_index.bin\n1699636911924:  [Info] Using -l 32 -m 36 as optimal for system RAM\n1699636911924:  [Info] MMAPPING db_build/Pseudomonas_Sphingomonas_db_optimized_db_mmer_bloom_36.bin\n1699636917599:  [Info] MMAPPING db_build/Pseudomonas_Sphingomonas_db_optimized_db_lmer_index_32.bin\n1699636940430:  [Info] Done with init for optimized DB with 5003356 kmers.  That took 28 seconds.\n1699636940430:  [Info] Will input reads from stdin and output snps to stdout.\n1699636940430:  [Info] Output will appear only after stdin reaches EOF.\nterminate called after throwing an instance of 'std::length_error'\n  what():  vector::_M_default_append\nAborted\n"

The reads are 150bp interleaved-paired-end fastq file

Thanks! Haim