blast
[jalview.git] / forester / java / src / org / forester / test / Test.java
index f01c864..5961561 100644 (file)
@@ -54,14 +54,17 @@ import org.forester.io.parsers.nhx.NHXParser;
 import org.forester.io.parsers.phyloxml.PhyloXmlParser;
 import org.forester.io.parsers.tol.TolParser;
 import org.forester.io.writers.PhylogenyWriter;
+import org.forester.msa.BasicMsa;
 import org.forester.msa.Mafft;
 import org.forester.msa.Msa;
 import org.forester.msa.MsaInferrer;
+import org.forester.msa.MsaMethods;
 import org.forester.pccx.TestPccx;
 import org.forester.phylogeny.Phylogeny;
 import org.forester.phylogeny.PhylogenyBranch;
 import org.forester.phylogeny.PhylogenyMethods;
 import org.forester.phylogeny.PhylogenyNode;
+import org.forester.phylogeny.PhylogenyNodeI.NH_CONVERSION_SUPPORT_VALUE_STYLE;
 import org.forester.phylogeny.data.BinaryCharacters;
 import org.forester.phylogeny.data.BranchWidth;
 import org.forester.phylogeny.data.Confidence;
@@ -70,6 +73,7 @@ import org.forester.phylogeny.data.DomainArchitecture;
 import org.forester.phylogeny.data.Event;
 import org.forester.phylogeny.data.Identifier;
 import org.forester.phylogeny.data.PhylogenyData;
+import org.forester.phylogeny.data.PhylogenyDataUtil;
 import org.forester.phylogeny.data.Polygon;
 import org.forester.phylogeny.data.PropertiesMap;
 import org.forester.phylogeny.data.Property;
@@ -79,14 +83,13 @@ import org.forester.phylogeny.data.Taxonomy;
 import org.forester.phylogeny.factories.ParserBasedPhylogenyFactory;
 import org.forester.phylogeny.factories.PhylogenyFactory;
 import org.forester.phylogeny.iterators.PhylogenyNodeIterator;
+import org.forester.protein.Protein;
 import org.forester.sdi.SDI;
 import org.forester.sdi.SDIR;
 import org.forester.sdi.SDIse;
-import org.forester.sdi.TaxonomyAssigner;
 import org.forester.sdi.TestGSDI;
 import org.forester.sequence.BasicSequence;
 import org.forester.sequence.Sequence;
-import org.forester.surfacing.Protein;
 import org.forester.surfacing.TestSurfacing;
 import org.forester.tools.ConfidenceAssessor;
 import org.forester.tools.SupportCount;
@@ -99,10 +102,10 @@ import org.forester.util.DescriptiveStatistics;
 import org.forester.util.ForesterConstants;
 import org.forester.util.ForesterUtil;
 import org.forester.util.GeneralTable;
-import org.forester.ws.uniprot.DatabaseTools;
-import org.forester.ws.uniprot.SequenceDatabaseEntry;
-import org.forester.ws.uniprot.UniProtTaxonomy;
-import org.forester.ws.uniprot.UniProtWsTools;
+import org.forester.util.SequenceIdParser;
+import org.forester.ws.seqdb.SequenceDatabaseEntry;
+import org.forester.ws.seqdb.SequenceDbWsTools;
+import org.forester.ws.seqdb.UniProtTaxonomy;
 import org.forester.ws.wabi.TxSearch;
 import org.forester.ws.wabi.TxSearch.RANKS;
 import org.forester.ws.wabi.TxSearch.TAX_NAME_CLASS;
@@ -167,6 +170,16 @@ public final class Test {
             System.exit( -1 );
         }
         final long start_time = new Date().getTime();
+        System.out.print( "Sequence id parsing: " );
+        if ( testSequenceIdParsing() ) {
+            System.out.println( "OK." );
+            succeeded++;
+        }
+        else {
+            System.out.println( "failed." );
+            System.exit( -1 ); //TODO FIXME remove me!! ~
+            failed++;
+        }
         System.out.print( "Hmmscan output parser: " );
         if ( testHmmscanOutputParser() ) {
             System.out.println( "OK." );
@@ -230,6 +243,15 @@ public final class Test {
             System.out.println( "failed." );
             failed++;
         }
+        System.out.print( "NHX parsing (MrBayes): " );
+        if ( Test.testNHXParsingMB() ) {
+            System.out.println( "OK." );
+            succeeded++;
+        }
+        else {
+            System.out.println( "failed." );
+            failed++;
+        }
         System.out.print( "Nexus characters parsing: " );
         if ( Test.testNexusCharactersParsing() ) {
             System.out.println( "OK." );
@@ -473,15 +495,6 @@ public final class Test {
             System.out.println( "failed." );
             failed++;
         }
-        System.out.print( "Taxonomy assigner: " );
-        if ( Test.testTaxonomyAssigner() ) {
-            System.out.println( "OK." );
-            succeeded++;
-        }
-        else {
-            System.out.println( "failed." );
-            failed++;
-        }
         System.out.print( "SDIunrooted: " );
         if ( Test.testSDIunrooted() ) {
             System.out.println( "OK." );
@@ -647,7 +660,6 @@ public final class Test {
             System.out.println( "failed." );
             failed++;
         }
-        
         System.out.print( "EMBL Entry Retrieval: " );
         if ( Test.testEmblEntryRetrieval() ) {
             System.out.println( "OK." );
@@ -657,7 +669,6 @@ public final class Test {
             System.out.println( "failed." );
             failed++;
         }
-        
         System.out.print( "Uniprot Entry Retrieval: " );
         if ( Test.testUniprotEntryRetrieval() ) {
             System.out.println( "OK." );
@@ -676,9 +687,27 @@ public final class Test {
             System.out.println( "failed." );
             failed++;
         }
-        if ( Mafft.isInstalled() ) {
+        //----
+        String path = "";
+        final String os = ForesterUtil.OS_NAME.toLowerCase();
+        if ( ( os.indexOf( "mac" ) >= 0 ) && ( os.indexOf( "os" ) > 0 ) ) {
+            path = "/usr/local/bin/mafft";
+        }
+        else if ( os.indexOf( "win" ) >= 0 ) {
+            path = "C:\\Program Files\\mafft-win\\mafft.bat";
+        }
+        else {
+            path = "/home/czmasek/bin/mafft";
+        }
+        if ( !Mafft.isInstalled( path ) ) {
+            path = "mafft";
+        }
+        if ( !Mafft.isInstalled( path ) ) {
+            path = "/usr/local/bin/mafft";
+        }
+        if ( Mafft.isInstalled( path ) ) {
             System.out.print( "MAFFT (external program): " );
-            if ( Test.testMafft() ) {
+            if ( Test.testMafft( path ) ) {
                 System.out.println( "OK." );
                 succeeded++;
             }
@@ -686,6 +715,25 @@ public final class Test {
                 System.out.println( "failed [will not count towards failed tests]" );
             }
         }
+        //----
+        System.out.print( "Next nodes with collapsed: " );
+        if ( Test.testNextNodeWithCollapsing() ) {
+            System.out.println( "OK." );
+            succeeded++;
+        }
+        else {
+            System.out.println( "failed." );
+            failed++;
+        }
+        System.out.print( "Simple MSA quality: " );
+        if ( Test.testMsaQualityMethod() ) {
+            System.out.println( "OK." );
+            succeeded++;
+        }
+        else {
+            System.out.println( "failed." );
+            failed++;
+        }
         //        System.out.print( "WABI TxSearch: " );
         //        if ( Test.testWabiTxSearch() ) {
         //            System.out.println( "OK." );
@@ -786,9 +834,12 @@ public final class Test {
                 return false;
             }
             final PhylogenyNode n1 = new PhylogenyNode();
-            final PhylogenyNode n2 = new PhylogenyNode( "", ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
-            final PhylogenyNode n3 = new PhylogenyNode( "n3", ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
-            final PhylogenyNode n4 = new PhylogenyNode( "n4:0.01", ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            final PhylogenyNode n2 = PhylogenyNode
+                    .createInstanceFromNhxString( "", PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            final PhylogenyNode n3 = PhylogenyNode
+                    .createInstanceFromNhxString( "n3", PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            final PhylogenyNode n4 = PhylogenyNode
+                    .createInstanceFromNhxString( "n4:0.01", PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
             if ( n1.isHasAssignedEvent() ) {
                 return false;
             }
@@ -1682,6 +1733,7 @@ public final class Test {
             final List<BasicTable<String>> tl = BasicTableParser.parse( source2.toString(),
                                                                         ";",
                                                                         false,
+                                                                        false,
                                                                         "comment:",
                                                                         false );
             if ( tl.size() != 2 ) {
@@ -1882,26 +1934,6 @@ public final class Test {
                 return false;
             }
             final PhylogenyNode n = t3.getNode( "ABC" );
-            PhylogenyNodeIterator it;
-            for( it = n.iterateChildNodesForward(); it.hasNext(); ) {
-                it.next();
-            }
-            for( it.reset(); it.hasNext(); ) {
-                it.next();
-            }
-            final PhylogenyNodeIterator it2 = n.iterateChildNodesForward();
-            if ( !it2.next().getName().equals( "A" ) ) {
-                return false;
-            }
-            if ( !it2.next().getName().equals( "B" ) ) {
-                return false;
-            }
-            if ( !it2.next().getName().equals( "C" ) ) {
-                return false;
-            }
-            if ( it2.hasNext() ) {
-                return false;
-            }
             final Phylogeny t4 = factory.create( "((A:1,B:2,C:10)ABC:1,(D:3,E:5)DE:3,(F,G,H,I))", new NHXParser() )[ 0 ];
             if ( t4.getNumberOfExternalNodes() != 9 ) {
                 return false;
@@ -2050,7 +2082,8 @@ public final class Test {
 
     private static boolean testCopyOfNodeData() {
         try {
-            final PhylogenyNode n1 = new PhylogenyNode( "n5:0.1[&&NHX:S=Ecoli:E=1.1.1.1:D=Y:Co=Y:B=56:T=1:O=22:SO=33:SN=44:W=2:C=10.20.30:XN=S=tag1=value1=unit1]" );
+            final PhylogenyNode n1 = PhylogenyNode
+                    .createInstanceFromNhxString( "n5:0.1[&&NHX:S=Ecoli:E=1.1.1.1:D=Y:Co=Y:B=56:T=1:O=22:SO=33:SN=44:W=2:C=10.20.30:XN=S=tag1=value1=unit1]" );
             final PhylogenyNode n2 = n1.copyNodeData();
             if ( !n1.toNewHampshireX().equals( n2.toNewHampshireX() ) ) {
                 return false;
@@ -2100,7 +2133,7 @@ public final class Test {
                 return false;
             }
             t2.setIdentifier( new Identifier( "ecoli" ) );
-            t2.setTaxonomyCode( "other" );
+            t2.setTaxonomyCode( "OTHER" );
             t2.setScientificName( "what" );
             t2.setCommonName( "something" );
             if ( !t1.isEqual( t2 ) ) {
@@ -2828,7 +2861,7 @@ public final class Test {
             dss3.addValue( 10 );
             final AsciiHistogram histo = new AsciiHistogram( dss3 );
             histo.toStringBuffer( 10, '=', 40, 5 );
-            histo.toStringBuffer( 3, 8, 10, '=', 40, 5 );
+            histo.toStringBuffer( 3, 8, 10, '=', 40, 5, null );
         }
         catch ( final Exception e ) {
             e.printStackTrace( System.out );
@@ -3178,200 +3211,199 @@ public final class Test {
             final PhylogenyFactory factory = ParserBasedPhylogenyFactory.getInstance();
             final Phylogeny p1 = factory.create( "((((((A,B)ab,C)abc,D)abcd,E)abcde,F)abcdef,(G,H)gh)abcdefgh",
                                                  new NHXParser() )[ 0 ];
-            final PhylogenyMethods pm = PhylogenyMethods.getInstance();
-            final PhylogenyNode A = pm.obtainLCA( p1.getNode( "A" ), p1.getNode( "A" ) );
+            final PhylogenyNode A = PhylogenyMethods.obtainLCA( p1.getNode( "A" ), p1.getNode( "A" ) );
             if ( !A.getName().equals( "A" ) ) {
                 return false;
             }
-            final PhylogenyNode gh = pm.obtainLCA( p1.getNode( "gh" ), p1.getNode( "gh" ) );
+            final PhylogenyNode gh = PhylogenyMethods.obtainLCA( p1.getNode( "gh" ), p1.getNode( "gh" ) );
             if ( !gh.getName().equals( "gh" ) ) {
                 return false;
             }
-            final PhylogenyNode ab = pm.obtainLCA( p1.getNode( "A" ), p1.getNode( "B" ) );
+            final PhylogenyNode ab = PhylogenyMethods.obtainLCA( p1.getNode( "A" ), p1.getNode( "B" ) );
             if ( !ab.getName().equals( "ab" ) ) {
                 return false;
             }
-            final PhylogenyNode ab2 = pm.obtainLCA( p1.getNode( "B" ), p1.getNode( "A" ) );
+            final PhylogenyNode ab2 = PhylogenyMethods.obtainLCA( p1.getNode( "B" ), p1.getNode( "A" ) );
             if ( !ab2.getName().equals( "ab" ) ) {
                 return false;
             }
-            final PhylogenyNode gh2 = pm.obtainLCA( p1.getNode( "H" ), p1.getNode( "G" ) );
+            final PhylogenyNode gh2 = PhylogenyMethods.obtainLCA( p1.getNode( "H" ), p1.getNode( "G" ) );
             if ( !gh2.getName().equals( "gh" ) ) {
                 return false;
             }
-            final PhylogenyNode gh3 = pm.obtainLCA( p1.getNode( "G" ), p1.getNode( "H" ) );
+            final PhylogenyNode gh3 = PhylogenyMethods.obtainLCA( p1.getNode( "G" ), p1.getNode( "H" ) );
             if ( !gh3.getName().equals( "gh" ) ) {
                 return false;
             }
-            final PhylogenyNode abc = pm.obtainLCA( p1.getNode( "C" ), p1.getNode( "A" ) );
+            final PhylogenyNode abc = PhylogenyMethods.obtainLCA( p1.getNode( "C" ), p1.getNode( "A" ) );
             if ( !abc.getName().equals( "abc" ) ) {
                 return false;
             }
-            final PhylogenyNode abc2 = pm.obtainLCA( p1.getNode( "A" ), p1.getNode( "C" ) );
+            final PhylogenyNode abc2 = PhylogenyMethods.obtainLCA( p1.getNode( "A" ), p1.getNode( "C" ) );
             if ( !abc2.getName().equals( "abc" ) ) {
                 return false;
             }
-            final PhylogenyNode abcd = pm.obtainLCA( p1.getNode( "A" ), p1.getNode( "D" ) );
+            final PhylogenyNode abcd = PhylogenyMethods.obtainLCA( p1.getNode( "A" ), p1.getNode( "D" ) );
             if ( !abcd.getName().equals( "abcd" ) ) {
                 return false;
             }
-            final PhylogenyNode abcd2 = pm.obtainLCA( p1.getNode( "D" ), p1.getNode( "A" ) );
+            final PhylogenyNode abcd2 = PhylogenyMethods.obtainLCA( p1.getNode( "D" ), p1.getNode( "A" ) );
             if ( !abcd2.getName().equals( "abcd" ) ) {
                 return false;
             }
-            final PhylogenyNode abcdef = pm.obtainLCA( p1.getNode( "A" ), p1.getNode( "F" ) );
+            final PhylogenyNode abcdef = PhylogenyMethods.obtainLCA( p1.getNode( "A" ), p1.getNode( "F" ) );
             if ( !abcdef.getName().equals( "abcdef" ) ) {
                 return false;
             }
-            final PhylogenyNode abcdef2 = pm.obtainLCA( p1.getNode( "F" ), p1.getNode( "A" ) );
+            final PhylogenyNode abcdef2 = PhylogenyMethods.obtainLCA( p1.getNode( "F" ), p1.getNode( "A" ) );
             if ( !abcdef2.getName().equals( "abcdef" ) ) {
                 return false;
             }
-            final PhylogenyNode abcdef3 = pm.obtainLCA( p1.getNode( "ab" ), p1.getNode( "F" ) );
+            final PhylogenyNode abcdef3 = PhylogenyMethods.obtainLCA( p1.getNode( "ab" ), p1.getNode( "F" ) );
             if ( !abcdef3.getName().equals( "abcdef" ) ) {
                 return false;
             }
-            final PhylogenyNode abcdef4 = pm.obtainLCA( p1.getNode( "F" ), p1.getNode( "ab" ) );
+            final PhylogenyNode abcdef4 = PhylogenyMethods.obtainLCA( p1.getNode( "F" ), p1.getNode( "ab" ) );
             if ( !abcdef4.getName().equals( "abcdef" ) ) {
                 return false;
             }
-            final PhylogenyNode abcde = pm.obtainLCA( p1.getNode( "A" ), p1.getNode( "E" ) );
+            final PhylogenyNode abcde = PhylogenyMethods.obtainLCA( p1.getNode( "A" ), p1.getNode( "E" ) );
             if ( !abcde.getName().equals( "abcde" ) ) {
                 return false;
             }
-            final PhylogenyNode abcde2 = pm.obtainLCA( p1.getNode( "E" ), p1.getNode( "A" ) );
+            final PhylogenyNode abcde2 = PhylogenyMethods.obtainLCA( p1.getNode( "E" ), p1.getNode( "A" ) );
             if ( !abcde2.getName().equals( "abcde" ) ) {
                 return false;
             }
-            final PhylogenyNode r = pm.obtainLCA( p1.getNode( "abcdefgh" ), p1.getNode( "abcdefgh" ) );
+            final PhylogenyNode r = PhylogenyMethods.obtainLCA( p1.getNode( "abcdefgh" ), p1.getNode( "abcdefgh" ) );
             if ( !r.getName().equals( "abcdefgh" ) ) {
                 return false;
             }
-            final PhylogenyNode r2 = pm.obtainLCA( p1.getNode( "A" ), p1.getNode( "H" ) );
+            final PhylogenyNode r2 = PhylogenyMethods.obtainLCA( p1.getNode( "A" ), p1.getNode( "H" ) );
             if ( !r2.getName().equals( "abcdefgh" ) ) {
                 return false;
             }
-            final PhylogenyNode r3 = pm.obtainLCA( p1.getNode( "H" ), p1.getNode( "A" ) );
+            final PhylogenyNode r3 = PhylogenyMethods.obtainLCA( p1.getNode( "H" ), p1.getNode( "A" ) );
             if ( !r3.getName().equals( "abcdefgh" ) ) {
                 return false;
             }
-            final PhylogenyNode abcde3 = pm.obtainLCA( p1.getNode( "E" ), p1.getNode( "abcde" ) );
+            final PhylogenyNode abcde3 = PhylogenyMethods.obtainLCA( p1.getNode( "E" ), p1.getNode( "abcde" ) );
             if ( !abcde3.getName().equals( "abcde" ) ) {
                 return false;
             }
-            final PhylogenyNode abcde4 = pm.obtainLCA( p1.getNode( "abcde" ), p1.getNode( "E" ) );
+            final PhylogenyNode abcde4 = PhylogenyMethods.obtainLCA( p1.getNode( "abcde" ), p1.getNode( "E" ) );
             if ( !abcde4.getName().equals( "abcde" ) ) {
                 return false;
             }
-            final PhylogenyNode ab3 = pm.obtainLCA( p1.getNode( "ab" ), p1.getNode( "B" ) );
+            final PhylogenyNode ab3 = PhylogenyMethods.obtainLCA( p1.getNode( "ab" ), p1.getNode( "B" ) );
             if ( !ab3.getName().equals( "ab" ) ) {
                 return false;
             }
-            final PhylogenyNode ab4 = pm.obtainLCA( p1.getNode( "B" ), p1.getNode( "ab" ) );
+            final PhylogenyNode ab4 = PhylogenyMethods.obtainLCA( p1.getNode( "B" ), p1.getNode( "ab" ) );
             if ( !ab4.getName().equals( "ab" ) ) {
                 return false;
             }
             final Phylogeny p2 = factory.create( "(a,b,(((c,d)cd,e)cde,f)cdef)r", new NHXParser() )[ 0 ];
-            final PhylogenyNode cd = pm.obtainLCA( p2.getNode( "c" ), p2.getNode( "d" ) );
+            final PhylogenyNode cd = PhylogenyMethods.obtainLCA( p2.getNode( "c" ), p2.getNode( "d" ) );
             if ( !cd.getName().equals( "cd" ) ) {
                 return false;
             }
-            final PhylogenyNode cd2 = pm.obtainLCA( p2.getNode( "d" ), p2.getNode( "c" ) );
+            final PhylogenyNode cd2 = PhylogenyMethods.obtainLCA( p2.getNode( "d" ), p2.getNode( "c" ) );
             if ( !cd2.getName().equals( "cd" ) ) {
                 return false;
             }
-            final PhylogenyNode cde = pm.obtainLCA( p2.getNode( "c" ), p2.getNode( "e" ) );
+            final PhylogenyNode cde = PhylogenyMethods.obtainLCA( p2.getNode( "c" ), p2.getNode( "e" ) );
             if ( !cde.getName().equals( "cde" ) ) {
                 return false;
             }
-            final PhylogenyNode cde2 = pm.obtainLCA( p2.getNode( "e" ), p2.getNode( "c" ) );
+            final PhylogenyNode cde2 = PhylogenyMethods.obtainLCA( p2.getNode( "e" ), p2.getNode( "c" ) );
             if ( !cde2.getName().equals( "cde" ) ) {
                 return false;
             }
-            final PhylogenyNode cdef = pm.obtainLCA( p2.getNode( "c" ), p2.getNode( "f" ) );
+            final PhylogenyNode cdef = PhylogenyMethods.obtainLCA( p2.getNode( "c" ), p2.getNode( "f" ) );
             if ( !cdef.getName().equals( "cdef" ) ) {
                 return false;
             }
-            final PhylogenyNode cdef2 = pm.obtainLCA( p2.getNode( "d" ), p2.getNode( "f" ) );
+            final PhylogenyNode cdef2 = PhylogenyMethods.obtainLCA( p2.getNode( "d" ), p2.getNode( "f" ) );
             if ( !cdef2.getName().equals( "cdef" ) ) {
                 return false;
             }
-            final PhylogenyNode cdef3 = pm.obtainLCA( p2.getNode( "f" ), p2.getNode( "d" ) );
+            final PhylogenyNode cdef3 = PhylogenyMethods.obtainLCA( p2.getNode( "f" ), p2.getNode( "d" ) );
             if ( !cdef3.getName().equals( "cdef" ) ) {
                 return false;
             }
-            final PhylogenyNode rt = pm.obtainLCA( p2.getNode( "c" ), p2.getNode( "a" ) );
+            final PhylogenyNode rt = PhylogenyMethods.obtainLCA( p2.getNode( "c" ), p2.getNode( "a" ) );
             if ( !rt.getName().equals( "r" ) ) {
                 return false;
             }
             final Phylogeny p3 = factory
                     .create( "((((a,(b,c)bc)abc,(d,e)de)abcde,f)abcdef,(((g,h)gh,(i,j)ij)ghij,k)ghijk,l)",
                              new NHXParser() )[ 0 ];
-            final PhylogenyNode bc_3 = pm.obtainLCA( p3.getNode( "b" ), p3.getNode( "c" ) );
+            final PhylogenyNode bc_3 = PhylogenyMethods.obtainLCA( p3.getNode( "b" ), p3.getNode( "c" ) );
             if ( !bc_3.getName().equals( "bc" ) ) {
                 return false;
             }
-            final PhylogenyNode ac_3 = pm.obtainLCA( p3.getNode( "a" ), p3.getNode( "c" ) );
+            final PhylogenyNode ac_3 = PhylogenyMethods.obtainLCA( p3.getNode( "a" ), p3.getNode( "c" ) );
             if ( !ac_3.getName().equals( "abc" ) ) {
                 return false;
             }
-            final PhylogenyNode ad_3 = pm.obtainLCA( p3.getNode( "a" ), p3.getNode( "d" ) );
+            final PhylogenyNode ad_3 = PhylogenyMethods.obtainLCA( p3.getNode( "a" ), p3.getNode( "d" ) );
             if ( !ad_3.getName().equals( "abcde" ) ) {
                 return false;
             }
-            final PhylogenyNode af_3 = pm.obtainLCA( p3.getNode( "a" ), p3.getNode( "f" ) );
+            final PhylogenyNode af_3 = PhylogenyMethods.obtainLCA( p3.getNode( "a" ), p3.getNode( "f" ) );
             if ( !af_3.getName().equals( "abcdef" ) ) {
                 return false;
             }
-            final PhylogenyNode ag_3 = pm.obtainLCA( p3.getNode( "a" ), p3.getNode( "g" ) );
+            final PhylogenyNode ag_3 = PhylogenyMethods.obtainLCA( p3.getNode( "a" ), p3.getNode( "g" ) );
             if ( !ag_3.getName().equals( "" ) ) {
                 return false;
             }
             if ( !ag_3.isRoot() ) {
                 return false;
             }
-            final PhylogenyNode al_3 = pm.obtainLCA( p3.getNode( "a" ), p3.getNode( "l" ) );
+            final PhylogenyNode al_3 = PhylogenyMethods.obtainLCA( p3.getNode( "a" ), p3.getNode( "l" ) );
             if ( !al_3.getName().equals( "" ) ) {
                 return false;
             }
             if ( !al_3.isRoot() ) {
                 return false;
             }
-            final PhylogenyNode kl_3 = pm.obtainLCA( p3.getNode( "k" ), p3.getNode( "l" ) );
+            final PhylogenyNode kl_3 = PhylogenyMethods.obtainLCA( p3.getNode( "k" ), p3.getNode( "l" ) );
             if ( !kl_3.getName().equals( "" ) ) {
                 return false;
             }
             if ( !kl_3.isRoot() ) {
                 return false;
             }
-            final PhylogenyNode fl_3 = pm.obtainLCA( p3.getNode( "f" ), p3.getNode( "l" ) );
+            final PhylogenyNode fl_3 = PhylogenyMethods.obtainLCA( p3.getNode( "f" ), p3.getNode( "l" ) );
             if ( !fl_3.getName().equals( "" ) ) {
                 return false;
             }
             if ( !fl_3.isRoot() ) {
                 return false;
             }
-            final PhylogenyNode gk_3 = pm.obtainLCA( p3.getNode( "g" ), p3.getNode( "k" ) );
+            final PhylogenyNode gk_3 = PhylogenyMethods.obtainLCA( p3.getNode( "g" ), p3.getNode( "k" ) );
             if ( !gk_3.getName().equals( "ghijk" ) ) {
                 return false;
             }
             final Phylogeny p4 = factory.create( "(a,b,c)r", new NHXParser() )[ 0 ];
-            final PhylogenyNode r_4 = pm.obtainLCA( p4.getNode( "b" ), p4.getNode( "c" ) );
+            final PhylogenyNode r_4 = PhylogenyMethods.obtainLCA( p4.getNode( "b" ), p4.getNode( "c" ) );
             if ( !r_4.getName().equals( "r" ) ) {
                 return false;
             }
             final Phylogeny p5 = factory.create( "((a,b),c,d)root", new NHXParser() )[ 0 ];
-            final PhylogenyNode r_5 = pm.obtainLCA( p5.getNode( "a" ), p5.getNode( "c" ) );
+            final PhylogenyNode r_5 = PhylogenyMethods.obtainLCA( p5.getNode( "a" ), p5.getNode( "c" ) );
             if ( !r_5.getName().equals( "root" ) ) {
                 return false;
             }
             final Phylogeny p6 = factory.create( "((a,b),c,d)rot", new NHXParser() )[ 0 ];
-            final PhylogenyNode r_6 = pm.obtainLCA( p6.getNode( "c" ), p6.getNode( "a" ) );
+            final PhylogenyNode r_6 = PhylogenyMethods.obtainLCA( p6.getNode( "c" ), p6.getNode( "a" ) );
             if ( !r_6.getName().equals( "rot" ) ) {
                 return false;
             }
             final Phylogeny p7 = factory.create( "(((a,b)x,c)x,d,e)rott", new NHXParser() )[ 0 ];
-            final PhylogenyNode r_7 = pm.obtainLCA( p7.getNode( "a" ), p7.getNode( "e" ) );
+            final PhylogenyNode r_7 = PhylogenyMethods.obtainLCA( p7.getNode( "a" ), p7.getNode( "e" ) );
             if ( !r_7.getName().equals( "rott" ) ) {
                 return false;
             }
@@ -3391,11 +3423,11 @@ public final class Test {
             parser1.parse();
             final HmmscanPerDomainTableParser parser2 = new HmmscanPerDomainTableParser( new File( test_dir
                     + ForesterUtil.getFileSeparator() + "hmmscan30b3_output_2" ), "MONBR", INDIVIDUAL_SCORE_CUTOFF.NONE );
-            final List<Protein> domain_collections = parser2.parse();
+            final List<Protein> proteins = parser2.parse();
             if ( parser2.getProteinsEncountered() != 4 ) {
                 return false;
             }
-            if ( domain_collections.size() != 4 ) {
+            if ( proteins.size() != 4 ) {
                 return false;
             }
             if ( parser2.getDomainsEncountered() != 69 ) {
@@ -3407,11 +3439,25 @@ public final class Test {
             if ( parser2.getDomainsIgnoredDueToEval() != 0 ) {
                 return false;
             }
-            final Protein p1 = domain_collections.get( 0 );
+            final Protein p1 = proteins.get( 0 );
             if ( p1.getNumberOfProteinDomains() != 15 ) {
                 return false;
             }
-            final Protein p4 = domain_collections.get( 3 );
+            if ( p1.getLength() != 850 ) {
+                return false;
+            }
+            final Protein p2 = proteins.get( 1 );
+            if ( p2.getNumberOfProteinDomains() != 51 ) {
+                return false;
+            }
+            if ( p2.getLength() != 1291 ) {
+                return false;
+            }
+            final Protein p3 = proteins.get( 2 );
+            if ( p3.getNumberOfProteinDomains() != 2 ) {
+                return false;
+            }
+            final Protein p4 = proteins.get( 3 );
             if ( p4.getNumberOfProteinDomains() != 1 ) {
                 return false;
             }
@@ -4198,7 +4244,7 @@ public final class Test {
                 return false;
             }
             final NHXParser nhxp = new NHXParser();
-            nhxp.setTaxonomyExtraction( ForesterUtil.TAXONOMY_EXTRACTION.NO );
+            nhxp.setTaxonomyExtraction( PhylogenyMethods.TAXONOMY_EXTRACTION.NO );
             nhxp.setReplaceUnderscores( true );
             final Phylogeny uc0 = factory.create( "(A__A_,_B_B)", nhxp )[ 0 ];
             if ( !uc0.getRoot().getChildNode( 0 ).getName().equals( "A A " ) ) {
@@ -4424,6 +4470,58 @@ public final class Test {
             if ( ( p46.length != 1 ) || !p46[ 0 ].isEmpty() ) {
                 return false;
             }
+            final Phylogeny p47 = factory.create( new StringBuffer( "((A,B)ab:2[0.44],C)" ), new NHXParser() )[ 0 ];
+            if ( !isEqual( 0.44, p47.getNode( "ab" ).getBranchData().getConfidence( 0 ).getValue() ) ) {
+                return false;
+            }
+            final Phylogeny p48 = factory.create( new StringBuffer( "((A,B)ab:2[88],C)" ), new NHXParser() )[ 0 ];
+            if ( !isEqual( 88, p48.getNode( "ab" ).getBranchData().getConfidence( 0 ).getValue() ) ) {
+                return false;
+            }
+            final Phylogeny p49 = factory
+                    .create( new StringBuffer( "((A,B)a[comment:a,b;(a)]b:2[0.44][comment(a,b,b);],C)" ),
+                             new NHXParser() )[ 0 ];
+            if ( !isEqual( 0.44, p49.getNode( "ab" ).getBranchData().getConfidence( 0 ).getValue() ) ) {
+                return false;
+            }
+            final Phylogeny p50 = factory.create( new StringBuffer( "((\"A\",B)ab:2[88],C)" ), new NHXParser() )[ 0 ];
+            if ( p50.getNode( "A" ) == null ) {
+                return false;
+            }
+            if ( !p50.toNewHampshire( false, NH_CONVERSION_SUPPORT_VALUE_STYLE.IN_SQUARE_BRACKETS )
+                    .equals( "((A,B)ab:2.0[88],C);" ) ) {
+                return false;
+            }
+            if ( !p50.toNewHampshire( false, NH_CONVERSION_SUPPORT_VALUE_STYLE.NONE ).equals( "((A,B)ab:2.0,C);" ) ) {
+                return false;
+            }
+            if ( !p50.toNewHampshire( false, NH_CONVERSION_SUPPORT_VALUE_STYLE.AS_INTERNAL_NODE_NAMES )
+                    .equals( "((A,B)88:2.0,C);" ) ) {
+                return false;
+            }
+            final Phylogeny p51 = factory.create( new StringBuffer( "((\"A(A\",B)ab:2[88],C)" ), new NHXParser() )[ 0 ];
+            if ( p51.getNode( "A(A" ) == null ) {
+                return false;
+            }
+            final Phylogeny p52 = factory.create( new StringBuffer( "(('A(A',B)ab:2[88],C)" ), new NHXParser() )[ 0 ];
+            if ( p52.getNode( "A(A" ) == null ) {
+                return false;
+            }
+            final Phylogeny p53 = factory
+                    .create( new StringBuffer( "(('A(A',\"B (x (a' ,b) f(x);\"[com])[ment]ab:2[88],C)" ),
+                             new NHXParser() )[ 0 ];
+            if ( p53.getNode( "B (x (a' ,b) f(x);" ) == null ) {
+                return false;
+            }
+            // 
+            final Phylogeny p54 = factory.create( new StringBuffer( "((A,B):[88],C)" ), new NHXParser() )[ 0 ];
+            if ( p54.getNode( "A" ) == null ) {
+                return false;
+            }
+            if ( !p54.toNewHampshire( false, NH_CONVERSION_SUPPORT_VALUE_STYLE.IN_SQUARE_BRACKETS )
+                    .equals( "((A,B)[88],C);" ) ) {
+                return false;
+            }
         }
         catch ( final Exception e ) {
             e.printStackTrace( System.out );
@@ -4435,11 +4533,13 @@ public final class Test {
     private static boolean testNHXconversion() {
         try {
             final PhylogenyNode n1 = new PhylogenyNode();
-            final PhylogenyNode n2 = new PhylogenyNode( "" );
-            final PhylogenyNode n3 = new PhylogenyNode( "n3" );
-            final PhylogenyNode n4 = new PhylogenyNode( "n4:0.01" );
-            final PhylogenyNode n5 = new PhylogenyNode( "n5:0.1[&&NHX:S=Ecoli:E=1.1.1.1:D=Y:Co=Y:B=56:T=1:W=2:C=10.20.30:XN=S=tag1=value1=unit1]" );
-            final PhylogenyNode n6 = new PhylogenyNode( "n6:0.000001[&&NHX:S=Ecoli:E=1.1.1.1:D=N:Co=N:B=100:T=1:W=2:C=0.0.0:XN=B=bool_tag=T]" );
+            final PhylogenyNode n2 = PhylogenyNode.createInstanceFromNhxString( "" );
+            final PhylogenyNode n3 = PhylogenyNode.createInstanceFromNhxString( "n3" );
+            final PhylogenyNode n4 = PhylogenyNode.createInstanceFromNhxString( "n4:0.01" );
+            final PhylogenyNode n5 = PhylogenyNode
+                    .createInstanceFromNhxString( "n5:0.1[&&NHX:S=Ecoli:E=1.1.1.1:D=Y:Co=Y:B=56:T=1:W=2:C=10.20.30:XN=S=tag1=value1=unit1]" );
+            final PhylogenyNode n6 = PhylogenyNode
+                    .createInstanceFromNhxString( "n6:0.000001[&&NHX:S=Ecoli:E=1.1.1.1:D=N:Co=N:B=100:T=1:W=2:C=0.0.0:XN=B=bool_tag=T]" );
             if ( !n1.toNewHampshireX().equals( "" ) ) {
                 return false;
             }
@@ -4453,11 +4553,10 @@ public final class Test {
                 return false;
             }
             if ( !n5.toNewHampshireX()
-                    .equals( "n5:0.1[&&NHX:T=1:S=Ecoli:D=Y:XN=S=tag1=value1=unit1:B=56.0:W=2.0:C=10.20.30]" ) ) {
+                    .equals( "n5:0.1[&&NHX:T=1:S=Ecoli:D=Y:XN=S=tag1=value1=unit1:B=56:W=2.0:C=10.20.30]" ) ) {
                 return false;
             }
-            if ( !n6.toNewHampshireX()
-                    .equals( "n6:1.0E-6[&&NHX:T=1:S=Ecoli:D=N:XN=B=bool_tag=T:B=100.0:W=2.0:C=0.0.0]" ) ) {
+            if ( !n6.toNewHampshireX().equals( "n6:1.0E-6[&&NHX:T=1:S=Ecoli:D=N:XN=B=bool_tag=T:B=100:W=2.0:C=0.0.0]" ) ) {
                 return false;
             }
         }
@@ -4471,14 +4570,15 @@ public final class Test {
     private static boolean testNHXNodeParsing() {
         try {
             final PhylogenyNode n1 = new PhylogenyNode();
-            final PhylogenyNode n2 = new PhylogenyNode( "" );
-            final PhylogenyNode n3 = new PhylogenyNode( "n3" );
-            final PhylogenyNode n4 = new PhylogenyNode( "n4:0.01" );
-            final PhylogenyNode n5 = new PhylogenyNode( "n5:0.1[&&NHX:S=Ecoli:E=1.1.1.1:D=Y:B=56:T=1:On=22:SOn=33:SNn=44:W=2:C=10.20.30:XN=S=tag1=value1=unit1:XN=S=tag3=value3=unit3]" );
+            final PhylogenyNode n2 = PhylogenyNode.createInstanceFromNhxString( "" );
+            final PhylogenyNode n3 = PhylogenyNode.createInstanceFromNhxString( "n3" );
+            final PhylogenyNode n4 = PhylogenyNode.createInstanceFromNhxString( "n4:0.01" );
+            final PhylogenyNode n5 = PhylogenyNode
+                    .createInstanceFromNhxString( "n5:0.1[&&NHX:S=Ecoli:E=1.1.1.1:D=Y:B=56:T=1:On=22:SOn=33:SNn=44:W=2:C=10.20.30:XN=S=tag1=value1=unit1:XN=S=tag3=value3=unit3]" );
             if ( !n3.getName().equals( "n3" ) ) {
                 return false;
             }
-            if ( n3.getDistanceToParent() != PhylogenyNode.DISTANCE_DEFAULT ) {
+            if ( n3.getDistanceToParent() != PhylogenyDataUtil.BRANCH_LENGTH_DEFAULT ) {
                 return false;
             }
             if ( n3.isDuplication() ) {
@@ -4520,136 +4620,180 @@ public final class Test {
             if ( n5.getNodeData().getProperties().getPropertyRefs().length != 2 ) {
                 return false;
             }
-            final PhylogenyNode n8 = new PhylogenyNode( "n8_ECOLI/12:0.01",
-                                                        ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            final PhylogenyNode n8 = PhylogenyNode
+                    .createInstanceFromNhxString( "n8_ECOLI/12:0.01",
+                                                  PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
             if ( !n8.getName().equals( "n8_ECOLI/12" ) ) {
                 return false;
             }
             if ( !PhylogenyMethods.getSpecies( n8 ).equals( "ECOLI" ) ) {
                 return false;
             }
-            final PhylogenyNode n9 = new PhylogenyNode( "n9_ECOLI/12=12:0.01",
-                                                        ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            final PhylogenyNode n9 = PhylogenyNode
+                    .createInstanceFromNhxString( "n9_ECOLI/12=12:0.01",
+                                                  PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
             if ( !n9.getName().equals( "n9_ECOLI/12=12" ) ) {
                 return false;
             }
             if ( !PhylogenyMethods.getSpecies( n9 ).equals( "ECOLI" ) ) {
                 return false;
             }
-            final PhylogenyNode n10 = new PhylogenyNode( "n10.ECOLI", ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            final PhylogenyNode n10 = PhylogenyNode
+                    .createInstanceFromNhxString( "n10.ECOLI", PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
             if ( !n10.getName().equals( "n10.ECOLI" ) ) {
                 return false;
             }
-            final PhylogenyNode n20 = new PhylogenyNode( "n20_ECOLI/1-2",
-                                                         ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            final PhylogenyNode n20 = PhylogenyNode
+                    .createInstanceFromNhxString( "n20_ECOLI/1-2", PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
             if ( !n20.getName().equals( "n20_ECOLI/1-2" ) ) {
                 return false;
             }
             if ( !PhylogenyMethods.getSpecies( n20 ).equals( "ECOLI" ) ) {
                 return false;
             }
-            final PhylogenyNode n20x = new PhylogenyNode( "n20_ECOL1/1-2", ForesterUtil.TAXONOMY_EXTRACTION.YES );
+            final PhylogenyNode n20x = PhylogenyNode
+                    .createInstanceFromNhxString( "n20_ECOL1/1-2", PhylogenyMethods.TAXONOMY_EXTRACTION.YES );
             if ( !n20x.getName().equals( "n20_ECOL1/1-2" ) ) {
                 return false;
             }
             if ( !PhylogenyMethods.getSpecies( n20x ).equals( "ECOL1" ) ) {
                 return false;
             }
-            final PhylogenyNode n20xx = new PhylogenyNode( "n20_eCOL1/1-2",
-                                                           ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            final PhylogenyNode n20xx = PhylogenyNode
+                    .createInstanceFromNhxString( "n20_eCOL1/1-2", PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
             if ( !n20xx.getName().equals( "n20_eCOL1/1-2" ) ) {
                 return false;
             }
             if ( PhylogenyMethods.getSpecies( n20xx ).length() > 0 ) {
                 return false;
             }
-            final PhylogenyNode n20xxx = new PhylogenyNode( "n20_ecoli/1-2",
-                                                            ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            final PhylogenyNode n20xxx = PhylogenyNode
+                    .createInstanceFromNhxString( "n20_ecoli/1-2", PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
             if ( !n20xxx.getName().equals( "n20_ecoli/1-2" ) ) {
                 return false;
             }
             if ( PhylogenyMethods.getSpecies( n20xxx ).length() > 0 ) {
                 return false;
             }
-            final PhylogenyNode n20xxxx = new PhylogenyNode( "n20_Ecoli/1-2",
-                                                             ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            final PhylogenyNode n20xxxx = PhylogenyNode
+                    .createInstanceFromNhxString( "n20_Ecoli/1-2", PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
             if ( !n20xxxx.getName().equals( "n20_Ecoli/1-2" ) ) {
                 return false;
             }
             if ( PhylogenyMethods.getSpecies( n20xxxx ).length() > 0 ) {
                 return false;
             }
-            final PhylogenyNode n21 = new PhylogenyNode( "n21_PIG", ForesterUtil.TAXONOMY_EXTRACTION.YES );
+            final PhylogenyNode n21 = PhylogenyNode
+                    .createInstanceFromNhxString( "n21_PIG", PhylogenyMethods.TAXONOMY_EXTRACTION.YES );
             if ( !n21.getName().equals( "n21_PIG" ) ) {
                 return false;
             }
             if ( !PhylogenyMethods.getSpecies( n21 ).equals( "PIG" ) ) {
                 return false;
             }
-            final PhylogenyNode n21x = new PhylogenyNode( "n21_PIG", ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            final PhylogenyNode n21x = PhylogenyNode
+                    .createInstanceFromNhxString( "n21_PIG", PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
             if ( !n21x.getName().equals( "n21_PIG" ) ) {
                 return false;
             }
             if ( PhylogenyMethods.getSpecies( n21x ).length() > 0 ) {
                 return false;
             }
-            final PhylogenyNode n22 = new PhylogenyNode( "n22/PIG", ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            final PhylogenyNode n22 = PhylogenyNode
+                    .createInstanceFromNhxString( "n22/PIG", PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
             if ( !n22.getName().equals( "n22/PIG" ) ) {
                 return false;
             }
             if ( PhylogenyMethods.getSpecies( n22 ).length() > 0 ) {
                 return false;
             }
-            final PhylogenyNode n23 = new PhylogenyNode( "n23/PIG_1", ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            final PhylogenyNode n23 = PhylogenyNode
+                    .createInstanceFromNhxString( "n23/PIG_1", PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
             if ( !n23.getName().equals( "n23/PIG_1" ) ) {
                 return false;
             }
             if ( PhylogenyMethods.getSpecies( n23 ).length() > 0 ) {
                 return false;
             }
-            if ( NHXParser.LIMIT_SPECIES_NAMES_TO_FIVE_CHARS ) {
-                final PhylogenyNode a = new PhylogenyNode( "n10_ECOLI/1-2",
-                                                           ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
-                if ( !a.getName().equals( "n10_ECOLI/1-2" ) ) {
-                    return false;
-                }
-                if ( !PhylogenyMethods.getSpecies( a ).equals( "ECOLI" ) ) {
-                    return false;
-                }
-                final PhylogenyNode b = new PhylogenyNode( "n10_ECOLI1/1-2",
-                                                           ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
-                if ( !b.getName().equals( "n10_ECOLI1/1-2" ) ) {
-                    return false;
-                }
-                if ( !PhylogenyMethods.getSpecies( b ).equals( "ECOLI" ) ) {
-                    return false;
-                }
-                final PhylogenyNode c = new PhylogenyNode( "n10_RATAF12/1000-2000",
-                                                           ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
-                if ( !c.getName().equals( "n10_RATAF12/1000-2000" ) ) {
-                    return false;
-                }
-                if ( !PhylogenyMethods.getSpecies( c ).equals( "RATAF" ) ) {
-                    return false;
-                }
-                final PhylogenyNode d = new PhylogenyNode( "n10_RAT1/1-2",
-                                                           ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
-                if ( !d.getName().equals( "n10_RAT1/1-2" ) ) {
-                    return false;
-                }
-                if ( !PhylogenyMethods.getSpecies( d ).equals( "RAT" ) ) {
-                    return false;
-                }
-                final PhylogenyNode e = new PhylogenyNode( "n10_RAT1", ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
-                if ( !e.getName().equals( "n10_RAT1" ) ) {
-                    return false;
-                }
-                if ( !ForesterUtil.isEmpty( PhylogenyMethods.getSpecies( e ) ) ) {
-                    return false;
-                }
+            final PhylogenyNode a = PhylogenyNode
+                    .createInstanceFromNhxString( "n10_ECOLI/1-2", PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            if ( !a.getName().equals( "n10_ECOLI/1-2" ) ) {
+                return false;
+            }
+            if ( !PhylogenyMethods.getSpecies( a ).equals( "ECOLI" ) ) {
+                return false;
+            }
+            final PhylogenyNode b = PhylogenyNode
+                    .createInstanceFromNhxString( "n10_ECOLI1/1-2",
+                                                  PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            if ( !b.getName().equals( "n10_ECOLI1/1-2" ) ) {
+                return false;
+            }
+            if ( !PhylogenyMethods.getSpecies( b ).equals( "ECOLI" ) ) {
+                return false;
+            }
+            final PhylogenyNode c = PhylogenyNode
+                    .createInstanceFromNhxString( "n10_RATAF12/1000-2000",
+                                                  PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            if ( !c.getName().equals( "n10_RATAF12/1000-2000" ) ) {
+                return false;
+            }
+            if ( !PhylogenyMethods.getSpecies( c ).equals( "RATAF" ) ) {
+                return false;
+            }
+            final PhylogenyNode c1 = PhylogenyNode
+                    .createInstanceFromNhxString( "n10_BOVIN_1/1000-2000",
+                                                  PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            if ( !c1.getName().equals( "n10_BOVIN_1/1000-2000" ) ) {
+                return false;
+            }
+            if ( !PhylogenyMethods.getSpecies( c1 ).equals( "BOVIN" ) ) {
+                return false;
+            }
+            final PhylogenyNode c2 = PhylogenyNode
+                    .createInstanceFromNhxString( "n10_Bovin_1/1000-2000",
+                                                  PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            if ( !c2.getName().equals( "n10_Bovin_1/1000-2000" ) ) {
+                return false;
+            }
+            if ( !PhylogenyMethods.getSpecies( c2 ).equals( "" ) ) {
+                return false;
+            }
+            final PhylogenyNode d = PhylogenyNode
+                    .createInstanceFromNhxString( "n10_RAT1/1-2", PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            if ( !d.getName().equals( "n10_RAT1/1-2" ) ) {
+                return false;
             }
-            final PhylogenyNode n11 = new PhylogenyNode( "n111111_ECOLI/jdj:0.4",
-                                                         ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            if ( !PhylogenyMethods.getSpecies( d ).equals( "RAT" ) ) {
+                return false;
+            }
+            final PhylogenyNode e = PhylogenyNode
+                    .createInstanceFromNhxString( "n10_RAT1", PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            if ( !e.getName().equals( "n10_RAT1" ) ) {
+                return false;
+            }
+            if ( !ForesterUtil.isEmpty( PhylogenyMethods.getSpecies( e ) ) ) {
+                return false;
+            }
+            final PhylogenyNode e2 = PhylogenyNode
+                    .createInstanceFromNhxString( "n10_RAT1", PhylogenyMethods.TAXONOMY_EXTRACTION.YES );
+            if ( !e2.getName().equals( "n10_RAT1" ) ) {
+                return false;
+            }
+            if ( !PhylogenyMethods.getSpecies( e2 ).equals( "RAT" ) ) {
+                return false;
+            }
+            final PhylogenyNode e3 = PhylogenyNode
+                    .createInstanceFromNhxString( "n10_RAT~", PhylogenyMethods.TAXONOMY_EXTRACTION.YES );
+            if ( !e3.getName().equals( "n10_RAT~" ) ) {
+                return false;
+            }
+            if ( !PhylogenyMethods.getSpecies( e3 ).equals( "RAT" ) ) {
+                return false;
+            }
+            final PhylogenyNode n11 = PhylogenyNode
+                    .createInstanceFromNhxString( "n111111_ECOLI/jdj:0.4",
+                                                  PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
             if ( !n11.getName().equals( "n111111_ECOLI/jdj" ) ) {
                 return false;
             }
@@ -4659,8 +4803,9 @@ public final class Test {
             if ( !PhylogenyMethods.getSpecies( n11 ).equals( "ECOLI" ) ) {
                 return false;
             }
-            final PhylogenyNode n12 = new PhylogenyNode( "n111111-ECOLI---/jdj:0.4",
-                                                         ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            final PhylogenyNode n12 = PhylogenyNode
+                    .createInstanceFromNhxString( "n111111-ECOLI---/jdj:0.4",
+                                                  PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
             if ( !n12.getName().equals( "n111111-ECOLI---/jdj" ) ) {
                 return false;
             }
@@ -4670,6 +4815,22 @@ public final class Test {
             if ( PhylogenyMethods.getSpecies( n12 ).length() > 0 ) {
                 return false;
             }
+            final PhylogenyNode m = PhylogenyNode
+                    .createInstanceFromNhxString( "n10_MOUSEa", PhylogenyMethods.TAXONOMY_EXTRACTION.YES );
+            if ( !m.getName().equals( "n10_MOUSEa" ) ) {
+                return false;
+            }
+            if ( !PhylogenyMethods.getSpecies( m ).equals( "MOUSE" ) ) {
+                return false;
+            }
+            final PhylogenyNode o = PhylogenyNode
+                    .createInstanceFromNhxString( "n10_MOUSE_", PhylogenyMethods.TAXONOMY_EXTRACTION.YES );
+            if ( !o.getName().equals( "n10_MOUSE_" ) ) {
+                return false;
+            }
+            if ( !PhylogenyMethods.getSpecies( o ).equals( "MOUSE" ) ) {
+                return false;
+            }
             final Property tvu1 = n5.getNodeData().getProperties().getProperty( "tag1" );
             final Property tvu3 = n5.getNodeData().getProperties().getProperty( "tag3" );
             if ( !tvu1.getRef().equals( "tag1" ) ) {
@@ -4702,7 +4863,7 @@ public final class Test {
             if ( PhylogenyMethods.getConfidenceValue( n1 ) != Confidence.CONFIDENCE_DEFAULT_VALUE ) {
                 return false;
             }
-            if ( n1.getDistanceToParent() != PhylogenyNode.DISTANCE_DEFAULT ) {
+            if ( n1.getDistanceToParent() != PhylogenyDataUtil.BRANCH_LENGTH_DEFAULT ) {
                 return false;
             }
             if ( n2.getName().compareTo( "" ) != 0 ) {
@@ -4711,10 +4872,11 @@ public final class Test {
             if ( PhylogenyMethods.getConfidenceValue( n2 ) != Confidence.CONFIDENCE_DEFAULT_VALUE ) {
                 return false;
             }
-            if ( n2.getDistanceToParent() != PhylogenyNode.DISTANCE_DEFAULT ) {
+            if ( n2.getDistanceToParent() != PhylogenyDataUtil.BRANCH_LENGTH_DEFAULT ) {
                 return false;
             }
-            final PhylogenyNode n00 = new PhylogenyNode( "n7:0.000001[&&NHX:GN=gene_name:AC=accession123:ID=node_identifier:S=Ecoli:D=N:Co=N:B=100:T=1:On=100:SOn=100:SNn=100:W=2:C=0.0.0:XN=U=url_tag=www.yahoo.com]" );
+            final PhylogenyNode n00 = PhylogenyNode
+                    .createInstanceFromNhxString( "n7:0.000001[&&NHX:GN=gene_name:AC=accession123:ID=node_identifier:S=Ecoli:D=N:Co=N:B=100:T=1:On=100:SOn=100:SNn=100:W=2:C=0.0.0:XN=U=url_tag=www.yahoo.com]" );
             if ( !n00.getNodeData().getNodeIdentifier().getValue().equals( "node_identifier" ) ) {
                 return false;
             }
@@ -4739,32 +4901,35 @@ public final class Test {
             if ( !n00.getNodeData().getProperties().getProperty( "url_tag" ).getUnit().equals( "" ) ) {
                 return false;
             }
-            final PhylogenyNode nx = new PhylogenyNode( "n5:0.1[&&NHX:S=Ecoli:GN=gene_1]" );
+            final PhylogenyNode nx = PhylogenyNode.createInstanceFromNhxString( "n5:0.1[&&NHX:S=Ecoli:GN=gene_1]" );
             if ( !nx.getNodeData().getSequence().getName().equals( "gene_1" ) ) {
                 return false;
             }
-            final PhylogenyNode nx2 = new PhylogenyNode( "n5:0.1[&&NHX:S=Ecoli:G=gene_2]" );
+            final PhylogenyNode nx2 = PhylogenyNode.createInstanceFromNhxString( "n5:0.1[&&NHX:S=Ecoli:G=gene_2]" );
             if ( !nx2.getNodeData().getSequence().getName().equals( "gene_2" ) ) {
                 return false;
             }
-            final PhylogenyNode n13 = new PhylogenyNode( "blah_12345/1-2",
-                                                         ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            final PhylogenyNode n13 = PhylogenyNode
+                    .createInstanceFromNhxString( "blah_12345/1-2",
+                                                  PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
             if ( !n13.getName().equals( "blah_12345/1-2" ) ) {
                 return false;
             }
-            if ( !PhylogenyMethods.getSpecies( n13 ).equals( "" ) ) {
+            if ( !PhylogenyMethods.getSpecies( n13 ).equals( "12345" ) ) {
                 return false;
             }
-            final PhylogenyNode n14 = new PhylogenyNode( "blah_12X45/1-2",
-                                                         ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            final PhylogenyNode n14 = PhylogenyNode
+                    .createInstanceFromNhxString( "blah_12X45/1-2",
+                                                  PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
             if ( !n14.getName().equals( "blah_12X45/1-2" ) ) {
                 return false;
             }
             if ( !PhylogenyMethods.getSpecies( n14 ).equals( "12X45" ) ) {
                 return false;
             }
-            final PhylogenyNode n15 = new PhylogenyNode( "something_wicked[123]",
-                                                         ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            final PhylogenyNode n15 = PhylogenyNode
+                    .createInstanceFromNhxString( "something_wicked[123]",
+                                                  PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
             if ( !n15.getName().equals( "something_wicked" ) ) {
                 return false;
             }
@@ -4774,8 +4939,9 @@ public final class Test {
             if ( !isEqual( n15.getBranchData().getConfidence( 0 ).getValue(), 123 ) ) {
                 return false;
             }
-            final PhylogenyNode n16 = new PhylogenyNode( "something_wicked2[9]",
-                                                         ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            final PhylogenyNode n16 = PhylogenyNode
+                    .createInstanceFromNhxString( "something_wicked2[9]",
+                                                  PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
             if ( !n16.getName().equals( "something_wicked2" ) ) {
                 return false;
             }
@@ -4785,15 +4951,17 @@ public final class Test {
             if ( !isEqual( n16.getBranchData().getConfidence( 0 ).getValue(), 9 ) ) {
                 return false;
             }
-            final PhylogenyNode n17 = new PhylogenyNode( "something_wicked3[a]",
-                                                         ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            final PhylogenyNode n17 = PhylogenyNode
+                    .createInstanceFromNhxString( "something_wicked3[a]",
+                                                  PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
             if ( !n17.getName().equals( "something_wicked3" ) ) {
                 return false;
             }
             if ( n17.getBranchData().getNumberOfConfidences() != 0 ) {
                 return false;
             }
-            final PhylogenyNode n18 = new PhylogenyNode( ":0.5[91]", ForesterUtil.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
+            final PhylogenyNode n18 = PhylogenyNode
+                    .createInstanceFromNhxString( ":0.5[91]", PhylogenyMethods.TAXONOMY_EXTRACTION.PFAM_STYLE_ONLY );
             if ( !isEqual( n18.getDistanceToParent(), 0.5 ) ) {
                 return false;
             }
@@ -4823,7 +4991,7 @@ public final class Test {
             if ( !p2[ 0 ].toNewHampshireX().equals( p2_S ) ) {
                 return false;
             }
-            final String p2b_S = "(((((((A:0.2[&NHX:S=qwerty]):0.2[&:S=uiop]):0.3[&NHX:S=asdf]):0.4[S=zxc]):0.5[]):0.6[&&NH:S=asd]):0.7[&&HX:S=za]):0.8[&&:S=zaq]";
+            final String p2b_S = "(((((((A:0.2[&NHX:S=qw,erty]):0.2[&:S=u(io)p]):0.3[&NHX:S=asdf]):0.4[S=zxc]):0.5[]):0.6[&&NH:S=asd]):0.7[&&HX:S=za]):0.8[&&:S=zaq]";
             final Phylogeny[] p2b = factory.create( p2b_S, new NHXParser() );
             if ( !p2b[ 0 ].toNewHampshireX().equals( "(((((((A:0.2):0.2):0.3):0.4):0.5):0.6):0.7):0.8" ) ) {
                 return false;
@@ -4865,13 +5033,13 @@ public final class Test {
                 return false;
             }
             final Phylogeny p9 = factory.create( "((A:0.2,B:0.3):0.5[91],C:0.1)root:0.1[100]", new NHXParser() )[ 0 ];
-            if ( !p9.toNewHampshireX().equals( "((A:0.2,B:0.3):0.5[&&NHX:B=91.0],C:0.1)root:0.1[&&NHX:B=100.0]" ) ) {
+            if ( !p9.toNewHampshireX().equals( "((A:0.2,B:0.3):0.5[&&NHX:B=91],C:0.1)root:0.1[&&NHX:B=100]" ) ) {
                 return false;
             }
             final Phylogeny p10 = factory
                     .create( " [79]   ( (A [co mment] :0 .2[comment],B:0.3[com])[com ment]: 0. 5 \t[ 9 1 ][ comment],C: 0.1)[comment]root:0.1[100] [comment]",
                              new NHXParser() )[ 0 ];
-            if ( !p10.toNewHampshireX().equals( "((A:0.2,B:0.3):0.5[&&NHX:B=91.0],C:0.1)root:0.1[&&NHX:B=100.0]" ) ) {
+            if ( !p10.toNewHampshireX().equals( "((A:0.2,B:0.3):0.5[&&NHX:B=91],C:0.1)root:0.1[&&NHX:B=100]" ) ) {
                 return false;
             }
         }
@@ -4946,7 +5114,7 @@ public final class Test {
             final Phylogeny p10 = factory
                     .create( " [79]   ( (\"A \n\tB \" [co mment] :0 .2[comment],'B':0.3[com])[com ment]: 0. 5 \t[ 9 1 ][ comment],'C (or D?\\//;,))': 0.1)[comment]'\nroot is here (cool,  was! ) ':0.1[100] [comment]",
                              new NHXParser() )[ 0 ];
-            final String p10_clean_str = "(('A B':0.2,B:0.3):0.5[&&NHX:B=91.0],'C (or D?\\//;,))':0.1)'root is here (cool,  was! )':0.1[&&NHX:B=100.0]";
+            final String p10_clean_str = "(('A B':0.2,B:0.3):0.5[&&NHX:B=91],'C (or D?\\//;,))':0.1)'root is here (cool,  was! )':0.1[&&NHX:B=100]";
             if ( !p10.toNewHampshireX().equals( p10_clean_str ) ) {
                 return false;
             }
@@ -4958,7 +5126,7 @@ public final class Test {
             final Phylogeny p12 = factory
                     .create( " [79]   ( (\"A \n\tB \" [[][] :0 .2[comment][\t&\t&\n N\tH\tX:S=mo\tnkey !],'\tB\t\b\t\n\f\rB B ':0.0\b3[])\t[com ment]: 0. 5 \t[ 9 1 ][ \ncomment],'C\t (or D?\\//;,))': 0.\b1)[comment]'\nroot \tis here (cool, \b\t\n\f\r was! ) ':0.1[100] [comment]",
                              new NHXParser() )[ 0 ];
-            final String p12_clean_str = "(('A B':0.2[&&NHX:S=monkey!],'BB B':0.03):0.5[&&NHX:B=91.0],'C (or D?\\//;,))':0.1)'root is here (cool,  was! )':0.1[&&NHX:B=100.0]";
+            final String p12_clean_str = "(('A B':0.2[&&NHX:S=monkey!],'BB B':0.03):0.5[&&NHX:B=91],'C (or D?\\//;,))':0.1)'root is here (cool,  was! )':0.1[&&NHX:B=100]";
             if ( !p12.toNewHampshireX().equals( p12_clean_str ) ) {
                 return false;
             }
@@ -4982,10 +5150,64 @@ public final class Test {
         return true;
     }
 
+    private static boolean testNHXParsingMB() {
+        try {
+            final PhylogenyFactory factory = ParserBasedPhylogenyFactory.getInstance();
+            final Phylogeny p1 = factory.create( "(1[&prob=0.9500000000000000e+00,prob_stddev=0.1100000000000000e+00,"
+                    + "prob_range={1.000000000000000e+00,1.000000000000000e+00},prob(percent)=\"100\","
+                    + "prob+-sd=\"100+-0\"]:4.129000000000000e-02[&length_mean=4.153987461671767e-02,"
+                    + "length_median=4.129000000000000e-02,length_95%HPD={3.217800000000000e-02,"
+                    + "5.026800000000000e-02}],2[&prob=0.810000000000000e+00,prob_stddev=0.000000000000000e+00,"
+                    + "prob_range={1.000000000000000e+00,1.000000000000000e+00},prob(percent)=\"100\","
+                    + "prob+-sd=\"100+-0\"]:6.375699999999999e-02[&length_mean=6.395210411945065e-02,"
+                    + "length_median=6.375699999999999e-02,length_95%HPD={5.388600000000000e-02,"
+                    + "7.369400000000000e-02}])", new NHXParser() )[ 0 ];
+            if ( !isEqual( p1.getNode( "1" ).getDistanceToParent(), 4.129e-02 ) ) {
+                return false;
+            }
+            if ( !isEqual( p1.getNode( "1" ).getBranchData().getConfidence( 0 ).getValue(), 0.9500000000000000e+00 ) ) {
+                return false;
+            }
+            if ( !isEqual( p1.getNode( "1" ).getBranchData().getConfidence( 0 ).getStandardDeviation(),
+                           0.1100000000000000e+00 ) ) {
+                return false;
+            }
+            if ( !isEqual( p1.getNode( "2" ).getDistanceToParent(), 6.375699999999999e-02 ) ) {
+                return false;
+            }
+            if ( !isEqual( p1.getNode( "2" ).getBranchData().getConfidence( 0 ).getValue(), 0.810000000000000e+00 ) ) {
+                return false;
+            }
+            final Phylogeny p2 = factory
+                    .create( "(1[something_else(?)s,prob=0.9500000000000000e+00{}(((,p)rob_stddev=0.110000000000e+00,"
+                                     + "prob_range={1.000000000000000e+00,1.000000000000000e+00},prob(percent)=\"100\","
+                                     + "prob+-sd=\"100+-0\"]:4.129000000000000e-02[&length_mean=4.153987461671767e-02,"
+                                     + "length_median=4.129000000000000e-02,length_95%HPD={3.217800000000000e-02,"
+                                     + "5.026800000000000e-02}],2[&prob=0.810000000000000e+00,prob_stddev=0.000000000000000e+00,"
+                                     + "prob_range={1.000000000000000e+00,1.000000000000000e+00},prob(percent)=\"100\","
+                                     + "prob+-sd=\"100+-0\"]:6.375699999999999e-02[&length_mean=6.395210411945065e-02,"
+                                     + "length_median=6.375699999999999e-02,length_95%HPD={5.388600000000000e-02,"
+                                     + "7.369400000000000e-02}])",
+                             new NHXParser() )[ 0 ];
+            if ( p2.getNode( "1" ) == null ) {
+                return false;
+            }
+            if ( p2.getNode( "2" ) == null ) {
+                return false;
+            }
+        }
+        catch ( final Exception e ) {
+            e.printStackTrace( System.out );
+            System.exit( -1 );
+            return false;
+        }
+        return true;
+    }
+
     private static boolean testPhylogenyBranch() {
         try {
-            final PhylogenyNode a1 = new PhylogenyNode( "a" );
-            final PhylogenyNode b1 = new PhylogenyNode( "b" );
+            final PhylogenyNode a1 = PhylogenyNode.createInstanceFromNhxString( "a" );
+            final PhylogenyNode b1 = PhylogenyNode.createInstanceFromNhxString( "b" );
             final PhylogenyBranch a1b1 = new PhylogenyBranch( a1, b1 );
             final PhylogenyBranch b1a1 = new PhylogenyBranch( b1, a1 );
             if ( !a1b1.equals( a1b1 ) ) {
@@ -6423,209 +6645,209 @@ public final class Test {
             final Phylogeny p0 = factory.create( "(((A,B,C),D),(E,(F,G)))R", new NHXParser() )[ 0 ];
             //Archaeopteryx.createApplication( p0 );
             final Set<PhylogenyNode> ex = new HashSet<PhylogenyNode>();
-            ex.add( new PhylogenyNode( "A" ) );
-            ex.add( new PhylogenyNode( "B" ) );
-            ex.add( new PhylogenyNode( "C" ) );
-            ex.add( new PhylogenyNode( "D" ) );
-            ex.add( new PhylogenyNode( "E" ) );
-            ex.add( new PhylogenyNode( "F" ) );
-            ex.add( new PhylogenyNode( "G" ) );
-            ex.add( new PhylogenyNode( "X" ) );
-            ex.add( new PhylogenyNode( "Y" ) );
+            ex.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            ex.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
+            ex.add( PhylogenyNode.createInstanceFromNhxString( "C" ) );
+            ex.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
+            ex.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            ex.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            ex.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
+            ex.add( PhylogenyNode.createInstanceFromNhxString( "X" ) );
+            ex.add( PhylogenyNode.createInstanceFromNhxString( "Y" ) );
             final TreeSplitMatrix s0 = new TreeSplitMatrix( p0, false, ex );
             // System.out.println( s0.toString() );
             //
             Set<PhylogenyNode> query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
-            query_nodes.add( new PhylogenyNode( "C" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "F" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "C" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
             if ( !s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
-            query_nodes.add( new PhylogenyNode( "C" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "C" ) );
             if ( !s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "D" ) );
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "F" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
             if ( !s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
-            query_nodes.add( new PhylogenyNode( "C" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "C" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
             if ( !s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "F" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
             if ( !s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "F" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
             if ( !s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
-            query_nodes.add( new PhylogenyNode( "C" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
-            query_nodes.add( new PhylogenyNode( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "C" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
             if ( !s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "F" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
-            query_nodes.add( new PhylogenyNode( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
             if ( !s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "F" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
             if ( !s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "F" ) );
-            query_nodes.add( new PhylogenyNode( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
-            query_nodes.add( new PhylogenyNode( "C" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "C" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "F" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
-            query_nodes.add( new PhylogenyNode( "C" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "C" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "C" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "C" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "F" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
-            query_nodes.add( new PhylogenyNode( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
@@ -6695,114 +6917,114 @@ public final class Test {
             //            }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "X" ) );
-            query_nodes.add( new PhylogenyNode( "Y" ) );
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "X" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "Y" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "X" ) );
-            query_nodes.add( new PhylogenyNode( "Y" ) );
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "X" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "Y" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             ///////////////////////////
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "X" ) );
-            query_nodes.add( new PhylogenyNode( "Y" ) );
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "X" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "Y" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "X" ) );
-            query_nodes.add( new PhylogenyNode( "Y" ) );
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "X" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "Y" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "X" ) );
-            query_nodes.add( new PhylogenyNode( "Y" ) );
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "C" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "X" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "Y" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "C" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "X" ) );
-            query_nodes.add( new PhylogenyNode( "Y" ) );
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "X" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "Y" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "X" ) );
-            query_nodes.add( new PhylogenyNode( "Y" ) );
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "X" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "Y" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "Y" ) );
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "Y" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "X" ) );
-            query_nodes.add( new PhylogenyNode( "Y" ) );
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "F" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "X" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "Y" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "X" ) );
-            query_nodes.add( new PhylogenyNode( "Y" ) );
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "X" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "Y" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "X" ) );
-            query_nodes.add( new PhylogenyNode( "Y" ) );
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
-            query_nodes.add( new PhylogenyNode( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "X" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "Y" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "X" ) );
-            query_nodes.add( new PhylogenyNode( "Y" ) );
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "X" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "Y" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
@@ -6819,205 +7041,205 @@ public final class Test {
             final PhylogenyFactory factory = ParserBasedPhylogenyFactory.getInstance();
             final Phylogeny p0 = factory.create( "(((A,B,C),D),(E,(F,G)))R", new NHXParser() )[ 0 ];
             final Set<PhylogenyNode> ex = new HashSet<PhylogenyNode>();
-            ex.add( new PhylogenyNode( "A" ) );
-            ex.add( new PhylogenyNode( "B" ) );
-            ex.add( new PhylogenyNode( "C" ) );
-            ex.add( new PhylogenyNode( "D" ) );
-            ex.add( new PhylogenyNode( "E" ) );
-            ex.add( new PhylogenyNode( "F" ) );
-            ex.add( new PhylogenyNode( "G" ) );
+            ex.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            ex.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
+            ex.add( PhylogenyNode.createInstanceFromNhxString( "C" ) );
+            ex.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
+            ex.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            ex.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            ex.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
             final TreeSplitMatrix s0 = new TreeSplitMatrix( p0, true, ex );
             Set<PhylogenyNode> query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
-            query_nodes.add( new PhylogenyNode( "C" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "F" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "C" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
             if ( !s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
-            query_nodes.add( new PhylogenyNode( "C" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "C" ) );
             if ( !s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "D" ) );
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "F" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
             if ( !s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
-            query_nodes.add( new PhylogenyNode( "C" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "C" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
             if ( !s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "F" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
             if ( !s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "F" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
             if ( !s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
-            query_nodes.add( new PhylogenyNode( "C" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
-            query_nodes.add( new PhylogenyNode( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "C" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
             if ( !s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "F" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
-            query_nodes.add( new PhylogenyNode( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
             if ( !s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "F" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
             if ( !s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "F" ) );
-            query_nodes.add( new PhylogenyNode( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
-            query_nodes.add( new PhylogenyNode( "C" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "C" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "F" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
-            query_nodes.add( new PhylogenyNode( "C" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "C" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "C" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "C" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "F" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "F" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "B" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "B" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
-            query_nodes.add( new PhylogenyNode( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
             //
             query_nodes = new HashSet<PhylogenyNode>();
-            query_nodes.add( new PhylogenyNode( "E" ) );
-            query_nodes.add( new PhylogenyNode( "D" ) );
-            query_nodes.add( new PhylogenyNode( "A" ) );
-            query_nodes.add( new PhylogenyNode( "G" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "E" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "D" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "A" ) );
+            query_nodes.add( PhylogenyNode.createInstanceFromNhxString( "G" ) );
             if ( s0.match( query_nodes ) ) {
                 return false;
             }
@@ -7284,450 +7506,97 @@ public final class Test {
         return true;
     }
 
-    private static boolean testTaxonomyAssigner() {
+    private static boolean testUniprotTaxonomySearch() {
         try {
-            String s0_str = "(((([&&NHX:S=A],[&&NHX:S=B])[&&NHX:S=AB],[&&NHX:S=C])[&&NHX:S=ABC],[&&NHX:S=D])[&&NHX:S=ABCD],[&&NHX:S=E])[&&NHX:S=ABCDE]";
-            String g0_str = "((([&&NHX:S=A],[&&NHX:S=A],[&&NHX:S=A])a,[&&NHX:S=B])b,[&&NHX:S=C])c";
-            Phylogeny s0 = ParserBasedPhylogenyFactory.getInstance().create( s0_str, new NHXParser() )[ 0 ];
-            Phylogeny g0 = ParserBasedPhylogenyFactory.getInstance().create( g0_str, new NHXParser() )[ 0 ];
-            s0.setRooted( true );
-            g0.setRooted( true );
-            TaxonomyAssigner.execute( g0, s0 );
-            if ( !g0.getNode( "a" ).getNodeData().getTaxonomy().getScientificName().equals( "A" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "b" ).getNodeData().getTaxonomy().getScientificName().equals( "AB" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "c" ).getNodeData().getTaxonomy().getScientificName().equals( "ABC" ) ) {
-                return false;
-            }
-            g0_str = "((([&&NHX:S=A],[&&NHX:S=A],[&&NHX:S=A])a,[&&NHX:S=A])b,[&&NHX:S=A])c";
-            g0 = ParserBasedPhylogenyFactory.getInstance().create( g0_str, new NHXParser() )[ 0 ];
-            g0.setRooted( true );
-            TaxonomyAssigner.execute( g0, s0 );
-            if ( !g0.getNode( "a" ).getNodeData().getTaxonomy().getScientificName().equals( "A" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "b" ).getNodeData().getTaxonomy().getScientificName().equals( "A" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "c" ).getNodeData().getTaxonomy().getScientificName().equals( "A" ) ) {
-                return false;
-            }
-            g0_str = "((([&&NHX:S=A],[&&NHX:S=A],[&&NHX:S=B])a,[&&NHX:S=A])b,[&&NHX:S=A])c";
-            g0 = ParserBasedPhylogenyFactory.getInstance().create( g0_str, new NHXParser() )[ 0 ];
-            g0.setRooted( true );
-            TaxonomyAssigner.execute( g0, s0 );
-            if ( !g0.getNode( "a" ).getNodeData().getTaxonomy().getScientificName().equals( "AB" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "b" ).getNodeData().getTaxonomy().getScientificName().equals( "AB" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "c" ).getNodeData().getTaxonomy().getScientificName().equals( "AB" ) ) {
-                return false;
-            }
-            g0_str = "((([&&NHX:S=A],[&&NHX:S=A],[&&NHX:S=B])a,[&&NHX:S=C])b,[&&NHX:S=A])c";
-            g0 = ParserBasedPhylogenyFactory.getInstance().create( g0_str, new NHXParser() )[ 0 ];
-            g0.setRooted( true );
-            TaxonomyAssigner.execute( g0, s0 );
-            if ( !g0.getNode( "a" ).getNodeData().getTaxonomy().getScientificName().equals( "AB" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "b" ).getNodeData().getTaxonomy().getScientificName().equals( "ABC" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "c" ).getNodeData().getTaxonomy().getScientificName().equals( "ABC" ) ) {
-                return false;
-            }
-            g0_str = "((([&&NHX:S=A],[&&NHX:S=A],[&&NHX:S=B])a,[&&NHX:S=C])b,[&&NHX:S=D])c";
-            g0 = ParserBasedPhylogenyFactory.getInstance().create( g0_str, new NHXParser() )[ 0 ];
-            g0.setRooted( true );
-            TaxonomyAssigner.execute( g0, s0 );
-            if ( !g0.getNode( "a" ).getNodeData().getTaxonomy().getScientificName().equals( "AB" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "b" ).getNodeData().getTaxonomy().getScientificName().equals( "ABC" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "c" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCD" ) ) {
-                return false;
-            }
-            g0_str = "((([&&NHX:S=A],[&&NHX:S=A],[&&NHX:S=E])a,[&&NHX:S=C])b,[&&NHX:S=D])c";
-            g0 = ParserBasedPhylogenyFactory.getInstance().create( g0_str, new NHXParser() )[ 0 ];
-            g0.setRooted( true );
-            TaxonomyAssigner.execute( g0, s0 );
-            if ( !g0.getNode( "a" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCDE" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "b" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCDE" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "c" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCDE" ) ) {
-                return false;
-            }
-            g0_str = "((([&&NHX:S=A],[&&NHX:S=A],[&&NHX:S=E])a,[&&NHX:S=A])b,[&&NHX:S=A])c";
-            g0 = ParserBasedPhylogenyFactory.getInstance().create( g0_str, new NHXParser() )[ 0 ];
-            g0.setRooted( true );
-            TaxonomyAssigner.execute( g0, s0 );
-            if ( !g0.getNode( "a" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCDE" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "b" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCDE" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "c" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCDE" ) ) {
-                return false;
-            }
-            s0_str = "(([&&NHX:S=A],[&&NHX:S=B],[&&NHX:S=C],[&&NHX:S=D])[&&NHX:S=ABCD],"
-                    + "([&&NHX:S=E],[&&NHX:S=F],[&&NHX:S=G],[&&NHX:S=H])[&&NHX:S=EFGH],"
-                    + "([&&NHX:S=I],[&&NHX:S=J],[&&NHX:S=K],[&&NHX:S=L])[&&NHX:S=IJKL], "
-                    + "([&&NHX:S=M],[&&NHX:S=N],[&&NHX:S=O],[&&NHX:S=P])[&&NHX:S=MNOP])[&&NHX:S=ROOT]";
-            s0 = ParserBasedPhylogenyFactory.getInstance().create( s0_str, new NHXParser() )[ 0 ];
-            s0.setRooted( true );
-            g0_str = "(([&&NHX:S=A],[&&NHX:S=B],[&&NHX:S=C],[&&NHX:S=D])a,"
-                    + "([&&NHX:S=E],[&&NHX:S=F],[&&NHX:S=G],[&&NHX:S=H])b,"
-                    + "([&&NHX:S=I],[&&NHX:S=J],[&&NHX:S=K],[&&NHX:S=L])c, "
-                    + "([&&NHX:S=M],[&&NHX:S=N],[&&NHX:S=O],[&&NHX:S=P])d)r";
-            g0 = ParserBasedPhylogenyFactory.getInstance().create( g0_str, new NHXParser() )[ 0 ];
-            g0.setRooted( true );
-            TaxonomyAssigner.execute( g0, s0 );
-            if ( !g0.getNode( "a" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCD" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "b" ).getNodeData().getTaxonomy().getScientificName().equals( "EFGH" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "c" ).getNodeData().getTaxonomy().getScientificName().equals( "IJKL" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "d" ).getNodeData().getTaxonomy().getScientificName().equals( "MNOP" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "r" ).getNodeData().getTaxonomy().getScientificName().equals( "ROOT" ) ) {
+            List<UniProtTaxonomy> results = SequenceDbWsTools.getTaxonomiesFromCommonNameStrict( "starlet sea anemone",
+                                                                                                 10 );
+            if ( results.size() != 1 ) {
                 return false;
             }
-            g0_str = "(([&&NHX:S=A],[&&NHX:S=B],[&&NHX:S=A],[&&NHX:S=B])a,"
-                    + "([&&NHX:S=E],[&&NHX:S=F],[&&NHX:S=F],[&&NHX:S=F])b,"
-                    + "([&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=I])c, "
-                    + "([&&NHX:S=M],[&&NHX:S=N],[&&NHX:S=O],[&&NHX:S=O])d)r";
-            g0 = ParserBasedPhylogenyFactory.getInstance().create( g0_str, new NHXParser() )[ 0 ];
-            g0.setRooted( true );
-            TaxonomyAssigner.execute( g0, s0 );
-            if ( !g0.getNode( "a" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCD" ) ) {
+            if ( !results.get( 0 ).getCode().equals( "NEMVE" ) ) {
                 return false;
             }
-            if ( !g0.getNode( "b" ).getNodeData().getTaxonomy().getScientificName().equals( "EFGH" ) ) {
+            if ( !results.get( 0 ).getCommonName().equalsIgnoreCase( "starlet sea anemone" ) ) {
                 return false;
             }
-            if ( !g0.getNode( "c" ).getNodeData().getTaxonomy().getScientificName().equals( "IJKL" ) ) {
+            if ( !results.get( 0 ).getId().equalsIgnoreCase( "45351" ) ) {
                 return false;
             }
-            if ( !g0.getNode( "d" ).getNodeData().getTaxonomy().getScientificName().equals( "MNOP" ) ) {
+            if ( !results.get( 0 ).getRank().equalsIgnoreCase( "species" ) ) {
                 return false;
             }
-            if ( !g0.getNode( "r" ).getNodeData().getTaxonomy().getScientificName().equals( "ROOT" ) ) {
+            if ( !results.get( 0 ).getScientificName().equals( "Nematostella vectensis" ) ) {
                 return false;
             }
-            g0_str = "(([&&NHX:S=A],[&&NHX:S=B],[&&NHX:S=A],[&&NHX:S=B])a,"
-                    + "([&&NHX:S=E],[&&NHX:S=F],[&&NHX:S=F],[&&NHX:S=F])b,"
-                    + "([&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L])c, "
-                    + "([&&NHX:S=M],[&&NHX:S=N],[&&NHX:S=A],[&&NHX:S=O])d)r";
-            g0 = ParserBasedPhylogenyFactory.getInstance().create( g0_str, new NHXParser() )[ 0 ];
-            g0.setRooted( true );
-            TaxonomyAssigner.execute( g0, s0 );
-            if ( !g0.getNode( "a" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCD" ) ) {
+            results = null;
+            results = SequenceDbWsTools.getTaxonomiesFromScientificNameStrict( "Nematostella vectensis", 10 );
+            if ( results.size() != 1 ) {
                 return false;
             }
-            if ( !g0.getNode( "b" ).getNodeData().getTaxonomy().getScientificName().equals( "EFGH" ) ) {
+            if ( !results.get( 0 ).getCode().equals( "NEMVE" ) ) {
                 return false;
             }
-            if ( !g0.getNode( "c" ).getNodeData().getTaxonomy().getScientificName().equals( "L" ) ) {
+            if ( !results.get( 0 ).getCommonName().equalsIgnoreCase( "starlet sea anemone" ) ) {
                 return false;
             }
-            if ( !g0.getNode( "d" ).getNodeData().getTaxonomy().getScientificName().equals( "ROOT" ) ) {
+            if ( !results.get( 0 ).getId().equalsIgnoreCase( "45351" ) ) {
                 return false;
             }
-            if ( !g0.getNode( "r" ).getNodeData().getTaxonomy().getScientificName().equals( "ROOT" ) ) {
+            if ( !results.get( 0 ).getRank().equalsIgnoreCase( "species" ) ) {
                 return false;
             }
-            g0_str = "(([&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L])a,"
-                    + "([&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L])b,"
-                    + "([&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L])c, "
-                    + "([&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L])d)r";
-            g0 = ParserBasedPhylogenyFactory.getInstance().create( g0_str, new NHXParser() )[ 0 ];
-            g0.setRooted( true );
-            TaxonomyAssigner.execute( g0, s0 );
-            if ( !g0.getNode( "a" ).getNodeData().getTaxonomy().getScientificName().equals( "L" ) ) {
+            if ( !results.get( 0 ).getScientificName().equals( "Nematostella vectensis" ) ) {
                 return false;
             }
-            if ( !g0.getNode( "b" ).getNodeData().getTaxonomy().getScientificName().equals( "L" ) ) {
+            results = null;
+            results = SequenceDbWsTools.getTaxonomiesFromId( "45351", 10 );
+            if ( results.size() != 1 ) {
                 return false;
             }
-            if ( !g0.getNode( "c" ).getNodeData().getTaxonomy().getScientificName().equals( "L" ) ) {
+            if ( !results.get( 0 ).getCode().equals( "NEMVE" ) ) {
                 return false;
             }
-            if ( !g0.getNode( "d" ).getNodeData().getTaxonomy().getScientificName().equals( "L" ) ) {
+            if ( !results.get( 0 ).getCommonName().equalsIgnoreCase( "starlet sea anemone" ) ) {
                 return false;
             }
-            if ( !g0.getNode( "r" ).getNodeData().getTaxonomy().getScientificName().equals( "L" ) ) {
+            if ( !results.get( 0 ).getId().equalsIgnoreCase( "45351" ) ) {
                 return false;
             }
-            g0_str = "((([&&NHX:S=A],[&&NHX:S=A],[&&NHX:S=A])a,[&&NHX:S=A])b,[&&NHX:S=A])c";
-            g0 = ParserBasedPhylogenyFactory.getInstance().create( g0_str, new NHXParser() )[ 0 ];
-            g0.setRooted( true );
-            TaxonomyAssigner.execute( g0, s0 );
-            if ( !g0.getNode( "a" ).getNodeData().getTaxonomy().getScientificName().equals( "A" ) ) {
+            if ( !results.get( 0 ).getRank().equalsIgnoreCase( "species" ) ) {
                 return false;
             }
-            if ( !g0.getNode( "b" ).getNodeData().getTaxonomy().getScientificName().equals( "A" ) ) {
+            if ( !results.get( 0 ).getScientificName().equals( "Nematostella vectensis" ) ) {
                 return false;
             }
-            if ( !g0.getNode( "c" ).getNodeData().getTaxonomy().getScientificName().equals( "A" ) ) {
+            results = null;
+            results = SequenceDbWsTools.getTaxonomiesFromTaxonomyCode( "NEMVE", 10 );
+            if ( results.size() != 1 ) {
                 return false;
             }
-            g0_str = "((([&&NHX:S=A],[&&NHX:S=A],[&&NHX:S=B])a,[&&NHX:S=I])b,[&&NHX:S=J])c";
-            g0 = ParserBasedPhylogenyFactory.getInstance().create( g0_str, new NHXParser() )[ 0 ];
-            g0.setRooted( true );
-            TaxonomyAssigner.execute( g0, s0 );
-            if ( !g0.getNode( "a" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCD" ) ) {
+            if ( !results.get( 0 ).getCode().equals( "NEMVE" ) ) {
                 return false;
             }
-            if ( !g0.getNode( "b" ).getNodeData().getTaxonomy().getScientificName().equals( "ROOT" ) ) {
+            if ( !results.get( 0 ).getCommonName().equalsIgnoreCase( "starlet sea anemone" ) ) {
                 return false;
             }
-            if ( !g0.getNode( "c" ).getNodeData().getTaxonomy().getScientificName().equals( "ROOT" ) ) {
+            if ( !results.get( 0 ).getId().equalsIgnoreCase( "45351" ) ) {
                 return false;
             }
-            g0_str = "(((([&&NHX:S=A],[&&NHX:S=B],[&&NHX:S=C],[&&NHX:S=D])a,"
-                    + "([&&NHX:S=D],[&&NHX:S=C],[&&NHX:S=B],[&&NHX:S=A])b)ab,"
-                    + "([&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L])c)abc, "
-                    + "([&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L])d)r";
-            g0 = ParserBasedPhylogenyFactory.getInstance().create( g0_str, new NHXParser() )[ 0 ];
-            g0.setRooted( true );
-            TaxonomyAssigner.execute( g0, s0 );
-            if ( !g0.getNode( "a" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCD" ) ) {
+            if ( !results.get( 0 ).getRank().equalsIgnoreCase( "species" ) ) {
                 return false;
             }
-            if ( !g0.getNode( "b" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCD" ) ) {
+            if ( !results.get( 0 ).getScientificName().equals( "Nematostella vectensis" ) ) {
                 return false;
             }
-            if ( !g0.getNode( "ab" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCD" ) ) {
+            if ( !results.get( 0 ).getLineage().get( 1 ).equals( "Eukaryota" ) ) {
                 return false;
             }
-            if ( !g0.getNode( "c" ).getNodeData().getTaxonomy().getScientificName().equals( "L" ) ) {
+            if ( !results.get( 0 ).getLineage().get( 2 ).equals( "Metazoa" ) ) {
                 return false;
             }
-            if ( !g0.getNode( "abc" ).getNodeData().getTaxonomy().getScientificName().equals( "ROOT" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "d" ).getNodeData().getTaxonomy().getScientificName().equals( "L" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "r" ).getNodeData().getTaxonomy().getScientificName().equals( "ROOT" ) ) {
-                return false;
-            }
-            g0_str = "(((([&&NHX:S=A],[&&NHX:S=A],[&&NHX:S=C],[&&NHX:S=D])a,"
-                    + "([&&NHX:S=D],[&&NHX:S=D],[&&NHX:S=B],[&&NHX:S=A])b)ab,"
-                    + "([&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L])c)abc, "
-                    + "([&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L])d)r";
-            g0 = ParserBasedPhylogenyFactory.getInstance().create( g0_str, new NHXParser() )[ 0 ];
-            g0.setRooted( true );
-            TaxonomyAssigner.execute( g0, s0 );
-            if ( !g0.getNode( "a" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCD" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "b" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCD" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "ab" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCD" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "c" ).getNodeData().getTaxonomy().getScientificName().equals( "L" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "abc" ).getNodeData().getTaxonomy().getScientificName().equals( "ROOT" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "d" ).getNodeData().getTaxonomy().getScientificName().equals( "L" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "r" ).getNodeData().getTaxonomy().getScientificName().equals( "ROOT" ) ) {
-                return false;
-            }
-            g0_str = "(((([&&NHX:S=A],[&&NHX:S=A],[&&NHX:S=C],[&&NHX:S=D])a,"
-                    + "([&&NHX:S=D],[&&NHX:S=D],[&&NHX:S=B],[&&NHX:S=A])b)ab,"
-                    + "([&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L])c)abc, "
-                    + "([&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=A])d)r";
-            g0 = ParserBasedPhylogenyFactory.getInstance().create( g0_str, new NHXParser() )[ 0 ];
-            g0.setRooted( true );
-            TaxonomyAssigner.execute( g0, s0 );
-            if ( !g0.getNode( "a" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCD" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "b" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCD" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "ab" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCD" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "c" ).getNodeData().getTaxonomy().getScientificName().equals( "L" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "abc" ).getNodeData().getTaxonomy().getScientificName().equals( "ROOT" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "d" ).getNodeData().getTaxonomy().getScientificName().equals( "ROOT" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "r" ).getNodeData().getTaxonomy().getScientificName().equals( "ROOT" ) ) {
-                return false;
-            }
-            g0_str = "(((([&&NHX:S=A],[&&NHX:S=A],[&&NHX:S=C],[&&NHX:S=D])a,"
-                    + "([&&NHX:S=D],[&&NHX:S=D],[&&NHX:S=B],[&&NHX:S=A])b)ab,"
-                    + "([&&NHX:S=A],[&&NHX:S=A],[&&NHX:S=A],[&&NHX:S=A])c)abc, "
-                    + "([&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=A])d)r";
-            g0 = ParserBasedPhylogenyFactory.getInstance().create( g0_str, new NHXParser() )[ 0 ];
-            g0.setRooted( true );
-            TaxonomyAssigner.execute( g0, s0 );
-            if ( !g0.getNode( "a" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCD" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "b" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCD" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "ab" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCD" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "c" ).getNodeData().getTaxonomy().getScientificName().equals( "A" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "abc" ).getNodeData().getTaxonomy().getScientificName().equals( "ABCD" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "d" ).getNodeData().getTaxonomy().getScientificName().equals( "ROOT" ) ) {
-                return false;
-            }
-            if ( !g0.getNode( "r" ).getNodeData().getTaxonomy().getScientificName().equals( "ROOT" ) ) {
-                return false;
-            }
-            s0_str = "(([&&NHX:S=A],[&&NHX:S=B],[&&NHX:S=C],[&&NHX:S=D]),"
-                    + "([&&NHX:S=E],[&&NHX:S=F],[&&NHX:S=G],[&&NHX:S=H]),"
-                    + "([&&NHX:S=I],[&&NHX:S=J],[&&NHX:S=K],[&&NHX:S=L]), "
-                    + "([&&NHX:S=M],[&&NHX:S=N],[&&NHX:S=O],[&&NHX:S=P]))";
-            s0 = ParserBasedPhylogenyFactory.getInstance().create( s0_str, new NHXParser() )[ 0 ];
-            s0.setRooted( true );
-            g0_str = "(((([&&NHX:S=A],[&&NHX:S=A],[&&NHX:S=C],[&&NHX:S=D])a,"
-                    + "([&&NHX:S=D],[&&NHX:S=D],[&&NHX:S=B],[&&NHX:S=A])b)ab,"
-                    + "([&&NHX:S=A],[&&NHX:S=A],[&&NHX:S=A],[&&NHX:S=A])c)abc, "
-                    + "([&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=L],[&&NHX:S=A])d)r";
-            g0 = ParserBasedPhylogenyFactory.getInstance().create( g0_str, new NHXParser() )[ 0 ];
-            g0.setRooted( true );
-            TaxonomyAssigner.execute( g0, s0 );
-            if ( g0.getNode( "a" ).getNodeData().isHasTaxonomy() ) {
-                return false;
-            }
-            if ( !g0.getNode( "c" ).getNodeData().getTaxonomy().getScientificName().equals( "A" ) ) {
-                return false;
-            }
-        }
-        catch ( final Exception e ) {
-            e.printStackTrace( System.out );
-            return false;
-        }
-        return true;
-    }
-
-    private static boolean testUniprotTaxonomySearch() {
-        try {
-            List<UniProtTaxonomy> results = UniProtWsTools
-                    .getTaxonomiesFromCommonNameStrict( "starlet sea anemone", 10 );
-            if ( results.size() != 1 ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getCode().equals( "NEMVE" ) ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getCommonName().equalsIgnoreCase( "starlet sea anemone" ) ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getId().equalsIgnoreCase( "45351" ) ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getRank().equalsIgnoreCase( "species" ) ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getScientificName().equals( "Nematostella vectensis" ) ) {
-                return false;
-            }
-            results = null;
-            results = UniProtWsTools.getTaxonomiesFromScientificNameStrict( "Nematostella vectensis", 10 );
-            if ( results.size() != 1 ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getCode().equals( "NEMVE" ) ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getCommonName().equalsIgnoreCase( "starlet sea anemone" ) ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getId().equalsIgnoreCase( "45351" ) ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getRank().equalsIgnoreCase( "species" ) ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getScientificName().equals( "Nematostella vectensis" ) ) {
-                return false;
-            }
-            results = null;
-            results = UniProtWsTools.getTaxonomiesFromId( "45351", 10 );
-            if ( results.size() != 1 ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getCode().equals( "NEMVE" ) ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getCommonName().equalsIgnoreCase( "starlet sea anemone" ) ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getId().equalsIgnoreCase( "45351" ) ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getRank().equalsIgnoreCase( "species" ) ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getScientificName().equals( "Nematostella vectensis" ) ) {
-                return false;
-            }
-            results = null;
-            results = UniProtWsTools.getTaxonomiesFromTaxonomyCode( "NEMVE", 10 );
-            if ( results.size() != 1 ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getCode().equals( "NEMVE" ) ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getCommonName().equalsIgnoreCase( "starlet sea anemone" ) ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getId().equalsIgnoreCase( "45351" ) ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getRank().equalsIgnoreCase( "species" ) ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getScientificName().equals( "Nematostella vectensis" ) ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getLineage()[ 0 ].equals( "Eukaryota" ) ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getLineage()[ 1 ].equals( "Metazoa" ) ) {
-                return false;
-            }
-            if ( !results.get( 0 ).getLineage()[ results.get( 0 ).getLineage().length - 1 ].equals( "Nematostella" ) ) {
+            if ( !results.get( 0 ).getLineage().get( results.get( 0 ).getLineage().size() - 1 )
+                    .equals( "Nematostella vectensis" ) ) {
+                System.out.println( results.get( 0 ).getLineage() );
                 return false;
             }
         }
@@ -7744,37 +7613,94 @@ public final class Test {
     }
 
     private static boolean testEmblEntryRetrieval() {
-        if ( !DatabaseTools.parseGenbankAccessor( "AY423861" ).equals( "AY423861" ) ) {
-            System.out.println( DatabaseTools.parseGenbankAccessor( "AY423861" ));
+        //The format for GenBank Accession numbers are:
+        //Nucleotide: 1 letter + 5 numerals OR 2 letters + 6 numerals
+        //Protein:    3 letters + 5 numerals
+        //http://www.ncbi.nlm.nih.gov/Sequin/acc.html
+        if ( !SequenceIdParser.parseGenbankAccessor( "AY423861" ).equals( "AY423861" ) ) {
+            return false;
+        }
+        if ( !SequenceIdParser.parseGenbankAccessor( ".AY423861." ).equals( "AY423861" ) ) {
+            return false;
+        }
+        if ( SequenceIdParser.parseGenbankAccessor( "AAY423861" ) != null ) {
+            return false;
+        }
+        if ( SequenceIdParser.parseGenbankAccessor( "AY4238612" ) != null ) {
+            return false;
+        }
+        if ( SequenceIdParser.parseGenbankAccessor( "AAY4238612" ) != null ) {
+            return false;
+        }
+        if ( SequenceIdParser.parseGenbankAccessor( "Y423861" ) != null ) {
+            return false;
+        }
+        if ( !SequenceIdParser.parseGenbankAccessor( "S12345" ).equals( "S12345" ) ) {
+            return false;
+        }
+        if ( !SequenceIdParser.parseGenbankAccessor( "|S12345|" ).equals( "S12345" ) ) {
+            return false;
+        }
+        if ( SequenceIdParser.parseGenbankAccessor( "|S123456" ) != null ) {
+            return false;
+        }
+        if ( SequenceIdParser.parseGenbankAccessor( "ABC123456" ) != null ) {
+            return false;
+        }
+        if ( !SequenceIdParser.parseGenbankAccessor( "ABC12345" ).equals( "ABC12345" ) ) {
+            return false;
+        }
+        if ( !SequenceIdParser.parseGenbankAccessor( "&ABC12345&" ).equals( "ABC12345" ) ) {
+            return false;
+        }
+        if ( SequenceIdParser.parseGenbankAccessor( "ABCD12345" ) != null ) {
             return false;
         }
         return true;
     }
-    
+
     private static boolean testUniprotEntryRetrieval() {
-        if ( !UniProtWsTools.parseUniProtAccessor( "P12345" ).equals( "P12345" ) ) {
+        if ( !SequenceDbWsTools.parseUniProtAccessor( "P12345" ).equals( "P12345" ) ) {
+            return false;
+        }
+        if ( SequenceDbWsTools.parseUniProtAccessor( "EP12345" ) != null ) {
+            return false;
+        }
+        if ( SequenceDbWsTools.parseUniProtAccessor( "3 4P12345" ) != null ) {
+            return false;
+        }
+        if ( SequenceDbWsTools.parseUniProtAccessor( "P12345E" ) != null ) {
+            return false;
+        }
+        if ( SequenceDbWsTools.parseUniProtAccessor( "P123455" ) != null ) {
+            return false;
+        }
+        if ( SequenceDbWsTools.parseUniProtAccessor( "EP12345E" ) != null ) {
+            return false;
+        }
+        if ( SequenceDbWsTools.parseUniProtAccessor( "AY423861" ) != null ) {
             return false;
         }
-        if ( UniProtWsTools.parseUniProtAccessor( "EP12345" ) != null ) {
+        if ( !SequenceDbWsTools.parseUniProtAccessor( "P1DDD5" ).equals( "P1DDD5" ) ) {
             return false;
         }
-        if ( !UniProtWsTools.parseUniProtAccessor( "P1DDD5" ).equals( "P1DDD5" ) ) {
+        if ( SequenceDbWsTools.parseUniProtAccessor( "P1DDDD" ) != null ) {
             return false;
         }
-        if ( UniProtWsTools.parseUniProtAccessor( "P1DDDD" ) != null ) {
+        if ( !SequenceDbWsTools.parseUniProtAccessor( "P1234X/P12345/12-42" ).equals( "P12345" ) ) {
             return false;
         }
-        if ( !UniProtWsTools.parseUniProtAccessor( "P1234X/P12345/12-42" ).equals( "P12345" ) ) {
+        if ( !SequenceDbWsTools.parseUniProtAccessor( "P1234X P12345 12-42" ).equals( "P12345" ) ) {
             return false;
         }
-        if ( !UniProtWsTools.parseUniProtAccessor( "P12345/12-42" ).equals( "P12345" ) ) {
+        if ( !SequenceDbWsTools.parseUniProtAccessor( "P12345/12-42" ).equals( "P12345" ) ) {
             return false;
         }
-        if ( !UniProtWsTools.parseUniProtAccessor( "P1234X/P12345" ).equals( "P12345" ) ) {
+        if ( !SequenceDbWsTools.parseUniProtAccessor( "P1234X/P12345" ).equals( "P12345" ) ) {
             return false;
         }
         try {
-            final SequenceDatabaseEntry entry = UniProtWsTools.obtainUniProtEntry( "P12345", 200 );
+            final SequenceDatabaseEntry entry = SequenceDbWsTools.obtainUniProtEntry( "P12345", 200 );
             if ( !entry.getAccession().equals( "P12345" ) ) {
                 return false;
             }
@@ -7957,12 +7883,48 @@ public final class Test {
         try {
             final String msa_str_0 = "seq1 abcd\n\nseq2 efgh\n";
             final Msa msa_0 = GeneralMsaParser.parse( new ByteArrayInputStream( msa_str_0.getBytes() ) );
-            final String msa_str_1 = "seq_1 abc\nseq2 ghi\nseq_1 def\nseq2 jkm\n";
+            final String msa_str_1 = "seq1 abc\nseq2 ghi\nseq1 def\nseq2 jkm\n";
             final Msa msa_1 = GeneralMsaParser.parse( new ByteArrayInputStream( msa_str_1.getBytes() ) );
             final String msa_str_2 = "seq1 abc\nseq2 ghi\n\ndef\njkm\n";
             final Msa msa_2 = GeneralMsaParser.parse( new ByteArrayInputStream( msa_str_2.getBytes() ) );
             final String msa_str_3 = "seq1 abc\n def\nseq2 ghi\n jkm\n";
             final Msa msa_3 = GeneralMsaParser.parse( new ByteArrayInputStream( msa_str_3.getBytes() ) );
+            if ( !msa_1.getSequenceAsString( 0 ).toString().equalsIgnoreCase( "abcdef" ) ) {
+                return false;
+            }
+            if ( !msa_1.getSequenceAsString( 1 ).toString().equalsIgnoreCase( "ghixkm" ) ) {
+                return false;
+            }
+            if ( !msa_1.getIdentifier( 0 ).toString().equals( "seq1" ) ) {
+                return false;
+            }
+            if ( !msa_1.getIdentifier( 1 ).toString().equals( "seq2" ) ) {
+                return false;
+            }
+            if ( !msa_2.getSequenceAsString( 0 ).toString().equalsIgnoreCase( "abcdef" ) ) {
+                return false;
+            }
+            if ( !msa_2.getSequenceAsString( 1 ).toString().equalsIgnoreCase( "ghixkm" ) ) {
+                return false;
+            }
+            if ( !msa_2.getIdentifier( 0 ).toString().equals( "seq1" ) ) {
+                return false;
+            }
+            if ( !msa_2.getIdentifier( 1 ).toString().equals( "seq2" ) ) {
+                return false;
+            }
+            if ( !msa_3.getSequenceAsString( 0 ).toString().equalsIgnoreCase( "abcdef" ) ) {
+                return false;
+            }
+            if ( !msa_3.getSequenceAsString( 1 ).toString().equalsIgnoreCase( "ghixkm" ) ) {
+                return false;
+            }
+            if ( !msa_3.getIdentifier( 0 ).toString().equals( "seq1" ) ) {
+                return false;
+            }
+            if ( !msa_3.getIdentifier( 1 ).toString().equals( "seq2" ) ) {
+                return false;
+            }
             final Msa msa_4 = GeneralMsaParser.parse( new FileInputStream( PATH_TO_TEST_DATA + "msa_1.txt" ) );
             if ( !msa_4.getSequenceAsString( 0 ).toString().equalsIgnoreCase( "abcdefeeeeeeeexx" ) ) {
                 return false;
@@ -8001,7 +7963,7 @@ public final class Test {
         return true;
     }
 
-    private static boolean testMafft() {
+    private static boolean testMafft( final String path ) {
         try {
             final List<String> opts = new ArrayList<String>();
             opts.add( "--maxiterate" );
@@ -8009,11 +7971,712 @@ public final class Test {
             opts.add( "--localpair" );
             opts.add( "--quiet" );
             Msa msa = null;
-            final MsaInferrer mafft = Mafft.createInstance();
-            msa = mafft.infer( new File( PATH_TO_TEST_DATA + "ncbi.fasta" ), opts );
-            if ( ( msa == null ) || ( msa.getLength() < 10 ) || ( msa.getNumberOfSequences() != 19 ) ) {
+            final MsaInferrer mafft = Mafft.createInstance( path );
+            msa = mafft.infer( new File( PATH_TO_TEST_DATA + "ncbi_sn.fasta" ), opts );
+            if ( ( msa == null ) || ( msa.getLength() < 20 ) || ( msa.getNumberOfSequences() != 19 ) ) {
+                return false;
+            }
+            if ( !msa.getIdentifier( 0 ).toString().equals( "a" ) ) {
+                return false;
+            }
+        }
+        catch ( final Exception e ) {
+            e.printStackTrace( System.out );
+            return false;
+        }
+        return true;
+    }
+
+    private static boolean testNextNodeWithCollapsing() {
+        try {
+            final PhylogenyFactory factory = ParserBasedPhylogenyFactory.getInstance();
+            PhylogenyNode n;
+            List<PhylogenyNode> ext = new ArrayList<PhylogenyNode>();
+            final StringBuffer sb0 = new StringBuffer( "((a,b)ab,(((c,d)cd,e)cde,(f,(g,h))fgh)cdefgh)abcdefgh" );
+            final Phylogeny t0 = factory.create( sb0, new NHXParser() )[ 0 ];
+            t0.getNode( "cd" ).setCollapse( true );
+            t0.getNode( "cde" ).setCollapse( true );
+            n = t0.getFirstExternalNode();
+            while ( n != null ) {
+                ext.add( n );
+                n = n.getNextExternalNodeWhileTakingIntoAccountCollapsedNodes();
+            }
+            if ( !ext.get( 0 ).getName().equals( "a" ) ) {
+                return false;
+            }
+            if ( !ext.get( 1 ).getName().equals( "b" ) ) {
+                return false;
+            }
+            if ( !ext.get( 2 ).getName().equals( "cde" ) ) {
+                return false;
+            }
+            if ( !ext.get( 3 ).getName().equals( "f" ) ) {
+                return false;
+            }
+            if ( !ext.get( 4 ).getName().equals( "g" ) ) {
+                return false;
+            }
+            if ( !ext.get( 5 ).getName().equals( "h" ) ) {
+                return false;
+            }
+            ext.clear();
+            final StringBuffer sb1 = new StringBuffer( "((a,b)ab,(((c,d)cd,e)cde,(f,(g,h))fgh)cdefgh)abcdefgh" );
+            final Phylogeny t1 = factory.create( sb1, new NHXParser() )[ 0 ];
+            t1.getNode( "ab" ).setCollapse( true );
+            t1.getNode( "cd" ).setCollapse( true );
+            t1.getNode( "cde" ).setCollapse( true );
+            n = t1.getNode( "ab" );
+            ext = new ArrayList<PhylogenyNode>();
+            while ( n != null ) {
+                ext.add( n );
+                n = n.getNextExternalNodeWhileTakingIntoAccountCollapsedNodes();
+            }
+            if ( !ext.get( 0 ).getName().equals( "ab" ) ) {
+                return false;
+            }
+            if ( !ext.get( 1 ).getName().equals( "cde" ) ) {
+                return false;
+            }
+            if ( !ext.get( 2 ).getName().equals( "f" ) ) {
+                return false;
+            }
+            if ( !ext.get( 3 ).getName().equals( "g" ) ) {
+                return false;
+            }
+            if ( !ext.get( 4 ).getName().equals( "h" ) ) {
+                return false;
+            }
+            //
+            //
+            ext.clear();
+            final StringBuffer sb2 = new StringBuffer( "((a,b)ab,(((c,d)cd,e)cde,(f,(g,h)gh)fgh)cdefgh)abcdefgh" );
+            final Phylogeny t2 = factory.create( sb2, new NHXParser() )[ 0 ];
+            t2.getNode( "ab" ).setCollapse( true );
+            t2.getNode( "cd" ).setCollapse( true );
+            t2.getNode( "cde" ).setCollapse( true );
+            t2.getNode( "c" ).setCollapse( true );
+            t2.getNode( "d" ).setCollapse( true );
+            t2.getNode( "e" ).setCollapse( true );
+            t2.getNode( "gh" ).setCollapse( true );
+            n = t2.getNode( "ab" );
+            ext = new ArrayList<PhylogenyNode>();
+            while ( n != null ) {
+                ext.add( n );
+                n = n.getNextExternalNodeWhileTakingIntoAccountCollapsedNodes();
+            }
+            if ( !ext.get( 0 ).getName().equals( "ab" ) ) {
+                return false;
+            }
+            if ( !ext.get( 1 ).getName().equals( "cde" ) ) {
+                return false;
+            }
+            if ( !ext.get( 2 ).getName().equals( "f" ) ) {
+                return false;
+            }
+            if ( !ext.get( 3 ).getName().equals( "gh" ) ) {
+                return false;
+            }
+            //
+            //
+            ext.clear();
+            final StringBuffer sb3 = new StringBuffer( "((a,b)ab,(((c,d)cd,e)cde,(f,(g,h)gh)fgh)cdefgh)abcdefgh" );
+            final Phylogeny t3 = factory.create( sb3, new NHXParser() )[ 0 ];
+            t3.getNode( "ab" ).setCollapse( true );
+            t3.getNode( "cd" ).setCollapse( true );
+            t3.getNode( "cde" ).setCollapse( true );
+            t3.getNode( "c" ).setCollapse( true );
+            t3.getNode( "d" ).setCollapse( true );
+            t3.getNode( "e" ).setCollapse( true );
+            t3.getNode( "gh" ).setCollapse( true );
+            t3.getNode( "fgh" ).setCollapse( true );
+            n = t3.getNode( "ab" );
+            ext = new ArrayList<PhylogenyNode>();
+            while ( n != null ) {
+                ext.add( n );
+                n = n.getNextExternalNodeWhileTakingIntoAccountCollapsedNodes();
+            }
+            if ( !ext.get( 0 ).getName().equals( "ab" ) ) {
+                return false;
+            }
+            if ( !ext.get( 1 ).getName().equals( "cde" ) ) {
+                return false;
+            }
+            if ( !ext.get( 2 ).getName().equals( "fgh" ) ) {
+                return false;
+            }
+            //
+            //
+            ext.clear();
+            final StringBuffer sb4 = new StringBuffer( "((a,b)ab,(((c,d)cd,e)cde,(f,(g,h)gh)fgh)cdefgh)abcdefgh" );
+            final Phylogeny t4 = factory.create( sb4, new NHXParser() )[ 0 ];
+            t4.getNode( "ab" ).setCollapse( true );
+            t4.getNode( "cd" ).setCollapse( true );
+            t4.getNode( "cde" ).setCollapse( true );
+            t4.getNode( "c" ).setCollapse( true );
+            t4.getNode( "d" ).setCollapse( true );
+            t4.getNode( "e" ).setCollapse( true );
+            t4.getNode( "gh" ).setCollapse( true );
+            t4.getNode( "fgh" ).setCollapse( true );
+            t4.getNode( "abcdefgh" ).setCollapse( true );
+            n = t4.getNode( "abcdefgh" );
+            if ( n.getNextExternalNodeWhileTakingIntoAccountCollapsedNodes() != null ) {
+                return false;
+            }
+            //
+            //
+            final StringBuffer sb5 = new StringBuffer( "((a,b)ab,(((c,d)cd,e)cde,(f,(g,h))fgh)cdefgh)abcdefgh" );
+            final Phylogeny t5 = factory.create( sb5, new NHXParser() )[ 0 ];
+            ext.clear();
+            n = t5.getFirstExternalNode();
+            while ( n != null ) {
+                ext.add( n );
+                n = n.getNextExternalNodeWhileTakingIntoAccountCollapsedNodes();
+            }
+            if ( ext.size() != 8 ) {
+                return false;
+            }
+            if ( !ext.get( 0 ).getName().equals( "a" ) ) {
+                return false;
+            }
+            if ( !ext.get( 1 ).getName().equals( "b" ) ) {
+                return false;
+            }
+            if ( !ext.get( 2 ).getName().equals( "c" ) ) {
+                return false;
+            }
+            if ( !ext.get( 3 ).getName().equals( "d" ) ) {
+                return false;
+            }
+            if ( !ext.get( 4 ).getName().equals( "e" ) ) {
+                return false;
+            }
+            if ( !ext.get( 5 ).getName().equals( "f" ) ) {
+                return false;
+            }
+            if ( !ext.get( 6 ).getName().equals( "g" ) ) {
+                return false;
+            }
+            if ( !ext.get( 7 ).getName().equals( "h" ) ) {
+                return false;
+            }
+            //
+            //
+            final StringBuffer sb6 = new StringBuffer( "((a,b)ab,(((c,d)cd,e)cde,(f,(g,h))fgh)cdefgh)abcdefgh" );
+            final Phylogeny t6 = factory.create( sb6, new NHXParser() )[ 0 ];
+            ext.clear();
+            t6.getNode( "ab" ).setCollapse( true );
+            n = t6.getNode( "ab" );
+            while ( n != null ) {
+                ext.add( n );
+                n = n.getNextExternalNodeWhileTakingIntoAccountCollapsedNodes();
+            }
+            if ( ext.size() != 7 ) {
+                return false;
+            }
+            if ( !ext.get( 0 ).getName().equals( "ab" ) ) {
+                return false;
+            }
+            if ( !ext.get( 1 ).getName().equals( "c" ) ) {
+                return false;
+            }
+            if ( !ext.get( 2 ).getName().equals( "d" ) ) {
+                return false;
+            }
+            if ( !ext.get( 3 ).getName().equals( "e" ) ) {
+                return false;
+            }
+            if ( !ext.get( 4 ).getName().equals( "f" ) ) {
+                return false;
+            }
+            if ( !ext.get( 5 ).getName().equals( "g" ) ) {
+                return false;
+            }
+            if ( !ext.get( 6 ).getName().equals( "h" ) ) {
+                return false;
+            }
+            //
+            //
+            final StringBuffer sb7 = new StringBuffer( "((a,b)ab,(((c,d)cd,e)cde,(f,(g,h))fgh)cdefgh)abcdefgh" );
+            final Phylogeny t7 = factory.create( sb7, new NHXParser() )[ 0 ];
+            ext.clear();
+            t7.getNode( "cd" ).setCollapse( true );
+            n = t7.getNode( "a" );
+            while ( n != null ) {
+                ext.add( n );
+                n = n.getNextExternalNodeWhileTakingIntoAccountCollapsedNodes();
+            }
+            if ( ext.size() != 7 ) {
+                return false;
+            }
+            if ( !ext.get( 0 ).getName().equals( "a" ) ) {
+                return false;
+            }
+            if ( !ext.get( 1 ).getName().equals( "b" ) ) {
+                return false;
+            }
+            if ( !ext.get( 2 ).getName().equals( "cd" ) ) {
+                return false;
+            }
+            if ( !ext.get( 3 ).getName().equals( "e" ) ) {
+                return false;
+            }
+            if ( !ext.get( 4 ).getName().equals( "f" ) ) {
+                return false;
+            }
+            if ( !ext.get( 5 ).getName().equals( "g" ) ) {
+                return false;
+            }
+            if ( !ext.get( 6 ).getName().equals( "h" ) ) {
+                return false;
+            }
+            //
+            //
+            final StringBuffer sb8 = new StringBuffer( "((a,b)ab,(((c,d)cd,e)cde,(f,(g,h))fgh)cdefgh)abcdefgh" );
+            final Phylogeny t8 = factory.create( sb8, new NHXParser() )[ 0 ];
+            ext.clear();
+            t8.getNode( "cd" ).setCollapse( true );
+            t8.getNode( "c" ).setCollapse( true );
+            t8.getNode( "d" ).setCollapse( true );
+            n = t8.getNode( "a" );
+            while ( n != null ) {
+                ext.add( n );
+                n = n.getNextExternalNodeWhileTakingIntoAccountCollapsedNodes();
+            }
+            if ( ext.size() != 7 ) {
+                return false;
+            }
+            if ( !ext.get( 0 ).getName().equals( "a" ) ) {
+                return false;
+            }
+            if ( !ext.get( 1 ).getName().equals( "b" ) ) {
+                return false;
+            }
+            if ( !ext.get( 2 ).getName().equals( "cd" ) ) {
+                System.out.println( "2 fail" );
+                return false;
+            }
+            if ( !ext.get( 3 ).getName().equals( "e" ) ) {
+                return false;
+            }
+            if ( !ext.get( 4 ).getName().equals( "f" ) ) {
+                return false;
+            }
+            if ( !ext.get( 5 ).getName().equals( "g" ) ) {
+                return false;
+            }
+            if ( !ext.get( 6 ).getName().equals( "h" ) ) {
+                return false;
+            }
+            //
+            //
+            final StringBuffer sb9 = new StringBuffer( "((a,b)ab,(((c,d)cd,e)cde,(f,(g,h)gh)fgh)cdefgh)abcdefgh" );
+            final Phylogeny t9 = factory.create( sb9, new NHXParser() )[ 0 ];
+            ext.clear();
+            t9.getNode( "gh" ).setCollapse( true );
+            n = t9.getNode( "a" );
+            while ( n != null ) {
+                ext.add( n );
+                n = n.getNextExternalNodeWhileTakingIntoAccountCollapsedNodes();
+            }
+            if ( ext.size() != 7 ) {
+                return false;
+            }
+            if ( !ext.get( 0 ).getName().equals( "a" ) ) {
+                return false;
+            }
+            if ( !ext.get( 1 ).getName().equals( "b" ) ) {
+                return false;
+            }
+            if ( !ext.get( 2 ).getName().equals( "c" ) ) {
+                return false;
+            }
+            if ( !ext.get( 3 ).getName().equals( "d" ) ) {
+                return false;
+            }
+            if ( !ext.get( 4 ).getName().equals( "e" ) ) {
+                return false;
+            }
+            if ( !ext.get( 5 ).getName().equals( "f" ) ) {
+                return false;
+            }
+            if ( !ext.get( 6 ).getName().equals( "gh" ) ) {
+                return false;
+            }
+            //
+            //
+            final StringBuffer sb10 = new StringBuffer( "((a,b)ab,(((c,d)cd,e)cde,(f,(g,h)gh)fgh)cdefgh)abcdefgh" );
+            final Phylogeny t10 = factory.create( sb10, new NHXParser() )[ 0 ];
+            ext.clear();
+            t10.getNode( "gh" ).setCollapse( true );
+            t10.getNode( "g" ).setCollapse( true );
+            t10.getNode( "h" ).setCollapse( true );
+            n = t10.getNode( "a" );
+            while ( n != null ) {
+                ext.add( n );
+                n = n.getNextExternalNodeWhileTakingIntoAccountCollapsedNodes();
+            }
+            if ( ext.size() != 7 ) {
+                return false;
+            }
+            if ( !ext.get( 0 ).getName().equals( "a" ) ) {
+                return false;
+            }
+            if ( !ext.get( 1 ).getName().equals( "b" ) ) {
+                return false;
+            }
+            if ( !ext.get( 2 ).getName().equals( "c" ) ) {
+                return false;
+            }
+            if ( !ext.get( 3 ).getName().equals( "d" ) ) {
+                return false;
+            }
+            if ( !ext.get( 4 ).getName().equals( "e" ) ) {
+                return false;
+            }
+            if ( !ext.get( 5 ).getName().equals( "f" ) ) {
+                return false;
+            }
+            if ( !ext.get( 6 ).getName().equals( "gh" ) ) {
+                return false;
+            }
+            //
+            //
+            final StringBuffer sb11 = new StringBuffer( "((a,b)ab,(((c,d)cd,e)cde,(f,(g,h)gh)fgh)cdefgh)abcdefgh" );
+            final Phylogeny t11 = factory.create( sb11, new NHXParser() )[ 0 ];
+            ext.clear();
+            t11.getNode( "gh" ).setCollapse( true );
+            t11.getNode( "fgh" ).setCollapse( true );
+            n = t11.getNode( "a" );
+            while ( n != null ) {
+                ext.add( n );
+                n = n.getNextExternalNodeWhileTakingIntoAccountCollapsedNodes();
+            }
+            if ( ext.size() != 6 ) {
+                return false;
+            }
+            if ( !ext.get( 0 ).getName().equals( "a" ) ) {
+                return false;
+            }
+            if ( !ext.get( 1 ).getName().equals( "b" ) ) {
+                return false;
+            }
+            if ( !ext.get( 2 ).getName().equals( "c" ) ) {
+                return false;
+            }
+            if ( !ext.get( 3 ).getName().equals( "d" ) ) {
+                return false;
+            }
+            if ( !ext.get( 4 ).getName().equals( "e" ) ) {
+                return false;
+            }
+            if ( !ext.get( 5 ).getName().equals( "fgh" ) ) {
+                return false;
+            }
+            //
+            //
+            final StringBuffer sb12 = new StringBuffer( "((a,b)ab,(((c,d)cd,e)cde,(f,(g,h)gh)fgh)cdefgh)abcdefgh" );
+            final Phylogeny t12 = factory.create( sb12, new NHXParser() )[ 0 ];
+            ext.clear();
+            t12.getNode( "gh" ).setCollapse( true );
+            t12.getNode( "fgh" ).setCollapse( true );
+            t12.getNode( "g" ).setCollapse( true );
+            t12.getNode( "h" ).setCollapse( true );
+            t12.getNode( "f" ).setCollapse( true );
+            n = t12.getNode( "a" );
+            while ( n != null ) {
+                ext.add( n );
+                n = n.getNextExternalNodeWhileTakingIntoAccountCollapsedNodes();
+            }
+            if ( ext.size() != 6 ) {
+                return false;
+            }
+            if ( !ext.get( 0 ).getName().equals( "a" ) ) {
+                return false;
+            }
+            if ( !ext.get( 1 ).getName().equals( "b" ) ) {
+                return false;
+            }
+            if ( !ext.get( 2 ).getName().equals( "c" ) ) {
+                return false;
+            }
+            if ( !ext.get( 3 ).getName().equals( "d" ) ) {
+                return false;
+            }
+            if ( !ext.get( 4 ).getName().equals( "e" ) ) {
+                return false;
+            }
+            if ( !ext.get( 5 ).getName().equals( "fgh" ) ) {
+                return false;
+            }
+            //
+            //
+            final StringBuffer sb13 = new StringBuffer( "((a,b)ab,(((c,d)cd,e)cde,(f,(g,h)gh)fgh)cdefgh)abcdefgh" );
+            final Phylogeny t13 = factory.create( sb13, new NHXParser() )[ 0 ];
+            ext.clear();
+            t13.getNode( "ab" ).setCollapse( true );
+            t13.getNode( "b" ).setCollapse( true );
+            t13.getNode( "fgh" ).setCollapse( true );
+            t13.getNode( "gh" ).setCollapse( true );
+            n = t13.getNode( "ab" );
+            while ( n != null ) {
+                ext.add( n );
+                n = n.getNextExternalNodeWhileTakingIntoAccountCollapsedNodes();
+            }
+            if ( ext.size() != 5 ) {
+                return false;
+            }
+            if ( !ext.get( 0 ).getName().equals( "ab" ) ) {
+                return false;
+            }
+            if ( !ext.get( 1 ).getName().equals( "c" ) ) {
+                return false;
+            }
+            if ( !ext.get( 2 ).getName().equals( "d" ) ) {
+                return false;
+            }
+            if ( !ext.get( 3 ).getName().equals( "e" ) ) {
+                return false;
+            }
+            if ( !ext.get( 4 ).getName().equals( "fgh" ) ) {
+                return false;
+            }
+            //
+            //
+            final StringBuffer sb14 = new StringBuffer( "((a,b,0)ab,(((c,d)cd,e)cde,(f,(g,h,1,2)gh,0)fgh)cdefgh)abcdefgh" );
+            final Phylogeny t14 = factory.create( sb14, new NHXParser() )[ 0 ];
+            ext.clear();
+            t14.getNode( "ab" ).setCollapse( true );
+            t14.getNode( "a" ).setCollapse( true );
+            t14.getNode( "fgh" ).setCollapse( true );
+            t14.getNode( "gh" ).setCollapse( true );
+            n = t14.getNode( "ab" );
+            while ( n != null ) {
+                ext.add( n );
+                n = n.getNextExternalNodeWhileTakingIntoAccountCollapsedNodes();
+            }
+            if ( ext.size() != 5 ) {
+                return false;
+            }
+            if ( !ext.get( 0 ).getName().equals( "ab" ) ) {
+                return false;
+            }
+            if ( !ext.get( 1 ).getName().equals( "c" ) ) {
+                return false;
+            }
+            if ( !ext.get( 2 ).getName().equals( "d" ) ) {
+                return false;
+            }
+            if ( !ext.get( 3 ).getName().equals( "e" ) ) {
+                return false;
+            }
+            if ( !ext.get( 4 ).getName().equals( "fgh" ) ) {
+                return false;
+            }
+            //
+            //
+            final StringBuffer sb15 = new StringBuffer( "((a,b,0)ab,(((c,d)cd,e)cde,x,(f,(g,h,1,2)gh,0)fgh)cdefgh)abcdefgh" );
+            final Phylogeny t15 = factory.create( sb15, new NHXParser() )[ 0 ];
+            ext.clear();
+            t15.getNode( "ab" ).setCollapse( true );
+            t15.getNode( "a" ).setCollapse( true );
+            t15.getNode( "fgh" ).setCollapse( true );
+            t15.getNode( "gh" ).setCollapse( true );
+            n = t15.getNode( "ab" );
+            while ( n != null ) {
+                ext.add( n );
+                n = n.getNextExternalNodeWhileTakingIntoAccountCollapsedNodes();
+            }
+            if ( ext.size() != 6 ) {
+                return false;
+            }
+            if ( !ext.get( 0 ).getName().equals( "ab" ) ) {
+                return false;
+            }
+            if ( !ext.get( 1 ).getName().equals( "c" ) ) {
+                return false;
+            }
+            if ( !ext.get( 2 ).getName().equals( "d" ) ) {
+                return false;
+            }
+            if ( !ext.get( 3 ).getName().equals( "e" ) ) {
+                return false;
+            }
+            if ( !ext.get( 4 ).getName().equals( "x" ) ) {
+                return false;
+            }
+            if ( !ext.get( 5 ).getName().equals( "fgh" ) ) {
+                return false;
+            }
+            //
+            //
+            final StringBuffer sb16 = new StringBuffer( "((a,b,0)ab,(((c,d)cd,e)cde,x,(f,(g,h,1,2)gh,0)fgh)cdefgh)abcdefgh" );
+            final Phylogeny t16 = factory.create( sb16, new NHXParser() )[ 0 ];
+            ext.clear();
+            t16.getNode( "ab" ).setCollapse( true );
+            t16.getNode( "a" ).setCollapse( true );
+            t16.getNode( "fgh" ).setCollapse( true );
+            t16.getNode( "gh" ).setCollapse( true );
+            t16.getNode( "cd" ).setCollapse( true );
+            t16.getNode( "cde" ).setCollapse( true );
+            t16.getNode( "d" ).setCollapse( true );
+            t16.getNode( "x" ).setCollapse( true );
+            n = t16.getNode( "ab" );
+            while ( n != null ) {
+                ext.add( n );
+                n = n.getNextExternalNodeWhileTakingIntoAccountCollapsedNodes();
+            }
+            if ( ext.size() != 4 ) {
+                return false;
+            }
+            if ( !ext.get( 0 ).getName().equals( "ab" ) ) {
+                return false;
+            }
+            if ( !ext.get( 1 ).getName().equals( "cde" ) ) {
+                return false;
+            }
+            if ( !ext.get( 2 ).getName().equals( "x" ) ) {
+                return false;
+            }
+            if ( !ext.get( 3 ).getName().equals( "fgh" ) ) {
+                return false;
+            }
+        }
+        catch ( final Exception e ) {
+            e.printStackTrace( System.out );
+            return false;
+        }
+        return true;
+    }
+
+    private static boolean testMsaQualityMethod() {
+        try {
+            final Sequence s0 = BasicSequence.createAaSequence( "a", "ABAXEFGHIJ" );
+            final Sequence s1 = BasicSequence.createAaSequence( "b", "ABBXEFGHIJ" );
+            final Sequence s2 = BasicSequence.createAaSequence( "c", "AXCXEFGHIJ" );
+            final Sequence s3 = BasicSequence.createAaSequence( "d", "AXDDEFGHIJ" );
+            final List<Sequence> l = new ArrayList<Sequence>();
+            l.add( s0 );
+            l.add( s1 );
+            l.add( s2 );
+            l.add( s3 );
+            final Msa msa = BasicMsa.createInstance( l );
+            if ( !isEqual( 1, MsaMethods.calculateIdentityRatio( msa, 0 ) ) ) {
+                return false;
+            }
+            if ( !isEqual( 0.5, MsaMethods.calculateIdentityRatio( msa, 1 ) ) ) {
+                return false;
+            }
+            if ( !isEqual( 0.25, MsaMethods.calculateIdentityRatio( msa, 2 ) ) ) {
+                return false;
+            }
+            if ( !isEqual( 0.75, MsaMethods.calculateIdentityRatio( msa, 3 ) ) ) {
+                return false;
+            }
+        }
+        catch ( final Exception e ) {
+            e.printStackTrace( System.out );
+            return false;
+        }
+        return true;
+    }
+
+    private static boolean testSequenceIdParsing() {
+        try {
+            Identifier id = SequenceIdParser.parse( "gb_ADF31344_segmented_worms_" );
+            if ( ( id == null ) || ForesterUtil.isEmpty( id.getValue() ) || ForesterUtil.isEmpty( id.getProvider() )
+                    || !id.getValue().equals( "ADF31344" ) || !id.getProvider().equals( "ncbi" ) ) {
+                if ( id != null ) {
+                    System.out.println( "value   =" + id.getValue() );
+                    System.out.println( "provider=" + id.getProvider() );
+                }
+                return false;
+            }
+            //
+            id = SequenceIdParser.parse( "segmented worms|gb_ADF31344" );
+            if ( ( id == null ) || ForesterUtil.isEmpty( id.getValue() ) || ForesterUtil.isEmpty( id.getProvider() )
+                    || !id.getValue().equals( "ADF31344" ) || !id.getProvider().equals( "ncbi" ) ) {
+                if ( id != null ) {
+                    System.out.println( "value   =" + id.getValue() );
+                    System.out.println( "provider=" + id.getProvider() );
+                }
+                return false;
+            }
+            //
+            id = SequenceIdParser.parse( "segmented worms gb_ADF31344 and more" );
+            if ( ( id == null ) || ForesterUtil.isEmpty( id.getValue() ) || ForesterUtil.isEmpty( id.getProvider() )
+                    || !id.getValue().equals( "ADF31344" ) || !id.getProvider().equals( "ncbi" ) ) {
+                if ( id != null ) {
+                    System.out.println( "value   =" + id.getValue() );
+                    System.out.println( "provider=" + id.getProvider() );
+                }
+                return false;
+            }
+            // 
+            id = SequenceIdParser.parse( "gb_AAA96518_1" );
+            if ( ( id == null ) || ForesterUtil.isEmpty( id.getValue() ) || ForesterUtil.isEmpty( id.getProvider() )
+                    || !id.getValue().equals( "AAA96518" ) || !id.getProvider().equals( "ncbi" ) ) {
+                if ( id != null ) {
+                    System.out.println( "value   =" + id.getValue() );
+                    System.out.println( "provider=" + id.getProvider() );
+                }
+                return false;
+            }
+            // 
+            id = SequenceIdParser.parse( "gb_EHB07727_1_rodents_" );
+            if ( ( id == null ) || ForesterUtil.isEmpty( id.getValue() ) || ForesterUtil.isEmpty( id.getProvider() )
+                    || !id.getValue().equals( "EHB07727" ) || !id.getProvider().equals( "ncbi" ) ) {
+                if ( id != null ) {
+                    System.out.println( "value   =" + id.getValue() );
+                    System.out.println( "provider=" + id.getProvider() );
+                }
+                return false;
+            }
+            // 
+            id = SequenceIdParser.parse( "dbj_BAF37827_1_turtles_" );
+            if ( ( id == null ) || ForesterUtil.isEmpty( id.getValue() ) || ForesterUtil.isEmpty( id.getProvider() )
+                    || !id.getValue().equals( "BAF37827" ) || !id.getProvider().equals( "ncbi" ) ) {
+                if ( id != null ) {
+                    System.out.println( "value   =" + id.getValue() );
+                    System.out.println( "provider=" + id.getProvider() );
+                }
+                return false;
+            }
+            // 
+            id = SequenceIdParser.parse( "emb_CAA73223_1_primates_" );
+            if ( ( id == null ) || ForesterUtil.isEmpty( id.getValue() ) || ForesterUtil.isEmpty( id.getProvider() )
+                    || !id.getValue().equals( "CAA73223" ) || !id.getProvider().equals( "ncbi" ) ) {
+                if ( id != null ) {
+                    System.out.println( "value   =" + id.getValue() );
+                    System.out.println( "provider=" + id.getProvider() );
+                }
+                return false;
+            }
+            // 
+            id = SequenceIdParser.parse( "mites|ref_XP_002434188_1" );
+            if ( ( id == null ) || ForesterUtil.isEmpty( id.getValue() ) || ForesterUtil.isEmpty( id.getProvider() )
+                    || !id.getValue().equals( "XP_002434188" ) || !id.getProvider().equals( "refseq" ) ) {
+                if ( id != null ) {
+                    System.out.println( "value   =" + id.getValue() );
+                    System.out.println( "provider=" + id.getProvider() );
+                }
+                return false;
+            }
+            // 
+            id = SequenceIdParser.parse( "mites_ref_XP_002434188_1_bla_XP_12345" );
+            if ( ( id == null ) || ForesterUtil.isEmpty( id.getValue() ) || ForesterUtil.isEmpty( id.getProvider() )
+                    || !id.getValue().equals( "XP_002434188" ) || !id.getProvider().equals( "refseq" ) ) {
+                if ( id != null ) {
+                    System.out.println( "value   =" + id.getValue() );
+                    System.out.println( "provider=" + id.getProvider() );
+                }
+                return false;
+            }
+            // 
+            id = SequenceIdParser.parse( "XP_12345" );
+            if ( id != null ) {
                 return false;
             }
+            // lcl_91970_unknown_
         }
         catch ( final Exception e ) {
             e.printStackTrace( System.out );