X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=forester%2Fjava%2Fsrc%2Forg%2Fforester%2Fmsa%2FMafft.java;h=ec72c0d01e5a068a7951ab6dd65f46f877976483;hb=ea402c168ecf73c1824d023a91331494939f2641;hp=644b5b52a469b00c5e1dce56c6ef1b297127594c;hpb=48f7a89be9d34f1930a1f863e608235cc27184c5;p=jalview.git diff --git a/forester/java/src/org/forester/msa/Mafft.java b/forester/java/src/org/forester/msa/Mafft.java index 644b5b5..ec72c0d 100644 --- a/forester/java/src/org/forester/msa/Mafft.java +++ b/forester/java/src/org/forester/msa/Mafft.java @@ -5,7 +5,7 @@ // Copyright (C) 2010 Christian M Zmasek // Copyright (C) 2010 Sanford-Burnham Medical Research Institute // All rights reserved -// +// // This library is free software; you can redistribute it and/or // modify it under the terms of the GNU Lesser General Public // License as published by the Free Software Foundation; either @@ -15,7 +15,7 @@ // but WITHOUT ANY WARRANTY; without even the implied warranty of // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU // Lesser General Public License for more details. -// +// // You should have received a copy of the GNU Lesser General Public // License along with this library; if not, write to the Free Software // Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA @@ -25,53 +25,33 @@ package org.forester.msa; +import java.io.BufferedWriter; import java.io.File; +import java.io.FileWriter; import java.io.IOException; import java.util.ArrayList; import java.util.List; import org.forester.io.parsers.FastaParser; -import org.forester.util.ForesterUtil; +import org.forester.io.writers.SequenceWriter; +import org.forester.io.writers.SequenceWriter.SEQ_FORMAT; +import org.forester.sequence.Sequence; import org.forester.util.SystemCommandExecutor; -public final class Mafft implements MsaInferrer { +public final class Mafft extends MsaInferrer { private final static String DEFAULT_PARAMETERS = "--maxiterate 1000 --localpair"; private String _error; private int _exit_code; private final String _path_to_prg; - public static MsaInferrer createInstance( final String path_to_prg ) { + public static MsaInferrer createInstance( final String path_to_prg ) throws IOException { return new Mafft( path_to_prg ); } - private static String getPathToCmd() { - //TODO this needs to come from env variable, etc. - 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/SOFTWARE/MSA/MAFFT/mafft-6.832-without-extensions/scripts/mafft"; - } - return path; - } - - public static boolean isInstalled() { - return SystemCommandExecutor.isExecuteableFile( new File( getPathToCmd() ) ); - } - - public static MsaInferrer createInstance() { - return createInstance( getPathToCmd() ); - } - - private Mafft( final String path_to_prg ) { - if ( !SystemCommandExecutor.isExecuteableFile( new File( path_to_prg ) ) ) { - throw new IllegalArgumentException( "cannot execute MAFFT via [" + path_to_prg + "]" ); + private Mafft( final String path_to_prg ) throws IOException { + if ( !isInstalled( path_to_prg ) ) { + throw new IOException( "cannot execute MAFFT with \"" + path_to_prg + "\"" ); } _path_to_prg = new String( path_to_prg ); init(); @@ -82,18 +62,28 @@ public final class Mafft implements MsaInferrer { } @Override - public Object clone() { - throw new NoSuchMethodError(); - } - public String getErrorDescription() { return _error; } + @Override public int getExitCode() { return _exit_code; } + @Override + public Msa infer( final List seqs, final List opts ) throws IOException, InterruptedException { + final File file = File.createTempFile( "__mafft_input_", ".fasta" ); + file.deleteOnExit(); + final BufferedWriter writer = new BufferedWriter( new FileWriter( file ) ); + SequenceWriter.writeSeqs( seqs, writer, SEQ_FORMAT.FASTA, 100 ); + writer.close(); + final Msa msa = infer( file, opts ); + file.delete(); + return msa; + } + + @Override public Msa infer( final File path_to_input_seqs, final List opts ) throws IOException, InterruptedException { init(); final List my_opts = new ArrayList(); @@ -102,17 +92,19 @@ public final class Mafft implements MsaInferrer { my_opts.add( opts.get( i ) ); } my_opts.add( path_to_input_seqs.getAbsolutePath() ); - final SystemCommandExecutor commandExecutor = new SystemCommandExecutor( my_opts ); - final int _exit_code = commandExecutor.executeCommand(); + final SystemCommandExecutor command_executor = new SystemCommandExecutor( my_opts ); + final int _exit_code = command_executor.executeCommand(); + final StringBuilder stderr = command_executor.getStandardErrorFromCommand(); + _error = stderr.toString(); if ( _exit_code != 0 ) { - throw new IOException( "MAFFT failed, exit code: " + _exit_code ); + throw new IOException( "MAFFT program failed, exit code: " + _exit_code + "\nCommand:\n" + my_opts + + "\nError:\n" + stderr ); + } + final StringBuilder stdout = command_executor.getStandardOutputFromCommand(); + if ( ( stdout == null ) || ( stdout.length() < 2 ) ) { + throw new IOException( "MAFFT program did not produce any output\nCommand:\n" + my_opts + "\nError:\n" + + stderr ); } - final StringBuilder stdout = commandExecutor.getStandardOutputFromCommand(); - final StringBuilder stderr = commandExecutor.getStandardErrorFromCommand(); - System.out.println( stdout ); - System.out.println(); - System.out.println( stderr ); - _error = stderr.toString(); final Msa msa = FastaParser.parseMsa( stdout.toString() ); return msa; }