Compare two genomes for similarity
Created: 2010-03-17 11:01:09
Last updated: 2010-03-17 11:05:00
The workflow takes as input two nucleotide/protein sequences in FastA format, requests a minimum percentage similarity threshold and then checks this against the similarity of the two sequences as determined by the Needleman-Wunsch global alignment.
Preview
Run
Run this Workflow in the Taverna Workbench...
Workflow Components
Authors (1)
Gregg Iceton
Newcastle University |
Titles (2)
Nucleotide sequence similarity checker |
Nucleotide or protein sequence similarity checker |
Descriptions (1)
The workflow takes as input two nucleotide/protein sequences in FastA format, requests a minimum percentage similarity threshold and then checks this against the similarity of the two sequences as determined by the Needleman-Wunsch global alignment. |
Dependencies (0)
Inputs (2)
Name |
Description |
A_Sequence |
Input the first sequence in FastA format
|
B_Sequence |
Input the second sequence in FastA format
|
Processors (11)
Name |
Type |
Description |
needle |
soaplab |
Endpointhttp://www.ebi.ac.uk/soaplab/services/alignment_global.needle |
Regex1 |
stringconstant |
Value\n |
Byte_Array_To_String |
localworker |
Scriptstring = new String(bytes);
|
Split_string_into_string_list_by_regular_expression |
localworker |
ScriptList split = new ArrayList();
if (!string.equals("")) {
String regexString = ",";
if (regex != void) {
regexString = regex;
}
String[] result = string.split(regexString);
for (int i = 0; i < result.length; i++) {
split.add(result[i]);
}
}
|
Regex2 |
stringconstant |
Value.*Similarity.+\((....) |
Filter_list_of_strings_extracting_match_to_a_regex |
localworker |
Scriptimport java.util.regex.*;
filteredlist = new ArrayList();
Pattern thePat = Pattern.compile(regex);
int theGroup = Integer.parseInt(group);
for (Iterator i = stringlist.iterator(); i.hasNext();) {
String item = (String) i.next();
Matcher matcher = thePat.matcher(item);
if (matcher.find()) {
filteredlist.add(matcher.group(theGroup));
}
}
|
Group |
stringconstant |
Value1 |
checkSimilarity |
beanshell |
Scriptimport java.lang.Double;
double seqSimDbl = Double.valueOf(seqSim.trim()).doubleValue();
double minSimDbl = Double.valueOf(minSim.trim()).doubleValue();
int i = Double.compare(seqSimDbl,minSimDbl);
if (i >= 0)
{
simBln = 1;
} else
{
simBln = 0;
} |
Ask |
localworker |
Scriptimport javax.swing.JOptionPane;
answer = JOptionPane.showInputDialog(null, (message == void ? null : message), (title == void ? null : title), JOptionPane.QUESTION_MESSAGE);
|
message_value |
stringconstant |
ValuePlease enter the minimum percentage sequence similarity required |
title_value |
stringconstant |
ValueSequence similarity |
Beanshells (1)
Name |
Description |
Inputs |
Outputs |
checkSimilarity |
|
seqSim
minSim
|
simBln
seqSimDbl
minSimDbl
|
Outputs (1)
Name |
Description |
Carry_on |
Boolean output (0 or 1) - 0 indicates the sequence similarity was lower than the input threshold, 1 indicates the similarity was equal to or higher than the threshold.
|
Datalinks (13)
Source |
Sink |
A_Sequence |
needle:asequence_direct_data |
B_Sequence |
needle:bsequence_direct_data |
needle:outfile |
Byte_Array_To_String:bytes |
Regex1:value |
Split_string_into_string_list_by_regular_expression:regex |
Byte_Array_To_String:string |
Split_string_into_string_list_by_regular_expression:string |
Regex2:value |
Filter_list_of_strings_extracting_match_to_a_regex:regex |
Split_string_into_string_list_by_regular_expression:split |
Filter_list_of_strings_extracting_match_to_a_regex:stringlist |
Group:value |
Filter_list_of_strings_extracting_match_to_a_regex:group |
Filter_list_of_strings_extracting_match_to_a_regex:filteredlist |
checkSimilarity:seqSim |
Ask:answer |
checkSimilarity:minSim |
message_value:value |
Ask:message |
title_value:value |
Ask:title |
checkSimilarity:simBln |
Carry_on |
Uploader
License
All versions of this Workflow are
licensed under:
Version 1 (earliest)
(of 4)
Credits (1)
(People/Groups)
Attributions (0)
(Workflows/Files)
None
Shared with Groups (1)
Featured In Packs (0)
None
Log in to add to one of your Packs
Attributed By (0)
(Workflows/Files)
None
Favourited By (0)
No one
Statistics
Other workflows that use similar services
(0)
There are no workflows in myExperiment that use similar services to this Workflow.
Comments (0)
No comments yet
Log in to make a comment