X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Futil%2FComparison.java;h=835a1b4c58c20712fbc58eddb22c08e6239a64d8;hb=17fe82c9d7f9da1b17f3091e4b3992a35e3e411e;hp=e224b710262d78ad5bb1a78cbeaa0b7100f0ad1e;hpb=7245a524318f1e426e0e0c68109a4291c4e77be5;p=jalview.git

diff --git a/src/jalview/util/Comparison.java b/src/jalview/util/Comparison.java
index e224b71..835a1b4 100644
--- a/src/jalview/util/Comparison.java
+++ b/src/jalview/util/Comparison.java
@@ -20,6 +20,9 @@
  */
 package jalview.util;
 
+import java.util.ArrayList;
+import java.util.List;
+
 import jalview.datamodel.SequenceI;
 
 /**
@@ -248,7 +251,7 @@ public class Comparison
   /**
    * Answers true if more than 85% of the sequence residues (ignoring gaps) are
    * A, G, C, T or U, else false. This is just a heuristic guess and may give a
-   * wrong answer (as AGCT are also animo acid codes).
+   * wrong answer (as AGCT are also amino acid codes).
    * 
    * @param seqs
    * @return
@@ -263,6 +266,12 @@ public class Comparison
     int aaCount = 0;
     for (SequenceI seq : seqs)
     {
+      if (seq == null)
+      {
+        continue;
+      }
+      // TODO could possibly make an informed guess just from the first sequence
+      // to save a lengthy calculation
       for (char c : seq.getSequence())
       {
         if ('a' <= c && c <= 'z')
@@ -295,4 +304,27 @@ public class Comparison
     }
 
   }
+
+  /**
+   * Convenience overload of isNucleotide
+   * 
+   * @param seqs
+   * @return
+   */
+  public static boolean isNucleotide(SequenceI[][] seqs)
+  {
+    if (seqs == null)
+    {
+      return false;
+    }
+    List<SequenceI> flattened = new ArrayList<SequenceI>();
+    for (SequenceI[] ss : seqs)
+    {
+      for (SequenceI s : ss) {
+      flattened.add(s);
+      }
+    }
+    final SequenceI[] oneDArray = flattened.toArray(new SequenceI[flattened.size()]);
+    return isNucleotide(oneDArray);
+  }
 }