dcjones / quip

Compressing next-generation sequencing data with extreme prejudice.
http://www.cs.washington.edu/homes/dcjones/quip/
BSD 3-Clause "New" or "Revised" License
78 stars 10 forks source link

segfault when creating a qp file when/after "assembling contigs ... done. (0 contigs, 0nt)" #3

Closed ghuls closed 12 years ago

ghuls commented 12 years ago

When creating a qp file, I see segfaults when the number of contigs is 0: $ zcat ./TLE93_read1.clean.fq.gz | head -n2320 | /software/quip-git/src/quip -v > /tmp/TLE93_read1.clean.fq-head.qp writing a block of 58580 compressed bases... counting k-mers ... done. assembling contigs ... done. (0 contigs, 0nt) Segmentation fault When the number of contigs is greater than 0, I see no segfault: $ zcat ./TLE93_read1.clean.fq.gz | head -n2330 | /software/quip-git/src/quip -v > /tmp/TLE93_read1.clean.fq-head.qp writing a block of 58782 compressed bases... counting k-mers ... done. assembling contigs ... done. (1 contigs, 278nt) indexing contigs ... done. aligning reads to contigs ... 9 / 582 [1.55%] reads aligned to contigs (0.00% aborted) done. indexing contigs ... done. id: 3476 / 31029 (11.20%) seq: 14609 / 58782 (24.85%) qual: 32300 / 58782 (54.95%)