JAL-2629 exported HMMs are now named correctly
[jalview.git] / src / jalview / hmmer / HMMBuildThread.java
index 48796e4..853b5f1 100644 (file)
@@ -2,6 +2,7 @@ package jalview.hmmer;
 
 import jalview.bin.Cache;
 import jalview.datamodel.AlignmentI;
+import jalview.datamodel.Sequence;
 import jalview.datamodel.SequenceGroup;
 import jalview.datamodel.SequenceI;
 import jalview.gui.AlignFrame;
@@ -23,9 +24,9 @@ public class HMMBuildThread implements Runnable
   HMMERCommands cmds = new HMMERCommands();
   AlignFrame af;
   AlignmentI alignment;
-
   SequenceGroup group;
 
+  boolean forGroup = false;
   boolean canRun = true;
 
   File hmmTemp = null;
@@ -40,6 +41,7 @@ public class HMMBuildThread implements Runnable
     if (af.getViewport().getSelectionGroup() != null)
     {
       group = af.getViewport().getSelectionGroup();
+      forGroup = true;
     }
     alignment = af.getViewport().getAlignment();
 
@@ -77,7 +79,7 @@ public class HMMBuildThread implements Runnable
         List<SequenceI> seqs = alignment
                 .getHMMConsensusSequences(true);
         cmds.setHmmSeqs(seqs);
-        if (group != null)
+        if (forGroup)
         {
           array = group.getSelectionAsNewSequences(alignment);
         }
@@ -89,8 +91,17 @@ public class HMMBuildThread implements Runnable
           }
           array = alignment.getSequencesArray();
         }
-        cmds.uniquifySequences(array);
-        cmds.exportData(array, stoTemp, null, null);
+
+        SequenceI[] newArr = new SequenceI[array.length];
+        int index = 0;
+        for (SequenceI seq : array)
+        {
+          newArr[index] = new Sequence(seq);
+          index++;
+        }
+
+        cmds.uniquifySequences(newArr);
+        cmds.exportData(newArr, stoTemp, null, null);
         jalview.analysis.SeqsetUtils.deuniquify(cmds.hash, array);
 
     } catch (FileNotFoundException e)
@@ -142,8 +153,21 @@ public class HMMBuildThread implements Runnable
     {
       return false;
     }
-    String command = cmds.HMMERFOLDER + cmds.HMMBUILD + cmds.NAME
-            + af.getName() + cmds.SPACE;
+    String command = cmds.HMMERFOLDER + cmds.HMMBUILD + cmds.NAME;
+    if (forGroup)
+    {
+      command += group.getName();
+    }
+    else
+    {
+      String name = af.getName();
+      if (name == null)
+      {
+        name = "Alignment";
+      }
+      command += name;
+    }
+    command += cmds.SPACE;
     if (!alignment.isNucleotide())
     {
       command += cmds.FORCEAMINO; // TODO check for rna