BiomartAndEMBOSSAnalysis
Created: 2009-12-15 22:12:00
Using Biomart and EMBOSS soaplab services, This workflow retrieves a number of sequences from 3 species: mouse, human, rat; align them, and returns a plot of the alignment result. Corresponding sequence ids are also returned.
Preview
Run
Run this Workflow in the Taverna Workbench...
Workflow Components
Authors (1)
Titles (1)
Descriptions (1)
Using Biomart and EMBOSS soaplab services, This workflow retrieves a number of sequences from 3 species: mouse, human, rat; align them, and returns a plot of the alignment result. Corresponding sequence ids are also returned. |
Dependencies (0)
Processors (10)
Name |
Type |
Description |
GetUniqueHomolog |
beanshell |
ScriptList HSOut = new ArrayList();
List RatOut = new ArrayList();
List MouseOut = new ArrayList();
Map hsToMouse = new HashMap();
Iterator j = MouseGeneIDs.iterator();
for (Iterator i = HSGeneIDs.iterator(); i.hasNext();) {
String id = (String)i.next();
hsToMouse.put(id, j.next());
}
Map hsToRat = new HashMap();
j = RatGeneIDs.iterator();
for (Iterator i = HSGeneIDs.iterator(); i.hasNext();) {
String id = (String)i.next();
hsToRat.put(id, j.next());
}
// Build the unique outputs
for (Iterator i = hsToRat.keySet().iterator(); i.hasNext();) {
String hsID = (String)i.next();
String ratID = (String)hsToRat.get(hsID);
// Remove version number
// ratID = (ratID.split("."))[0];
String mouseID = (String)hsToMouse.get(hsID);
// Remove version number
//mouseId = (mouseID.split("."))[0];
if (ratID != null && mouseID != null && ratID.equals("")==false && mouseID.equals("")==false) {
HSOut.add(hsID);
RatOut.add(ratID.split("\\.")[0]);
MouseOut.add(mouseID.split("\\.")[0]);
}
} |
hsapiensGeneEnsembl |
biomart |
|
getHSapSequence |
biomart |
|
CreateFasta |
beanshell |
Scriptfasta = ">Human\n"+hsSeq+"\n>Mouse\n"+mmSeq+"\n>Rat\n"+rnSeq; |
getRNorSequence |
biomart |
|
getMMusSequence |
biomart |
|
FlattenImageList |
localworker |
Scriptflatten(inputs, outputs, depth) {
for (i = inputs.iterator(); i.hasNext();) {
element = i.next();
if (element instanceof Collection && depth > 0) {
flatten(element, outputs, depth - 1);
} else {
outputs.add(element);
}
}
}
outputlist = new ArrayList();
flatten(inputlist, outputlist, 1); |
seqret |
soaplab |
Endpointhttp://www.ebi.ac.uk/soaplab/emboss4/services/edit.seqret |
plot |
soaplab |
Endpointhttp://www.ebi.ac.uk/soaplab/emboss4/services/alignment_multiple.prettyplot |
emma |
soaplab |
Endpointhttp://www.ebi.ac.uk/soaplab/emboss4/services/alignment_multiple.emma |
Beanshells (2)
Name |
Description |
Inputs |
Outputs |
GetUniqueHomolog |
|
HSGeneIDs
MouseGeneIDs
RatGeneIDs
|
HSOut
RatOut
MouseOut
|
CreateFasta |
|
hsSeq
mmSeq
rnSeq
|
fasta
|
Outputs (4)
Name |
Description |
outputPlot |
|
HSapIDs |
|
MMusIDs |
|
RNorIDs |
|
Datalinks (17)
Source |
Sink |
hsapiensGeneEnsembl:hsapiens_gene_ensembl.ensembl_gene_id |
GetUniqueHomolog:HSGeneIDs |
hsapiensGeneEnsembl:hsapiens_gene_ensembl.rat_ensembl_gene |
GetUniqueHomolog:RatGeneIDs |
hsapiensGeneEnsembl:hsapiens_gene_ensembl.mouse_ensembl_gene |
GetUniqueHomolog:MouseGeneIDs |
GetUniqueHomolog:HSOut |
getHSapSequence:hsapiens_gene_ensembl.ensembl_gene_id_filter |
getHSapSequence:hsapiens_gene_ensembl.coding_gene_flank |
CreateFasta:hsSeq |
getRNorSequence:rnorvegicus_gene_ensembl.coding_gene_flank |
CreateFasta:rnSeq |
getMMusSequence:mmusculus_gene_ensembl.coding_gene_flank |
CreateFasta:mmSeq |
GetUniqueHomolog:RatOut |
getRNorSequence:rnorvegicus_gene_ensembl.ensembl_gene_id_filter |
GetUniqueHomolog:MouseOut |
getMMusSequence:mmusculus_gene_ensembl.ensembl_gene_id_filter |
plot:Graphics_in_PNG |
FlattenImageList:inputlist |
CreateFasta:fasta |
seqret:sequence_direct_data |
emma:outseq |
plot:sequences_direct_data |
seqret:outseq |
emma:sequence_direct_data |
FlattenImageList:outputlist |
outputPlot |
GetUniqueHomolog:HSOut |
HSapIDs |
GetUniqueHomolog:MouseOut |
MMusIDs |
GetUniqueHomolog:RatOut |
RNorIDs |
Uploader
License
All versions of this Workflow are
licensed under:
Version 1
(of 1)
Credits (1)
(People/Groups)
Attributions (0)
(Workflows/Files)
None
Shared with Groups (1)
Featured In Packs (1)
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