patch for JAL-556 - more robust mapping of sequence to structure sequence with correc...
authorjprocter <Jim Procter>
Fri, 11 Jun 2010 09:54:26 +0000 (09:54 +0000)
committerjprocter <Jim Procter>
Fri, 11 Jun 2010 09:54:26 +0000 (09:54 +0000)
src/MCview/AppletPDBCanvas.java
src/jalview/structure/StructureSelectionManager.java

index 6be8dd5..cdad671 100755 (executable)
@@ -176,8 +176,9 @@ public class AppletPDBCanvas extends Panel implements MouseListener,
       // Now lets compare the sequences to get
       // the start and end points.
       // Align the sequence to the pdb
+      // TODO: DNa/Pep switch
       AlignSeq as = new AlignSeq(sequence, ((PDBChain) pdb.chains
-              .elementAt(i)).sequence, "pep");
+              .elementAt(i)).sequence, ((PDBChain) pdb.chains.elementAt(i)).isNa ? AlignSeq.DNA : AlignSeq.PEP);
       as.calcScoreMatrix();
       as.traceAlignment();
       PrintStream ps = new PrintStream(System.out)
index 478706a..483f1a4 100644 (file)
@@ -133,8 +133,10 @@ public class StructureSelectionManager
 
       for (int i = 0; i < pdb.chains.size(); i++)
       {
+        
+        // TODO: correctly determine sequence type for mixed na/peptide structures
         AlignSeq as = new AlignSeq(sequence[s], ((PDBChain) pdb.chains
-                .elementAt(i)).sequence, AlignSeq.PEP);
+                .elementAt(i)).sequence, ((PDBChain)pdb.chains.elementAt(i)).isNa ? AlignSeq.DNA : AlignSeq.PEP);
         as.calcScoreMatrix();
         as.traceAlignment();
         PDBChain chain = ((PDBChain) pdb.chains.elementAt(i));