bustools count segmentation fault

[rmandla]

Hi,

I am running bustools 0.39.3 on Manjaro 5.4 and have recently got a segmentation fault when running the command:

sudo ~/Downloads/bustools/bustools count -o eqcount/tcc -g ../../transcripts_to_genes.txt -e matrix.ec -t transcripts.txt output.correct.sort.bus

I get the following error:

[1] 2595173 segmentation fault sudo ~/Downloads/bustools/bustools count -o eqcount/tcc -g -e matrix.ec -t

I have previously used bustools count successfully. The only difference this time is that my pseudoalignment includes 3 pairs of fq files rather than a single pair which I tested previously. Also, genecounts seem to work fine, and it is just the TCC counting which I am having trouble with.

Thanks!

And here are the logs for all my upstream kallisto/bustools commands and outputs for reference:

$ sudo ~/Downloads/kallisto/kallisto bus -i /mnt/labshare/ravi/genomes/encode/Mus_musculus.GRCm38.cdna.all.idx -o s1-fastp_bus_output/ -x 10xv2 -t 4 /mnt/labshare/DATASETS/chen-scrna/S1_first_sequencing_R1.fastq.gz /mnt/labshare/DATASETS/chen-scrna/S1_first_sequencing_R2.fastq.gz /mnt/labshare/DATASETS/chen-scrna/S5_first_sequencing_R1.fastq.gz /mnt/labshare/DATASETS/chen-scrna/S5_first_sequencing_R2.fastq.gz /mnt/labshare/DATASETS/chen-scrna/S5_deeper_sequencing_R1.fastq.gz /mnt/labshare/DATASETS/chen-scrna/S5_deeper_sequencing_R2.fastq.gz

[index] k-mer length: 31
[index] number of targets: 118,489
[index] number of k-mers: 100,614,952
[index] number of equivalence classes: 433,624
[quant] will process sample 1: /mnt/labshare/DATASETS/chen-scrna/S1_first_sequencing_R1.fastq.gz
/mnt/labshare/DATASETS/chen-scrna/S1_first_sequencing_R2.fastq.gz
[quant] will process sample 2: /mnt/labshare/DATASETS/chen-scrna/S5_first_sequencing_R1.fastq.gz
/mnt/labshare/DATASETS/chen-scrna/S5_first_sequencing_R2.fastq.gz
[quant] will process sample 3: /mnt/labshare/DATASETS/chen-scrna/S5_deeper_sequencing_R1.fastq.gz
/mnt/labshare/DATASETS/chen-scrna/S5_deeper_sequencing_R2.fastq.gz
[quant] finding pseudoalignments for the reads ... done
[quant] processed 400,009,570 reads, 308,864,287 reads pseudoaligned

$ sudo ~/Downloads/bustools/bustools correct -w ../10xv2_whitelist.txt -o output.correct.bus output.bus
Found 737280 barcodes in the whitelist
Number of hamming dist 1 barcodes = 20550336
Processed 308864287 bus records
In whitelist = 297255718
Corrected = 3547753

$ sudo ~/Downloads/bustools/bustools sort -t 4 -o output.correct.sort.bus output.correct.bus
Read in 300803471 BUS records