From 3fabec8a0391ddc43b1f0549a4ffd33e068cbb57 Mon Sep 17 00:00:00 2001
From: Amandine PERRIN <amandine.perrin@pasteur.fr>
Date: Tue, 5 Sep 2023 09:15:17 +0200
Subject: [PATCH] Change --anysymbol to --amino for prot alignment
---
PanACoTA/align_module/alignment.py | 2 +-
test/test_unit/test_align/test_alignment.py | 10 +++++-----
2 files changed, 6 insertions(+), 6 deletions(-)
diff --git a/PanACoTA/align_module/alignment.py b/PanACoTA/align_module/alignment.py
index 1618f4a8..0344c44d 100755
--- a/PanACoTA/align_module/alignment.py
+++ b/PanACoTA/align_module/alignment.py
@@ -557,7 +557,7 @@ def mafft_align(num_fam, prt_file, mafft_file, nbfprt, logger):
False otherwise
"""
logger.log(utils.detail_lvl(), f"Aligning family {num_fam}")
- cmd = f"mafft --auto --anysymbol {prt_file}"
+ cmd = f"mafft --auto --amino {prt_file}"
error = f"Problem while trying to align fam {num_fam}"
stdout = open(mafft_file, "w")
stderr = open(mafft_file + ".log", "w")
diff --git a/test/test_unit/test_align/test_alignment.py b/test/test_unit/test_align/test_alignment.py
index 8260be8a..1102e12a 100755
--- a/test/test_unit/test_align/test_alignment.py
+++ b/test/test_unit/test_align/test_alignment.py
@@ -470,7 +470,7 @@ def test_family_align_mafftempty_btrempty(caplog):
"in test/data/align/exp_files/exp_aldir-pers/current.8.prt (3) and proteins aligned "
"in existing test/data/align/generated_by_unit-tests/test_fam_align.8.aln") in caplog.text
assert "Aligning family 8" in caplog.text
- assert "Mafft command: mafft --auto --anysymbol test/data/align/exp_files/exp_aldir-pers/current.8.prt" in caplog.text
+ assert "Mafft command: mafft --auto --amino test/data/align/exp_files/exp_aldir-pers/current.8.prt" in caplog.text
assert "Back-translating family 8" in caplog.text
# Check content of mafft and btr files
exp_mafft = os.path.join(EXPPATH, "exp_aldir-pers", "mafft-align.8.aln")
@@ -577,7 +577,7 @@ def test_family_align_nomafft_btrempty_errormafft(caplog):
ngenomes, logger)
assert "Checking extractions for family 8" in caplog.text
assert "Aligning family 8" in caplog.text
- assert ("command '>mafft --auto --anysymbol test/data/align/exp_files/exp_aldir-pers/current.8.prt' "
+ assert ("command '>mafft --auto --amino test/data/align/exp_files/exp_aldir-pers/current.8.prt' "
"is not possible") in caplog.text
# Check content of mafft and btr files
assert not os.path.isfile(mafft_file)
@@ -808,7 +808,7 @@ def test_handle_family_true():
q.put(None)
assert "Checking extractions for family 8" in q.get().message
assert "Aligning family 8" in q.get().message
- assert ("mafft --auto --anysymbol test/data/align/generated_by_unit-tests/aldir/"
+ assert ("mafft --auto --amino test/data/align/generated_by_unit-tests/aldir/"
"TESThandlefam-current.8.prt") in q.get().message
assert "Back-translating family 8" in q.get().message
# Writes it twice: 1 for nucleic and one for aa
@@ -849,7 +849,7 @@ def test_handle_family_true_nomiss():
q.put(None)
assert "Checking extractions for family 8" in q.get().message
assert "Aligning family 8" in q.get().message
- assert ("mafft --auto --anysymbol test/data/align/generated_by_unit-tests/aldir/"
+ assert ("mafft --auto --amino test/data/align/generated_by_unit-tests/aldir/"
"TESThandlefam-current.8.prt") in q.get().message
assert "Back-translating family 8" in q.get().message
assert not q.get()
@@ -894,7 +894,7 @@ def test_handle_family_emptyaln_true():
"and proteins aligned in existing "
"test/data/align/generated_by_unit-tests/aldir/TESThandlefam-mafft-align.8.aln") in q.get().message
assert "Aligning family 8" in q.get().message
- assert ("mafft --auto --anysymbol test/data/align/generated_by_unit-tests/aldir/"
+ assert ("mafft --auto --amino test/data/align/generated_by_unit-tests/aldir/"
"TESThandlefam-current.8.prt") in q.get().message
assert "Back-translating family 8" in q.get().message
assert "Adding missing genomes for family 8 in protein alignment" in q.get().message
--
GitLab