Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

An error occurred in convertalis when outputting alntmscore、fident and prob #340

Open
TigerWindWood opened this issue Aug 29, 2024 · 0 comments

Comments

@TigerWindWood
Copy link
Contributor

Thank you for creating such an efficient tool.

Expected Behavior

I would like to obtain a foldseek result that includes alntmscore and the protein sequences of the hits

Current Behavior

$foldseek createdb query.pdb queryDB
$foldseek search queryDB targetDB aln_DB tmp
$foldseek convertalis queryDB targetDB aln_DB result.tsv --format-mode 4 --format-output query,target,alntmscore,fident,evalue,bits,prob,lddt,qstart,qend,tstart,tend,qcov,tcov,tseq
An error occurred: "Backtrace cigar is missing in the alignment result. Please recompute the alignment with the -a flag."
I followed the instructions and executed 'mmseqs align queryDB targetDB aln_DB alnNew -a', but the same error still occurred.
When I executed “$foldseek convertalis queryDB targetDB aln_DB result.tsv --format-mode 4 --format-output query,target,fident,evalue,bits,qstart,qend,tstart,tend,qcov,tcov,tseq,” it successfully output the results.
It seems that convertalis encountered an issue when outputting alntmscore, fident, and prob.

Foldssek Output (for bugs)

MMseqs Version: 9.427df8a
Substitution matrix aa:3di.out,nucl:3di.out
Alignment format 4
Format alignment output query,target,alntmscore,fident,evalue,bits,prob,lddt,qstart,qend,tstart,tend,qcov,tcov,tseq
Gap open cost aa:10,nucl:10
Gap extension cost aa:1,nucl:1
Database output false
Preload mode 0
Threads 150
Compressed 0
Verbosity 3
Exact TMscore 0

[=================================================================] 100.00% 1 eta -
Backtrace cigar is missing in the alignment result. Please recompute the alignment with the -a flag.
Command: mmseqs align queryDB targetDB aln_DB alnNew -a

Your Environment

foldseek Version: 9.427df8a

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant