+ public void write( final Writer w, final MSA_FORMAT format ) throws IOException {
+ switch ( format ) {
+ case PHYLIP:
+ writeToPhylip( w );
+ break;
+ case FASTA:
+ writeToFasta( w );
+ break;
+ case NEXUS:
+ writeToNexus( w );
+ break;
+ default:
+ throw new RuntimeException( "unknown format " + format );
+ }
+ }
+
+ private short determineMaxIdLength() {
+ short max = 0;
+ for( int row = 0; row < getNumberOfSequences(); ++row ) {
+ final short l = ( short ) getIdentifier( row ).length();
+ if ( l > max ) {
+ max = l;
+ }
+ }
+ return max;
+ }
+
+ private void writeToFasta( final Writer w ) throws IOException {
+ SequenceWriter.writeSeqs( asSequenceList(), w, SEQ_FORMAT.FASTA, 100 );
+ }
+
+ private void writeToNexus( final Writer w ) throws IOException {
+ final int max = determineMaxIdLength() + 1;
+ w.write( "Begin Data;" );
+ w.write( ForesterUtil.LINE_SEPARATOR );
+ w.write( " Dimensions NTax=" + getNumberOfSequences() );
+ w.write( " NChar=" + getLength() );
+ w.write( ";" );
+ w.write( ForesterUtil.LINE_SEPARATOR );
+ w.write( " Format DataType=Protein Interleave=No gap=-;" );
+ w.write( ForesterUtil.LINE_SEPARATOR );
+ w.write( " Matrix" );
+ w.write( ForesterUtil.LINE_SEPARATOR );
+ for( int row = 0; row < getNumberOfSequences(); ++row ) {
+ final MolecularSequence seq = getSequence( row );
+ final String s = seq.getMolecularSequenceAsString();
+ w.write( " " );
+ w.write( ForesterUtil.pad( getIdentifier( row ).replace( ' ', '_' ), max, ' ', false ).toString() );
+ w.write( " " );
+ w.write( s );
+ w.write( ForesterUtil.LINE_SEPARATOR );
+ }
+ w.write( " ;" );
+ w.write( ForesterUtil.LINE_SEPARATOR );
+ w.write( "End;" );
+ w.write( ForesterUtil.LINE_SEPARATOR );
+ }
+
+ private void writeToPhylip( final Writer w ) throws IOException {