Fix issue #393

steineggerlab · Dec 8, 2024 · dd579d9 · dd579d9
1 parent 0dd4b7f
commit dd579d9
Showing 1 changed file with 16 additions and 0 deletions.
diff --git a/src/strucclustutils/structureconvertalis.cpp b/src/strucclustutils/structureconvertalis.cpp
@@ -694,6 +694,20 @@ R"html(<!DOCTYPE html>
                                     Sequence::extractProfileConsensus(targetSeqData, targetEntryLen, *subMat, targetProfData);
                                 }
                             }
+                            for(size_t i = 0; i < outcodes.size(); i++) {
+                                bool computed = false;
+                                switch (outcodes[i]) {
+                                    case LocalParameters::OUTFMT_U:
+                                    case LocalParameters::OUTFMT_T:
+                                        tmres = tmaligner->computeTMscore(targetCaData, &targetCaData[res.dbLen], &targetCaData[res.dbLen+res.dbLen], res.dbLen,
+                                                                          res.qStartPos, res.dbStartPos, tmpBt, res.dbLen);
+                                        computed = true;
+                                        break;
+                                }
+                                if(computed){
+                                    break;
+                                }
+                            }
                             for(size_t i = 0; i < outcodes.size(); i++) {
                                 switch (outcodes[i]) {
                                     case Parameters::OUTFMT_QUERY:
@@ -1022,6 +1036,8 @@ R"html(<!DOCTYPE html>
                         result.append(" ");
                         result.append(targetId);
                         result.append("\n");
+                        tmres = tmaligner->computeTMscore(targetCaData, &targetCaData[res.dbLen], &targetCaData[res.dbLen+res.dbLen], res.dbLen,
+                                                          res.qStartPos, res.dbStartPos, tmpBt, res.dbLen);
                         for(unsigned int tpos = 0; tpos < res.dbLen; tpos++){
                             size_t tId = tDbr->sequenceReader->getId(res.dbKey);
                             char* targetSeqData  = (char*) tDbr->sequenceReader->getData(tId, thread_idx);