realign commands doesn't seem to work

[Eric-Pelletier]

I want to apply the realign method on a cluster, but it seems to not function :

diamond makedb --in test.faa -d test

diamond deepclust -d test -o tset.dc -p 8
...
Clustering round 4 complete. #Input sequences: 18032 #Clusters: 16458 Time: 4s
Generating output...  [0.007s]

diamond realign -d test --clusters tset.dc -o tset.al -f 6 qseqid qlen qstart qend sseqid slen sstart send length approx_pident ppos score evalue bitscore
diamond v2.1.1.155 (C) Max Planck Society for the Advancement of Science
Documentation, support and updates available at http://www.diamondsearch.org
Please cite: http://dx.doi.org/10.1038/s41592-021-01101-x Nature Methods (2021)

Scoring parameters: (Matrix=BLOSUM62 Lambda=0.267 K=0.041 Penalties=11/1)
Opening the database...  [0.022s]
#Database sequences: 20834, #Letters: 4117251
Found 16458 centroids, 20834 mappings in input file.
Block size: 8000000000 byte.
Loading centroid block...  [0.005s]
Processing centroid block 1, member block 1... terminate called after throwing an instance of 'std::runtime_error'
what():  Traceback error. MRLQHAAVA...

Any idea ?

-- Eric Pelletier

[bbuchfink]

I can reproduce the problem and am working on a fix. In the meantime, it should work without the ppos field.

[Eric-Pelletier]

I can reproduce the problem and am working on a fix. In the meantime, it should work without the ppos field.

It works, indeed. :)

[Eric-Pelletier]

I can reproduce the problem and am working on a fix. In the meantime, it should work without the ppos field.

It works, indeed. :)

Except that the qlen field always states 0 (while slen gives the correct value)

[bbuchfink]

This should be fixed in the latest release.