X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;ds=sidebyside;f=src%2Fjalview%2Fanalysis%2FDna.java;h=5dcf212299933186ca544d7d23777cbbfeff9d9f;hb=57738a1f3c19b1c3a00bd3ac5108f8cd0af32f99;hp=4b817858fe5106cf54789fb011687ebbc9dc7935;hpb=26b115b0a77d521da92a06572d9b7819c2d0d49a;p=jalview.git

diff --git a/src/jalview/analysis/Dna.java b/src/jalview/analysis/Dna.java
index 4b81785..5dcf212 100644
--- a/src/jalview/analysis/Dna.java
+++ b/src/jalview/analysis/Dna.java
@@ -163,7 +163,8 @@ public class Dna
    * @param ac2
    * @return
    */
-  public static final int compareCodonPos(AlignedCodon ac1, AlignedCodon ac2)
+  public static final int compareCodonPos(AlignedCodon ac1,
+          AlignedCodon ac2)
   {
     return comparator.compare(ac1, ac2);
     // return jalview_2_8_2compare(ac1, ac2);
@@ -194,10 +195,11 @@ public class Dna
   }
 
   /**
+   * Translates cDNA using the specified code table
    * 
    * @return
    */
-  public AlignmentI translateCdna()
+  public AlignmentI translateCdna(GeneticCodeI codeTable)
   {
     AlignedCodonFrame acf = new AlignedCodonFrame();
 
@@ -209,7 +211,7 @@ public class Dna
     for (s = 0; s < sSize; s++)
     {
       SequenceI newseq = translateCodingRegion(selection.get(s),
-              seqstring[s], acf, pepseqs);
+              seqstring[s], acf, pepseqs, codeTable);
 
       if (newseq != null)
       {
@@ -269,7 +271,7 @@ public class Dna
             mappedrefs.add(ref); // add translated protein maps
           }
         }
-        dnarefs = mappedrefs;//.toArray(new DBRefEntry[mappedrefs.size()]);
+        dnarefs = mappedrefs;// .toArray(new DBRefEntry[mappedrefs.size()]);
         for (int d = 0, nd = dnarefs.size(); d < nd; d++)
         {
           Mapping mp = dnarefs.get(d).getMap();
@@ -430,11 +432,12 @@ public class Dna
    * @param acf
    *          Definition of global ORF alignment reference frame
    * @param proteinSeqs
+   * @param codeTable
    * @return sequence ready to be added to alignment.
    */
   protected SequenceI translateCodingRegion(SequenceI selection,
           String seqstring, AlignedCodonFrame acf,
-          List<SequenceI> proteinSeqs)
+          List<SequenceI> proteinSeqs, GeneticCodeI codeTable)
   {
     List<int[]> skip = new ArrayList<>();
     int[] skipint = null;
@@ -469,7 +472,7 @@ public class Dna
          */
         AlignedCodon alignedCodon = new AlignedCodon(cdp[0], cdp[1],
                 cdp[2]);
-        String aa = ResidueProperties.codonTranslate(new String(codon));
+        String aa = codeTable.translate(new String(codon));
         rf = 0;
         final String gapString = String.valueOf(gapChar);
         if (aa == null)
@@ -797,29 +800,29 @@ public class Dna
   private static void transferCodedFeatures(SequenceI dna, SequenceI pep,
           MapList map)
   {
-	//  BH 2019.01.25 nop?
-//    List<DBRefEntry> dnarefs = DBRefUtils.selectRefs(dna.getDBRefs(),
-//            DBRefSource.DNACODINGDBS);
-//    if (dnarefs != null)
-//    {
-//      // intersect with pep
-//      for (int d = 0, nd = dnarefs.size(); d < nd; d++)
-//      {
-//        Mapping mp = dnarefs.get(d).getMap();
-//        if (mp != null)
-//        {
-//        }
-//      }
-//    }
+    // BH 2019.01.25 nop?
+    // List<DBRefEntry> dnarefs = DBRefUtils.selectRefs(dna.getDBRefs(),
+    // DBRefSource.DNACODINGDBS);
+    // if (dnarefs != null)
+    // {
+    // // intersect with pep
+    // for (int d = 0, nd = dnarefs.size(); d < nd; d++)
+    // {
+    // Mapping mp = dnarefs.get(d).getMap();
+    // if (mp != null)
+    // {
+    // }
+    // }
+    // }
     for (SequenceFeature sf : dna.getFeatures().getAllFeatures())
     {
-        if (FeatureProperties.isCodingFeature(null, sf.getType()))
+      if (FeatureProperties.isCodingFeature(null, sf.getType()))
+      {
+        // if (map.intersectsFrom(sf[f].begin, sf[f].end))
         {
-          // if (map.intersectsFrom(sf[f].begin, sf[f].end))
-          {
 
-          }
         }
+      }
     }
   }