PASSIONLab / BELLA

BELLA: a Computationally-Efficient and Highly-Accurate Long-Read to Long-Read Aligner and Overlapper
Other
48 stars 8 forks source link

Freeing invalid pointer when using non-zero window size for minimizer selection #30

Closed gabe-raulet closed 3 years ago

gabe-raulet commented 3 years ago

Hi Giulia,

I've run into a problem when I try to run Bella using minimizers. Specifically, when I run the command

"bella -f input.txt -o output --paf -k 31 -w 17 -e 0.005 -l 2 -u 100",

I get the following output log:

INFO: src/../include/kmercount.hpp(108) InputFile = metagenome-anonymous.fastq INFO: src/../include/kmercount.hpp(109) InputSize = 67.693069 MB INFO: src/main.cpp(200) OutputFile = output.out INFO: src/main.cpp(203) kmerSize = 31 INFO: src/main.cpp(206) GPUs = DISABLED INFO: src/main.cpp(209) UserDefinedMemory = 8000.000000 MB INFO: src/main.cpp(212) OutputPAF = 1 INFO: src/main.cpp(215) BinSize = 500 INFO: src/main.cpp(218) DeltaChernoff = 0.100000 INFO: src/main.cpp(221) RunPairwiseAlignment = 1 INFO: src/main.cpp(231) HOPC = DISABLED INFO: src/main.cpp(235) xDrop = 7 INFO: src/main.cpp(238) KmerSplitCount = 1 INFO: src/main.cpp(243) useMinimizer = ENABLED INFO: src/main.cpp(245) minimizerWindow = 17 INFO: src/main.cpp(276) numThreads = 64 INFO: src/../include/kmercount.hpp(711) ReadingFASTQ = metagenome-anonymous.fastq Error in `bella': free(): invalid pointer: 0x00002aaacc076480 Aborted

When I run without the "-w 17" option, it runs perfectly.

Thanks, Gabe