proteinaln2nucl did not use the threads parameter

soedinglab · Apr 16, 2024 · 0e2ef17 · 0e2ef17
1 parent ff8fc36
commit 0e2ef17
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 1 deletion.
diff --git a/data/guidedNuclAssemble.sh b/data/guidedNuclAssemble.sh
@@ -103,7 +103,7 @@ while [ $STEP -lt $NUM_IT ]; do
 
  # 3. Ungapped alignment protein 2 nucl
  if notExists "${TMP_PATH_GUIDED_ASSEMBLY}/aln_nucl_$STEP.done"; then
- "$MMSEQS" proteinaln2nucl "$INPUT_NUCL" "$INPUT_NUCL" "$INPUT_AA" "$INPUT_AA" "${TMP_PATH_GUIDED_ASSEMBLY}/aln_$STEP" "${TMP_PATH_GUIDED_ASSEMBLY}/aln_nucl_$STEP" \
+ "$MMSEQS" proteinaln2nucl "$INPUT_NUCL" "$INPUT_NUCL" "$INPUT_AA" "$INPUT_AA" "${TMP_PATH_GUIDED_ASSEMBLY}/aln_$STEP" "${TMP_PATH_GUIDED_ASSEMBLY}/aln_nucl_$STEP" ${PROTEIN_ALN_2_NUCL_PAR} \
  || fail "Ungapped alignment 2 nucl step died"
  deleteIncremental "$PREV_ALN_NUCL"
  touch "${TMP_PATH_GUIDED_ASSEMBLY}/aln_nucl_${STEP}.done"

diff --git a/src/workflow/GuidedNuclassembler.cpp b/src/workflow/GuidedNuclassembler.cpp
@@ -157,6 +157,8 @@ int guidedNuclAssemble(int argc, const char **argv, const Command &command) {
  bool addBacktrace = par.addBacktrace;
  par.addBacktrace = true;
  cmd.addVariable("UNGAPPED_ALN_PAR", par.createParameterString(par.rescorediagonal).c_str());
+ cmd.addVariable("PROTEIN_ALN_2_NUCL_PAR", par.createParameterString(par.proteinaln2nucl).c_str());
+
 
  // # 3. Assembly: Extend by left and right extension
  par.seqIdThr = par.multiSeqIdThr.nucleotides;