Removed * aa from dna translation at the end of sequence

TeselaGen · Nov 27, 2023 · 436b190 · 436b190
1 parent 860dc2f
commit 436b190
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 1 deletion.
diff --git a/packages/sequence-utils/src/getAminoAcidStringFromSequenceString.js b/packages/sequence-utils/src/getAminoAcidStringFromSequenceString.js
@@ -7,10 +7,14 @@ export default function getAminoAcidStringFromSequenceString(sequenceString) {
   );
   const aaArray = [];
   let aaString = "";
-  aminoAcidsPerBase.forEach(aa => {
+  aminoAcidsPerBase.forEach((aa, index) => {
     if (!aa.fullCodon) {
       return;
     }
+    // Check if the current amino acid is the last in the sequence and is a stop codon
+    if (index === aminoAcidsPerBase.length - 1 && aa.aminoAcid.value === '*') {
+      return;
+    }
     aaArray[aa.aminoAcidIndex] = aa.aminoAcid.value;
   });
   aaString = aaArray.join("");

diff --git a/packages/sequence-utils/src/getAminoAcidStringFromSequenceString.test.js b/packages/sequence-utils/src/getAminoAcidStringFromSequenceString.test.js
@@ -14,5 +14,6 @@ describe("getAminoAcidStringFromSequenceString", () => {
     assert.equal("MM", getAminoAcidStringFromSequenceString("atgatg"));
     assert.equal("M--", getAminoAcidStringFromSequenceString("atg------"));
     assert.equal("", getAminoAcidStringFromSequenceString("at"));
+    assert.equal("MTNYNQKNEN", getAminoAcidStringFromSequenceString("atgactaattataatcaaaaaaatgaaaattaa"));
   });
 });