MMseqs2 icon indicating copy to clipboard operation
MMseqs2 copied to clipboard

Including `taln` column in output segfaults

Open rickbeeloo opened this issue 2 months ago • 2 comments

When running:

mmseqs easy-search  q.fasta   t.fasta   mmseqs_tmp  mmseqs_results   --format-output "query,target,pident,alnlen,mismatch,gapopen,qstart,qend,tstart,tend,evalue,bits,qcov,tcov,qlen,qaln"  --threads 60 --search-type 3

it works fine, but adding taln to the output column breaks it:

mmseqs easy-search  q.fasta   t.fasta   mmseqs_tmp  mmseqs_results   --format-output "query,target,pident,alnlen,mismatch,gapopen,qstart,qend,tstart,tend,evalue,bits,qcov,tcov,qlen,qaln,taln"  --threads 60 --search-type 3

With error:

Segmentation fault (core dumped)                                  ] 6.32% 7 eta 0s       
Error: Convert Alignments died

Target and query Fastas are valid (no empty sequences, no invalid chars, no empty headers, valid header starts (>), etc)

rickbeeloo avatar Oct 16 '25 12:10 rickbeeloo

Default max-seq length is 65535, the longest sequence in both q/t.fasta is 51777, but increasing the max len to --max-seq-len 66000 does work for some reason?

rickbeeloo avatar Oct 16 '25 12:10 rickbeeloo

That's weird. Could you share the two fasta files so I can try to reproduce the issue

milot-mirdita avatar Oct 16 '25 12:10 milot-mirdita