MMseqs2
MMseqs2 copied to clipboard
Including `taln` column in output segfaults
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)
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?
That's weird. Could you share the two fasta files so I can try to reproduce the issue