<?xml version="1.0" encoding="UTF-8"?>
<s:scufl xmlns:s="http://org.embl.ebi.escience/xscufl/0.1alpha" version="0.2" log="0">
  <s:workflowdescription lsid="urn:lsid:www.mygrid.org.uk:operation:Z6160QD7NK59" author="Tom Oinn, Katy Wolstencroft" title="Find homologues, discover surrounding SNPs and perform multiple alignment">Finds homologues for the given human gene (Ensembl gene ID) in R. Norvegicus and M. Musculus, search for SNPs in the region 1000bp upstream and downstream from all genes (including original search and homologues) producing a SNP report and in addition perform a multiple alignment on all homologues against the original gene sequence.

Human gene id is hardcoded to ENSG00000086506</s:workflowdescription>
  <s:processor name="Add1000">
    <s:beanshell>
      <s:scriptvalue>i = Integer.parseInt(input);
output = ""+(i+1000);</s:scriptvalue>
      <s:beanshellinputlist>
        <s:beanshellinput s:syntactictype="'text/plain'">input</s:beanshellinput>
      </s:beanshellinputlist>
      <s:beanshelloutputlist>
        <s:beanshelloutput s:syntactictype="'text/plain'">output</s:beanshelloutput>
      </s:beanshelloutputlist>
    </s:beanshell>
  </s:processor>
  <s:processor name="Subtract1000">
    <s:beanshell>
      <s:scriptvalue>i = Integer.parseInt(input);
output = ""+(i - 1000);</s:scriptvalue>
      <s:beanshellinputlist>
        <s:beanshellinput s:syntactictype="'text/plain'">input</s:beanshellinput>
      </s:beanshellinputlist>
      <s:beanshelloutputlist>
        <s:beanshelloutput s:syntactictype="'text/plain'">output</s:beanshelloutput>
      </s:beanshelloutputlist>
    </s:beanshell>
  </s:processor>
  <s:processor name="mouseReport">
    <s:beanshell>
      <s:scriptvalue>result = strand + "," + chromStart + "," + chrom + "," + snpID;</s:scriptvalue>
      <s:beanshellinputlist>
        <s:beanshellinput s:syntactictype="'text/plain'">strand</s:beanshellinput>
        <s:beanshellinput s:syntactictype="'text/plain'">chromStart</s:beanshellinput>
        <s:beanshellinput s:syntactictype="'text/plain'">chrom</s:beanshellinput>
        <s:beanshellinput s:syntactictype="'text/plain'">snpID</s:beanshellinput>
      </s:beanshellinputlist>
      <s:beanshelloutputlist>
        <s:beanshelloutput s:syntactictype="'text/plain'">result</s:beanshelloutput>
      </s:beanshelloutputlist>
    </s:beanshell>
    <s:iterationstrategy>
      <i:dot xmlns:i="http://org.embl.ebi.escience/xscufliteration/0.1beta10">
        <i:iterator name="strand" />
        <i:iterator name="chromStart" />
        <i:iterator name="chrom" />
        <i:iterator name="snpID" />
      </i:dot>
    </s:iterationstrategy>
  </s:processor>
  <s:processor name="mouseSubtract1000">
    <s:beanshell>
      <s:scriptvalue>i = Integer.parseInt(input);
output = ""+(i-1000);</s:scriptvalue>
      <s:beanshellinputlist>
        <s:beanshellinput s:syntactictype="'text/plain'">input</s:beanshellinput>
      </s:beanshellinputlist>
      <s:beanshelloutputlist>
        <s:beanshelloutput s:syntactictype="'text/plain'">output</s:beanshelloutput>
      </s:beanshelloutputlist>
    </s:beanshell>
  </s:processor>
  <s:processor name="ratAdd1000">
    <s:beanshell>
      <s:scriptvalue>i = Integer.parseInt(input);
output = ""+(i+1000);</s:scriptvalue>
      <s:beanshellinputlist>
        <s:beanshellinput s:syntactictype="'text/plain'">input</s:beanshellinput>
      </s:beanshellinputlist>
      <s:beanshelloutputlist>
        <s:beanshelloutput s:syntactictype="'text/plain'">output</s:beanshelloutput>
      </s:beanshelloutputlist>
    </s:beanshell>
  </s:processor>
  <s:processor name="humanGene" boring="true">
    <s:stringconstant>ENSG00000086506</s:stringconstant>
  </s:processor>
  <s:processor name="ratSubtract1000">
    <s:beanshell>
      <s:scriptvalue>i = Integer.parseInt(input);
output = ""+(i-1000);</s:scriptvalue>
      <s:beanshellinputlist>
        <s:beanshellinput s:syntactictype="'text/plain'">input</s:beanshellinput>
      </s:beanshellinputlist>
      <s:beanshelloutputlist>
        <s:beanshelloutput s:syntactictype="'text/plain'">output</s:beanshelloutput>
      </s:beanshelloutputlist>
    </s:beanshell>
  </s:processor>
  <s:processor name="mouseAdd1000">
    <s:beanshell>
      <s:scriptvalue>i = Integer.parseInt(inout);
output = ""+(i+1000);</s:scriptvalue>
      <s:beanshellinputlist>
        <s:beanshellinput s:syntactictype="'text/plain'">inout</s:beanshellinput>
      </s:beanshellinputlist>
      <s:beanshelloutputlist>
        <s:beanshelloutput s:syntactictype="'text/plain'">output</s:beanshelloutput>
      </s:beanshelloutputlist>
    </s:beanshell>
  </s:processor>
  <s:processor name="ratReport">
    <s:beanshell>
      <s:scriptvalue>result = strand + "," + chromStart + "," + chrom + "," + snpID;</s:scriptvalue>
      <s:beanshellinputlist>
        <s:beanshellinput s:syntactictype="'text/plain'">strand</s:beanshellinput>
        <s:beanshellinput s:syntactictype="'text/plain'">chrom</s:beanshellinput>
        <s:beanshellinput s:syntactictype="'text/plain'">chromStart</s:beanshellinput>
        <s:beanshellinput s:syntactictype="'text/plain'">snpID</s:beanshellinput>
      </s:beanshellinputlist>
      <s:beanshelloutputlist>
        <s:beanshelloutput s:syntactictype="'text/plain'">report</s:beanshelloutput>
      </s:beanshelloutputlist>
    </s:beanshell>
    <s:iterationstrategy>
      <i:dot xmlns:i="http://org.embl.ebi.escience/xscufliteration/0.1beta10">
        <i:iterator name="strand" />
        <i:iterator name="chrom" />
        <i:iterator name="chromStart" />
        <i:iterator name="snpID" />
      </i:dot>
    </s:iterationstrategy>
  </s:processor>
  <s:processor name="seqFormat">
    <s:beanshell>
      <s:scriptvalue>fastaSeq = "&gt;" + seqID + '\n' + seq + '\n';</s:scriptvalue>
      <s:beanshellinputlist>
        <s:beanshellinput s:syntactictype="'text/plain'">seqID</s:beanshellinput>
        <s:beanshellinput s:syntactictype="'text/plain'">seq</s:beanshellinput>
      </s:beanshellinputlist>
      <s:beanshelloutputlist>
        <s:beanshelloutput s:syntactictype="'text/plain'">fastaSeq</s:beanshelloutput>
      </s:beanshelloutputlist>
    </s:beanshell>
    <s:iterationstrategy>
      <i:dot xmlns:i="http://org.embl.ebi.escience/xscufliteration/0.1beta10">
        <i:iterator name="seqID" />
        <i:iterator name="seq" />
      </i:dot>
    </s:iterationstrategy>
  </s:processor>
  <s:processor name="ratSeqFormat">
    <s:beanshell>
      <s:scriptvalue>fastaSeq = "&gt;" + seqID + '\n' + seq + '\n';</s:scriptvalue>
      <s:beanshellinputlist>
        <s:beanshellinput s:syntactictype="'text/plain'">seqID</s:beanshellinput>
        <s:beanshellinput s:syntactictype="'text/plain'">seq</s:beanshellinput>
      </s:beanshellinputlist>
      <s:beanshelloutputlist>
        <s:beanshelloutput s:syntactictype="'text/plain'">fastaSeq</s:beanshelloutput>
      </s:beanshelloutputlist>
    </s:beanshell>
    <s:iterationstrategy>
      <i:dot xmlns:i="http://org.embl.ebi.escience/xscufliteration/0.1beta10">
        <i:iterator name="seqID" />
        <i:iterator name="seq" />
      </i:dot>
    </s:iterationstrategy>
  </s:processor>
  <s:processor name="mouseSeqFormat">
    <s:beanshell>
      <s:scriptvalue>fastaSeq = "&gt;" + seqID + '\n' + seq + '\n';</s:scriptvalue>
      <s:beanshellinputlist>
        <s:beanshellinput s:syntactictype="'text/plain'">seqID</s:beanshellinput>
        <s:beanshellinput s:syntactictype="'text/plain'">seq</s:beanshellinput>
      </s:beanshellinputlist>
      <s:beanshelloutputlist>
        <s:beanshelloutput s:syntactictype="'text/plain'">fastaSeq</s:beanshelloutput>
      </s:beanshelloutputlist>
    </s:beanshell>
    <s:iterationstrategy>
      <i:dot xmlns:i="http://org.embl.ebi.escience/xscufliteration/0.1beta10">
        <i:iterator name="seqID" />
        <i:iterator name="seq" />
      </i:dot>
    </s:iterationstrategy>
  </s:processor>
  <s:processor name="ratSeq">
    <s:description>Rattus norvegicus genes (RGSC3.4)</s:description>
    <s:biomart>
      <biomart:MartQuery xmlns:biomart="http://org.embl.ebi.escience/xscufl-biomart/0.1alpha">
        <biomart:MartService location="http://www.biomart.org/biomart/martservice" />
        <biomart:MartDataset displayName="Rattus norvegicus genes (RGSC3.4)" name="rnorvegicus_gene_ensembl" type="TableSet" initialBatchSize="100" maximumBatchSize="50000" visible="false">
          <biomart:MartURLLocation default="1" displayName="ENSEMBL 39 (SANGER)" host="www.biomart.org" name="ensembl" port="80" serverVirtualSchema="default" virtualSchema="default" visible="1" />
        </biomart:MartDataset>
        <biomart:Query virtualSchemaName="default" count="0">
          <biomart:Dataset name="rnorvegicus_gene_ensembl">
            <biomart:Filter name="ensembl_gene_id" value="" list="true" />
          </biomart:Dataset>
          <biomart:Dataset name="rnorvegicus_genomic_sequence">
            <biomart:Attribute name="coding" />
          </biomart:Dataset>
          <biomart:Dataset name="rnorvegicus_gene_ensembl_structure">
            <biomart:Attribute name="gene_stable_id" />
          </biomart:Dataset>
        </biomart:Query>
      </biomart:MartQuery>
    </s:biomart>
  </s:processor>
  <s:processor name="mouseSeq">
    <s:description>Mus musculus genes (NCBIM36)</s:description>
    <s:biomart>
      <biomart:MartQuery xmlns:biomart="http://org.embl.ebi.escience/xscufl-biomart/0.1alpha">
        <biomart:MartService location="http://www.biomart.org/biomart/martservice" />
        <biomart:MartDataset displayName="Mus musculus genes (NCBIM36)" name="mmusculus_gene_ensembl" type="TableSet" initialBatchSize="100" maximumBatchSize="50000" visible="false">
          <biomart:MartURLLocation default="1" displayName="ENSEMBL 39 (SANGER)" host="www.biomart.org" name="ensembl" port="80" serverVirtualSchema="default" virtualSchema="default" visible="1" />
        </biomart:MartDataset>
        <biomart:Query virtualSchemaName="default" count="0">
          <biomart:Dataset name="mmusculus_gene_ensembl">
            <biomart:Filter name="ensembl_gene_id" value="" list="true" />
          </biomart:Dataset>
          <biomart:Dataset name="mmusculus_genomic_sequence">
            <biomart:Attribute name="coding" />
          </biomart:Dataset>
          <biomart:Dataset name="mmusculus_gene_ensembl_structure">
            <biomart:Attribute name="gene_stable_id" />
          </biomart:Dataset>
        </biomart:Query>
      </biomart:MartQuery>
    </s:biomart>
  </s:processor>
  <s:processor name="findHomologues">
    <s:description>Homo sapiens genes (NCBI36)</s:description>
    <s:biomart>
      <biomart:MartQuery xmlns:biomart="http://org.embl.ebi.escience/xscufl-biomart/0.1alpha">
        <biomart:MartService location="http://www.biomart.org/biomart/martservice" />
        <biomart:MartDataset displayName="Homo sapiens genes (NCBI36)" name="hsapiens_gene_ensembl" type="TableSet" initialBatchSize="100" maximumBatchSize="50000" visible="false">
          <biomart:MartURLLocation default="1" displayName="ENSEMBL 39 (SANGER)" host="www.biomart.org" name="ensembl" port="80" serverVirtualSchema="default" virtualSchema="default" visible="1" />
        </biomart:MartDataset>
        <biomart:Query virtualSchemaName="default" count="0">
          <biomart:Dataset name="hsapiens_gene_ensembl">
            <biomart:Attribute name="mouse_ensembl_gene" />
            <biomart:Attribute name="rat_ensembl_gene" />
            <biomart:Filter name="ensembl_gene_id" value="" list="true" />
          </biomart:Dataset>
        </biomart:Query>
      </biomart:MartQuery>
    </s:biomart>
  </s:processor>
  <s:processor name="rnorvegicus_snp">
    <s:description>Rattus norvegicus SNPs (dbSNP 125)</s:description>
    <s:biomart>
      <biomart:MartQuery xmlns:biomart="http://org.embl.ebi.escience/xscufl-biomart/0.1alpha">
        <biomart:MartService location="http://www.biomart.org/biomart/martservice" />
        <biomart:MartDataset displayName="Rattus norvegicus SNPs (dbSNP 125)" name="rnorvegicus_snp" type="TableSet" initialBatchSize="100" maximumBatchSize="50000" visible="false">
          <biomart:MartURLLocation default="0" displayName="SNP 39 (SANGER)" host="www.biomart.org" name="snp" port="80" serverVirtualSchema="default" virtualSchema="default" visible="1" />
        </biomart:MartDataset>
        <biomart:Query virtualSchemaName="default" count="0">
          <biomart:Dataset name="rnorvegicus_snp">
            <biomart:Attribute name="chr_name" />
            <biomart:Attribute name="chrom_start" />
            <biomart:Attribute name="chrom_strand" />
            <biomart:Attribute name="refsnp_id" />
            <biomart:Filter name="chr_name" value="1" />
            <biomart:Filter name="chrom_end" value="" />
            <biomart:Filter name="chrom_start" value="" />
          </biomart:Dataset>
        </biomart:Query>
      </biomart:MartQuery>
    </s:biomart>
    <s:iterationstrategy>
      <i:dot xmlns:i="http://org.embl.ebi.escience/xscufliteration/0.1beta10">
        <i:iterator name="rnorvegicus_snp.chr_name_filter" />
        <i:iterator name="rnorvegicus_snp.chrom_end_filter" />
        <i:iterator name="rnorvegicus_snp.chrom_start_filter" />
      </i:dot>
    </s:iterationstrategy>
  </s:processor>
  <s:processor name="findRatGene">
    <s:description>Rattus norvegicus genes (RGSC3.4)</s:description>
    <s:biomart>
      <biomart:MartQuery xmlns:biomart="http://org.embl.ebi.escience/xscufl-biomart/0.1alpha">
        <biomart:MartService location="http://www.biomart.org/biomart/martservice" />
        <biomart:MartDataset displayName="Rattus norvegicus genes (RGSC3.4)" name="rnorvegicus_gene_ensembl" type="TableSet" initialBatchSize="100" maximumBatchSize="50000" visible="false">
          <biomart:MartURLLocation default="1" displayName="ENSEMBL 39 (SANGER)" host="www.biomart.org" name="ensembl" port="80" serverVirtualSchema="default" virtualSchema="default" visible="1" />
        </biomart:MartDataset>
        <biomart:Query virtualSchemaName="default" count="0">
          <biomart:Dataset name="rnorvegicus_gene_ensembl">
            <biomart:Attribute name="chromosome_name" />
            <biomart:Attribute name="end_position" />
            <biomart:Attribute name="start_position" />
            <biomart:Filter name="ensembl_gene_id" value="" list="true" />
          </biomart:Dataset>
        </biomart:Query>
      </biomart:MartQuery>
    </s:biomart>
  </s:processor>
  <s:processor name="findMouseGene">
    <s:description>Mus musculus genes (NCBIM36)</s:description>
    <s:biomart>
      <biomart:MartQuery xmlns:biomart="http://org.embl.ebi.escience/xscufl-biomart/0.1alpha">
        <biomart:MartService location="http://www.biomart.org/biomart/martservice" />
        <biomart:MartDataset displayName="Mus musculus genes (NCBIM36)" name="mmusculus_gene_ensembl" type="TableSet" initialBatchSize="100" maximumBatchSize="50000" visible="false">
          <biomart:MartURLLocation default="1" displayName="ENSEMBL 39 (SANGER)" host="www.biomart.org" name="ensembl" port="80" serverVirtualSchema="default" virtualSchema="default" visible="1" />
        </biomart:MartDataset>
        <biomart:Query virtualSchemaName="default" count="0">
          <biomart:Dataset name="mmusculus_gene_ensembl">
            <biomart:Attribute name="chromosome_name" />
            <biomart:Attribute name="start_position" />
            <biomart:Attribute name="end_position" />
            <biomart:Filter name="ensembl_gene_id" value="" list="true" />
          </biomart:Dataset>
        </biomart:Query>
      </biomart:MartQuery>
    </s:biomart>
  </s:processor>
  <s:processor name="mmusculus_snp">
    <s:description>Mus musculus SNPs (dbSNP 126; Sanger)</s:description>
    <s:biomart>
      <biomart:MartQuery xmlns:biomart="http://org.embl.ebi.escience/xscufl-biomart/0.1alpha">
        <biomart:MartService location="http://www.biomart.org/biomart/martservice" />
        <biomart:MartDataset displayName="Mus musculus SNPs (dbSNP 126; Sanger)" name="mmusculus_snp" type="TableSet" initialBatchSize="100" maximumBatchSize="50000" visible="false">
          <biomart:MartURLLocation default="0" displayName="SNP 39 (SANGER)" host="www.biomart.org" name="snp" port="80" serverVirtualSchema="default" virtualSchema="default" visible="1" />
        </biomart:MartDataset>
        <biomart:Query virtualSchemaName="default" count="0">
          <biomart:Dataset name="mmusculus_snp">
            <biomart:Attribute name="chr_name" />
            <biomart:Attribute name="chrom_strand" />
            <biomart:Attribute name="chrom_start" />
            <biomart:Attribute name="refsnp_id" />
            <biomart:Filter name="chr_name" value="1" />
            <biomart:Filter name="chrom_end" value="" />
            <biomart:Filter name="chrom_start" value="" />
          </biomart:Dataset>
        </biomart:Query>
      </biomart:MartQuery>
    </s:biomart>
    <s:iterationstrategy>
      <i:dot xmlns:i="http://org.embl.ebi.escience/xscufliteration/0.1beta10">
        <i:iterator name="mmusculus_snp.chr_name_filter" />
        <i:iterator name="mmusculus_snp.chrom_end_filter" />
        <i:iterator name="mmusculus_snp.chrom_start_filter" />
      </i:dot>
    </s:iterationstrategy>
  </s:processor>
  <s:processor name="NCBI36">
    <s:description>Homo sapiens genes (NCBI36)</s:description>
    <s:biomart>
      <biomart:MartQuery xmlns:biomart="http://org.embl.ebi.escience/xscufl-biomart/0.1alpha">
        <biomart:MartService location="http://www.biomart.org/biomart/martservice" />
        <biomart:MartDataset displayName="Homo sapiens genes (NCBI36)" name="hsapiens_gene_ensembl" type="TableSet" initialBatchSize="100" maximumBatchSize="50000" visible="false">
          <biomart:MartURLLocation default="1" displayName="ENSEMBL 38 (SANGER)" host="www.biomart.org" name="ensembl" port="80" serverVirtualSchema="default" virtualSchema="default" visible="1" />
        </biomart:MartDataset>
        <biomart:Query virtualSchemaName="default" count="0">
          <biomart:Dataset name="hsapiens_gene_ensembl">
            <biomart:Attribute name="chromosome_name" />
            <biomart:Attribute name="end_position" />
            <biomart:Attribute name="start_position" />
            <biomart:Filter name="ensembl_gene_id" value="" list="true" />
          </biomart:Dataset>
        </biomart:Query>
      </biomart:MartQuery>
    </s:biomart>
  </s:processor>
  <s:processor name="humanSeq">
    <s:description>Homo sapiens genes (NCBI36)</s:description>
    <s:biomart>
      <biomart:MartQuery xmlns:biomart="http://org.embl.ebi.escience/xscufl-biomart/0.1alpha">
        <biomart:MartService location="http://www.biomart.org/biomart/martservice" />
        <biomart:MartDataset displayName="Homo sapiens genes (NCBI36)" name="hsapiens_gene_ensembl" type="TableSet" initialBatchSize="100" maximumBatchSize="50000" visible="false">
          <biomart:MartURLLocation default="1" displayName="ENSEMBL 39 (SANGER)" host="www.biomart.org" name="ensembl" port="80" serverVirtualSchema="default" virtualSchema="default" visible="1" />
        </biomart:MartDataset>
        <biomart:Query virtualSchemaName="default" count="0">
          <biomart:Dataset name="hsapiens_genomic_sequence">
            <biomart:Attribute name="coding" />
          </biomart:Dataset>
          <biomart:Dataset name="hsapiens_gene_ensembl_structure">
            <biomart:Attribute name="gene_stable_id" />
          </biomart:Dataset>
          <biomart:Dataset name="hsapiens_gene_ensembl">
            <biomart:Filter name="ensembl_gene_id" value="" list="true" />
          </biomart:Dataset>
        </biomart:Query>
      </biomart:MartQuery>
    </s:biomart>
  </s:processor>
  <s:processor name="dbSNP">
    <s:description>Homo sapiens SNPs (dbSNP 125)</s:description>
    <s:biomart>
      <biomart:MartQuery xmlns:biomart="http://org.embl.ebi.escience/xscufl-biomart/0.1alpha">
        <biomart:MartService location="http://www.biomart.org/biomart/martservice" />
        <biomart:MartDataset displayName="Homo sapiens SNPs (dbSNP 125)" name="hsapiens_snp" type="TableSet" initialBatchSize="100" maximumBatchSize="50000" visible="false">
          <biomart:MartURLLocation default="0" displayName="SNP 38 (SANGER)" host="www.biomart.org" name="snp" port="80" serverVirtualSchema="default" virtualSchema="default" visible="1" />
        </biomart:MartDataset>
        <biomart:Query virtualSchemaName="default" count="0">
          <biomart:Dataset name="hsapiens_snp">
            <biomart:Attribute name="refsnp_id" />
            <biomart:Attribute name="chr_name" />
            <biomart:Attribute name="chrom_start" />
            <biomart:Attribute name="chrom_strand" />
            <biomart:Filter name="chr_name" value="1" />
            <biomart:Filter name="chrom_end" value="" />
            <biomart:Filter name="chrom_start" value="" />
          </biomart:Dataset>
        </biomart:Query>
      </biomart:MartQuery>
    </s:biomart>
    <s:iterationstrategy>
      <i:dot xmlns:i="http://org.embl.ebi.escience/xscufliteration/0.1beta10">
        <i:iterator name="hsapiens_snp.chrom_start_filter" />
        <i:iterator name="hsapiens_snp.chrom_end_filter" />
        <i:iterator name="hsapiens_snp.chr_name_filter" />
      </i:dot>
    </s:iterationstrategy>
  </s:processor>
  <s:processor name="CreateReport">
    <s:beanshell>
      <s:scriptvalue>result = strand + "," + chromStart + "," + chrom + "," + snpID;</s:scriptvalue>
      <s:beanshellinputlist>
        <s:beanshellinput s:syntactictype="'text/plain'">strand</s:beanshellinput>
        <s:beanshellinput s:syntactictype="'text/plain'">chromStart</s:beanshellinput>
        <s:beanshellinput s:syntactictype="'text/plain'">chrom</s:beanshellinput>
        <s:beanshellinput s:syntactictype="'text/plain'">snpID</s:beanshellinput>
      </s:beanshellinputlist>
      <s:beanshelloutputlist>
        <s:beanshelloutput s:syntactictype="'text/plain'">result</s:beanshelloutput>
      </s:beanshelloutputlist>
    </s:beanshell>
    <s:iterationstrategy>
      <i:dot xmlns:i="http://org.embl.ebi.escience/xscufliteration/0.1beta10">
        <i:iterator name="snpID" />
        <i:iterator name="chrom" />
        <i:iterator name="chromStart" />
        <i:iterator name="strand" />
      </i:dot>
    </s:iterationstrategy>
  </s:processor>
  <s:processor name="nucleotide_Alignment">
    <s:defaults>
      <s:default name="in1">ALIGNED</s:default>
      <s:default name="in3">15.00</s:default>
      <s:default name="in4">6.66</s:default>
      <s:default name="in5">15.00</s:default>
      <s:default name="in6">6.66</s:default>
      <s:default name="in7">30</s:default>
      <s:default name="in8">0.50</s:default>
      <s:default name="in9">iub</s:default>
      <s:default name="in10">iub</s:default>
    </s:defaults>
    <s:arbitrarywsdl>
      <s:wsdl>http://www.vbi.vt.edu/~pathport/services/wsdls/beta/msagt.wsdl</s:wsdl>
      <s:operation>nucleotide_Alignment</s:operation>
    </s:arbitrarywsdl>
  </s:processor>
  <s:processor name="CreateFasta">
    <s:workflow>
      <s:scufl version="0.2" log="0">
        <s:workflowdescription lsid="urn:lsid:www.mygrid.org.uk:operation:Z6160QD7NK30" author="" title="" />
        <s:processor name="mergeMouse">
          <s:local>org.embl.ebi.escience.scuflworkers.java.StringListMerge</s:local>
        </s:processor>
        <s:processor name="mergeRat">
          <s:local>org.embl.ebi.escience.scuflworkers.java.StringListMerge</s:local>
        </s:processor>
        <s:processor name="mergeHuman">
          <s:local>org.embl.ebi.escience.scuflworkers.java.StringListMerge</s:local>
        </s:processor>
        <s:processor name="combineAll">
          <s:local>org.embl.ebi.escience.scuflworkers.java.StringListMerge</s:local>
          <s:mergemode input="stringlist" mode="merge" />
        </s:processor>
        <s:link source="humanList" sink="mergeHuman:stringlist" />
        <s:link source="mouseList" sink="mergeMouse:stringlist" />
        <s:link source="ratList" sink="mergeRat:stringlist" />
        <s:link source="combineAll:concatenated" sink="fasta" />
        <s:link source="mergeHuman:concatenated" sink="combineAll:stringlist" />
        <s:link source="mergeMouse:concatenated" sink="combineAll:stringlist" />
        <s:link source="mergeRat:concatenated" sink="combineAll:stringlist" />
        <s:source name="ratList" />
        <s:source name="humanList" />
        <s:source name="mouseList" />
        <s:sink name="fasta" />
      </s:scufl>
    </s:workflow>
  </s:processor>
  <s:link source="Add1000:output" sink="dbSNP:hsapiens_snp.chrom_end_filter" />
  <s:link source="CreateFasta:fasta" sink="nucleotide_Alignment:in0" />
  <s:link source="CreateFasta:fasta" sink="fastaSequences" />
  <s:link source="humanSeq:hsapiens_gene_ensembl_structure.gene_stable_id" sink="seqFormat:seqID" />
  <s:link source="humanSeq:hsapiens_genomic_sequence.coding" sink="seqFormat:seq" />
  <s:link source="mmusculus_snp:mmusculus_snp.chr_name" sink="mouseReport:chrom" />
  <s:link source="mmusculus_snp:mmusculus_snp.chrom_start" sink="mouseReport:chromStart" />
  <s:link source="mmusculus_snp:mmusculus_snp.chrom_strand" sink="mouseReport:strand" />
  <s:link source="mmusculus_snp:mmusculus_snp.refsnp_id" sink="mouseReport:snpID" />
  <s:link source="mouseAdd1000:output" sink="mmusculus_snp:mmusculus_snp.chrom_end_filter" />
  <s:link source="mouseSeq:mmusculus_gene_ensembl_structure.gene_stable_id" sink="mouseSeqFormat:seqID" />
  <s:link source="mouseSeq:mmusculus_genomic_sequence.coding" sink="mouseSeqFormat:seq" />
  <s:link source="mouseSeqFormat:fastaSeq" sink="CreateFasta:mouseList" />
  <s:link source="mouseSubtract1000:output" sink="mmusculus_snp:mmusculus_snp.chrom_start_filter" />
  <s:link source="ratSeqFormat:fastaSeq" sink="CreateFasta:ratList" />
  <s:link source="seqFormat:fastaSeq" sink="CreateFasta:humanList" />
  <s:link source="CreateReport:result" sink="HSapSNPs" />
  <s:link source="NCBI36:hsapiens_gene_ensembl.chromosome_name" sink="dbSNP:hsapiens_snp.chr_name_filter" />
  <s:link source="NCBI36:hsapiens_gene_ensembl.end_position" sink="Add1000:input" />
  <s:link source="NCBI36:hsapiens_gene_ensembl.start_position" sink="Subtract1000:input" />
  <s:link source="Subtract1000:output" sink="dbSNP:hsapiens_snp.chrom_start_filter" />
  <s:link source="dbSNP:hsapiens_snp.chr_name" sink="CreateReport:chrom" />
  <s:link source="dbSNP:hsapiens_snp.chrom_start" sink="CreateReport:chromStart" />
  <s:link source="dbSNP:hsapiens_snp.chrom_strand" sink="CreateReport:strand" />
  <s:link source="dbSNP:hsapiens_snp.refsnp_id" sink="CreateReport:snpID" />
  <s:link source="findHomologues:hsapiens_gene_ensembl.mouse_ensembl_gene" sink="findMouseGene:mmusculus_gene_ensembl.ensembl_gene_id_filter" />
  <s:link source="findHomologues:hsapiens_gene_ensembl.mouse_ensembl_gene" sink="mouseSeq:mmusculus_gene_ensembl.ensembl_gene_id_filter" />
  <s:link source="findHomologues:hsapiens_gene_ensembl.rat_ensembl_gene" sink="findRatGene:rnorvegicus_gene_ensembl.ensembl_gene_id_filter" />
  <s:link source="findHomologues:hsapiens_gene_ensembl.rat_ensembl_gene" sink="ratSeq:rnorvegicus_gene_ensembl.ensembl_gene_id_filter" />
  <s:link source="findMouseGene:mmusculus_gene_ensembl.chromosome_name" sink="mmusculus_snp:mmusculus_snp.chr_name_filter" />
  <s:link source="findMouseGene:mmusculus_gene_ensembl.end_position" sink="mouseAdd1000:inout" />
  <s:link source="findMouseGene:mmusculus_gene_ensembl.start_position" sink="mouseSubtract1000:input" />
  <s:link source="findRatGene:rnorvegicus_gene_ensembl.chromosome_name" sink="rnorvegicus_snp:rnorvegicus_snp.chr_name_filter" />
  <s:link source="findRatGene:rnorvegicus_gene_ensembl.end_position" sink="ratAdd1000:input" />
  <s:link source="findRatGene:rnorvegicus_gene_ensembl.start_position" sink="ratSubtract1000:input" />
  <s:link source="humanGene:value" sink="NCBI36:hsapiens_gene_ensembl.ensembl_gene_id_filter" />
  <s:link source="humanGene:value" sink="findHomologues:hsapiens_gene_ensembl.ensembl_gene_id_filter" />
  <s:link source="humanGene:value" sink="humanSeq:hsapiens_gene_ensembl.ensembl_gene_id_filter" />
  <s:link source="mouseReport:result" sink="MMusSNPs" />
  <s:link source="mouseSeqFormat:fastaSeq" sink="MMusFasta" />
  <s:link source="nucleotide_Alignment:nucleotide_AlignmentReturn" sink="MultipleAlignment" />
  <s:link source="ratAdd1000:output" sink="rnorvegicus_snp:rnorvegicus_snp.chrom_end_filter" />
  <s:link source="ratSeq:rnorvegicus_gene_ensembl_structure.gene_stable_id" sink="ratSeqFormat:seqID" />
  <s:link source="ratSeq:rnorvegicus_genomic_sequence.coding" sink="ratSeqFormat:seq" />
  <s:link source="ratReport:report" sink="RNorSNPs" />
  <s:link source="ratSeqFormat:fastaSeq" sink="RNorFasta" />
  <s:link source="ratSubtract1000:output" sink="rnorvegicus_snp:rnorvegicus_snp.chrom_start_filter" />
  <s:link source="rnorvegicus_snp:rnorvegicus_snp.chr_name" sink="ratReport:chrom" />
  <s:link source="rnorvegicus_snp:rnorvegicus_snp.chrom_start" sink="ratReport:chromStart" />
  <s:link source="rnorvegicus_snp:rnorvegicus_snp.chrom_strand" sink="ratReport:strand" />
  <s:link source="rnorvegicus_snp:rnorvegicus_snp.refsnp_id" sink="ratReport:snpID" />
  <s:link source="seqFormat:fastaSeq" sink="HSapFasta" />
  <s:sink name="HSapSNPs" />
  <s:sink name="MMusSNPs" />
  <s:sink name="RNorSNPs" />
  <s:sink name="HSapFasta" />
  <s:sink name="RNorFasta" />
  <s:sink name="MMusFasta" />
  <s:sink name="MultipleAlignment">
    <s:metadata>
      <s:mimeTypes>
        <s:mimeType>text/xml</s:mimeType>
      </s:mimeTypes>
    </s:metadata>
  </s:sink>
  <s:sink name="fastaSequences" />
</s:scufl>

