import jalview.ws.params.WsParamSetI;
import jalview.ws.params.simple.BooleanOption;
import jalview.ws.params.simple.DoubleParameter;
+import jalview.ws.params.simple.FileParameter;
import jalview.ws.params.simple.IntegerParameter;
import jalview.ws.params.simple.LogarithmicParameter;
import jalview.ws.params.simple.Option;
+import jalview.ws.params.simple.RadioChoiceParameter;
import jalview.ws.params.simple.StringParameter;
+import java.io.File;
import java.io.IOException;
import java.util.ArrayList;
import java.util.Arrays;
-import java.util.Collection;
import java.util.List;
import java.util.Scanner;
private static final String HMMSEARCH = "hmmsearch";
+ private static final String JACKHMMER = "jackhmmer";
+
private String name;
private List<WsParamSetI> presets = new ArrayList<>();
return new HMMERParamStore(HMMSEARCH, viewport);
}
+ public static HMMERParamStore forJackhmmer(AlignmentViewport viewport)
+ {
+ return new HMMERParamStore(JACKHMMER, viewport);
+ }
+
@Override
public List<WsParamSetI> getPresets()
{
case HMMBUILD:
getHMMBuildParams(args);
break;
+ case JACKHMMER:
+ getJackhmmerParams(args);
default:
}
/*
* 'Options'
*/
- addChoiceOfHmm(args);
-
- String names = Cache.getProperty(Preferences.HMMSEARCH_DBS);
- if (names != null && !names.isEmpty())
- {
- List<String> databases = new ArrayList<>();
- databases.add(MessageManager.getString(HMMSearch.THIS_ALIGNMENT_KEY));
- Scanner nameScanner = new Scanner(names);
-
- if (nameScanner.hasNext())
- {
- while (nameScanner.hasNext())
- {
- String next = nameScanner.next();
- if ("null".equals(next))
- {
- Cache.setProperty(Preferences.HMMSEARCH_DBS, "");
- Cache.setProperty(Preferences.HMMSEARCH_DB_PATHS, "");
- }
- else
- {
- databases.add(next);
- }
- }
- }
- nameScanner.close();
- args.add(new StringParameter(
- MessageManager.getString(HMMSearch.DATABASE_KEY),
- MessageManager.getString("label.database_for_hmmsearch"),
- true, MessageManager.getString(HMMSearch.THIS_ALIGNMENT_KEY),
- MessageManager.getString(HMMSearch.THIS_ALIGNMENT_KEY),
- databases));
- }
args.add(new BooleanOption(
MessageManager.getString(HMMSearch.AUTO_ALIGN_SEQS_KEY),
MessageManager.getString("label.auto_align_seqs_desc"), false,
- false, true, null));
+ false, false, null));
args.add(new BooleanOption(
MessageManager.getString(HMMSearch.USE_ACCESSIONS_KEY),
MessageManager.getString("label.use_accessions_desc"), false,
MessageManager.getString(HMMSearch.TRIM_TERMINI_KEY),
MessageManager.getString("label.trim_termini_desc"), false,
false, true, null));
+ args.add(new BooleanOption(
+ MessageManager.getString(HMMSearch.RETURN_N_NEW_SEQ),
+ MessageManager.getString("label.check_for_new_sequences_desc"),
+ false, false, false, null));
/*
* 'Parameters'
*/
+ addChoiceOfHmm(args);
+
+ // addChoiceOfDatabase(args);
+
+ String thisAlignment = MessageManager
+ .getString(HMMSearch.THIS_ALIGNMENT_KEY);
+ String database = MessageManager.getString("label.database");
+ args.add(new FileParameter(database, "", false, "", ""));
args.add(new IntegerParameter(
MessageManager.getString(HMMSearch.NUMBER_OF_RESULTS_KEY),
MessageManager.getString("label.number_of_results_desc"), true,
100, 0, 100000));
- args.add(new StringParameter(
+ args.add(new RadioChoiceParameter(
MessageManager.getString(HMMSearch.REPORTING_CUTOFF_KEY), null,
- true, HMMSearch.CUTOFF_NONE, HMMSearch.CUTOFF_NONE,
Arrays.asList(HMMSearch.CUTOFF_NONE, HMMSearch.CUTOFF_EVALUE,
- HMMSearch.CUTOFF_SCORE)));
+ HMMSearch.CUTOFF_SCORE),
+ HMMSearch.CUTOFF_EVALUE));
args.add(new LogarithmicParameter(
MessageManager.getString(HMMSearch.SEQ_EVALUE_KEY),
MessageManager.getString("label.seq_e_value_desc"), false, 1D,
- 1E-38, 10D));
+ 1E-100, 10D));
args.add(new LogarithmicParameter(
MessageManager.getString(HMMSearch.DOM_EVALUE_KEY),
MessageManager.getString("label.dom_e_value_desc"), false, 1D,
- 1E-38, 10D));
+ 1E-100, 10D));
args.add(
new DoubleParameter(
MessageManager.getString(HMMSearch.SEQ_SCORE_KEY),
}
/**
+ * Answers default parameters for jackhmmer, taking into account any configured
+ * as user preferences
+ *
+ * @param args
+ */
+ private void getJackhmmerParams(List<ArgumentI> args)
+ {
+
+ /*
+ * 'Parameters'
+ */
+ addChoiceOfSequence(args);
+
+ // addChoiceOfDatabase(args);
+
+ String database = MessageManager.getString("label.database");
+ args.add(new FileParameter(database, "", false, "", ""));
+ args.add(new RadioChoiceParameter(
+ MessageManager.getString(JackHMMER.REPORTING_CUTOFF_KEY), null,
+ Arrays.asList(MessageManager.getString(JackHMMER.CUTOFF_NONE),
+ MessageManager.getString(JackHMMER.CUTOFF_EVALUE),
+ MessageManager.getString(JackHMMER.CUTOFF_SCORE)),
+ MessageManager.getString(JackHMMER.CUTOFF_EVALUE)));
+ args.add(new LogarithmicParameter(
+ MessageManager.getString(JackHMMER.SEQ_EVALUE_KEY),
+ MessageManager.getString("label.seq_e_value_desc"), false, 1D,
+ 1E-38, 10D));
+ args.add(new LogarithmicParameter(
+ MessageManager.getString(JackHMMER.DOM_EVALUE_KEY),
+ MessageManager.getString("label.dom_e_value_desc"), false, 1D,
+ 1E-38, 10D));
+ args.add(new DoubleParameter(
+ MessageManager.getString(JackHMMER.SEQ_SCORE_KEY),
+ MessageManager.getString("label.seq_score_desc"), false, 0d, 0d,
+ 1000d));
+ args.add(new DoubleParameter(
+ MessageManager.getString(JackHMMER.DOM_SCORE_KEY),
+ MessageManager.getString("label.dom_score_desc"), false, 0d, 0d,
+ 1000d));
+ }
+
+ /**
+ * Constructs a choice parameter for database to search; always includes 'this
+ * alignment', and also includes any databases held under user preferences key
+ * "HMMSEARCH_DBS" as a comma-delimited list
+ *
+ * @param args
+ */
+ protected void addChoiceOfDatabase(List<ArgumentI> args)
+ {
+ String names = Cache.getProperty(Preferences.HMMSEARCH_DBS);
+ if (names == null || names.isEmpty())
+ {
+ return;
+ }
+
+ List<String> filePaths = new ArrayList<>();
+ List<String> fileNames = new ArrayList<>();
+
+ String thisAlignment = MessageManager.getString(HMMSearch.THIS_ALIGNMENT_KEY);
+ filePaths.add(thisAlignment);
+ fileNames.add(thisAlignment);
+
+ Scanner nameScanner = new Scanner(names);
+ nameScanner.useDelimiter(Preferences.COMMA);
+
+ while (nameScanner.hasNext())
+ {
+ String next = nameScanner.next();
+ if ("null".equals(next))
+ {
+ Cache.setProperty(Preferences.HMMSEARCH_DBS, "");
+ }
+ else
+ {
+ filePaths.add(next);
+ int pos = next.lastIndexOf(File.separator);
+ String fileName = next.substring(pos + 1);
+ fileNames.add(fileName);
+ }
+ }
+ nameScanner.close();
+ ArgumentI databasesOption = new StringParameter(
+ MessageManager.getString(HMMSearch.DATABASE_KEY),
+ MessageManager.getString("label.database_for_hmmsearch"), true,
+ thisAlignment,
+ thisAlignment,
+ filePaths, fileNames);
+ args.add(databasesOption);
+ }
+
+ /**
* Answers default parameters for hmmalign, taking into account any configured
* as user preferences
*
/**
* Adds an argument representing the choice of HMM sequences (profiles)
- * against which to perform align or search, provided more than one is found
+ * against which to perform align or search, provided at least one is found
*
* @param args
*/
protected void addChoiceOfHmm(List<ArgumentI> args)
{
List<SequenceI> hmms = viewport.getAlignment().getHmmSequences();
- if (hmms.size() > 1)
+ if (!hmms.isEmpty())
{
List<String> options = new ArrayList<>();
for (SequenceI hmmSeq : hmms)
String defseq = options.get(0);
ArgumentI arg = new StringParameter(
MessageManager.getString("label.use_hmm"), null, true, defseq,
- defseq, options);
+ defseq, options, null);
args.add(arg);
}
}
/**
+ * Adds an argument representing the choice of sequence against which to perform
+ * jackhmmer
+ *
+ * @param args
+ */
+ protected void addChoiceOfSequence(List<ArgumentI> args)
+ {
+ List<SequenceI> sequences = viewport.getAlignment().getSequences();
+
+ List<String> options = new ArrayList<>();
+
+ for (SequenceI seq : sequences)
+ {
+ options.add(seq.getName());
+ }
+
+ String defseq = options.get(0);
+ ArgumentI arg = new StringParameter(
+ MessageManager.getString("label.use_sequence"), null, true,
+ defseq,
+ defseq, options, null);
+ args.add(arg);
+ }
+
+ /**
* Answers default parameters for hmmbuild, taking into account any configured
* as user preferences
*
*/
if (!viewport.getAlignment().getGroups().isEmpty())
{
- Collection<String> options = new ArrayList<>();
+ List<String> options = new ArrayList<>();
options.add(MessageManager.getString("label.alignment"));
options.add(MessageManager.getString("label.groups_and_alignment"));
options.add(MessageManager.getString("label.groups"));