Services - BLASTatlas (SOAP) [Provided by: Center for Biological Sequence Analysis (CBS)]

Provider:
Center for Biological Sequence Analysis (CBS)

Location:
Denmark

Submitter / Source:
Franck Tanoh (about 1 year ago)

Base URL:
http://ws.cbs.dtu.dk/cgi-bin/soap/ws/quasi.cgi

WSDL Location:
http://www.cbs.dtu.dk/ws/BLASTatlas/BLASTatlas_1_0_ws2.wsdl(download last cached WSDL file)

Documentation URL(s): Franck Tanoh (about 1 year ago)http://www.cbs.dtu.dk/ws/ws.php?entry=BLASTatlas

pfhallin (9 months ago)http://www.cbs.dtu.dk/ws/BLASTatlas

This service creates maps of genome homology of a list of sequences against a reference
genome using either blastp, blastn, tblastn, or blastx.

In the following, ‘database’ refers to the sequences being searched whereas
‘reference’ refers to the genome on to which the search is mapped. A map will always
have the same length as the reference genome on to which the search is mapped.

For DNA blast, there will be a 1:1 map of alignment:genome, and values can attain either
0 (gap/mismatch) or 1 (match). For protein searches, there will be a 1:3 map of alignment:genome.
For a given amino acid match, all three bases in the genome that encodes this residue, will
attain the value ‘1’. For conservative mismatches all three positions will attain the
value of 0.5, whereas gaps and non-conservative matches will attain the value ‘0’. The
different methods require different input, as described in the table below:

The result of the service is a PDF and a PostScript document containing the BLASTatlas

1. runBLASTatlas
This operation reads in a reference genome and its annotation together with a
number of databases to search again. A standard BLASTatlas is constructed, returned
as a PostScript files. This is here you want to get started!

INPUT
|- main : Main title of the atlas
|- sub : Sub title of the atlas
|- stamp : Type of atlas, printed on lower right corner
|+ window (optional) The window boundaries when zooming
||- begin : Begin of zoom window
|.- end : End of zoom window
|
|- modus : ‘circle’ or ‘linear’: Defines show the lanes should be printed
|+ ann’ (optional, array) This element is used to label annotations individually
|| on the perimeter/top of the atlas. Be aware, that the
|| annotation must be present in the ‘atlasannotations’
|| section (see below), in order for the coordinates to
|| be known to GeneWiz (see the provided Perl example)
||- type : Feature type, e.g. tRNA, rRNA, CDS, misc etc.
||- label : Label, must be the same as given in ‘atlasannotations’
|.- dir : Either ‘pos’ or ‘neg’
|
|- DNAproperties (optional, array) These elements allows a list of DNA properties (like
| Intrinsic Curvature, Stacking Energy etc.) to be added to the lower/inner
| part of the atlas. Please refer to the XSD to see which properties are
| supported, or see the Genome Atlas documentation for details on these
| properties: http://www.cbs.dtu.dk/ws/GenomeAtlas
|
|+ customMap (optional, array) This element allows you to add a custom
|| map, a list of floating point numbers, separated by comma.
|| Tools like SIDDbase (http://www.cbs.dtu.dk/ws/SIDDbase) returns
|| calculations in this format, and can be easily inserted to the atlas.
|| Each of the custom lanes submitted will be printed between the
|| DNAproperties (innermost/bottom) and the BLASTlanes (outermost/top)
||- boxfilter : (optional) The smoothing of the data
||- legend : Legend text for this lane
|| (there is a choice between ‘byrange’ and ‘byaverage’)
||+ byrange
|||- bottom : Lower boundary of the color range
||.- top : Upper boundary of the color range
||+ byaverage
||.- stddev : Color range will go this many standard deviations around
|| the average
|.+ color
| |+ from Lower color of range
| ||- r : Red component
| ||- g : Green component
| |.- b : Blue component
| |+ to Upper color of range
| ||- r : Red component
| ||- g : Green component
| |.- b : Blue component
| .+ via : (optional), mid part of the color range
| |- r : Red component
| |- g : Green component
| .- b : Blue component
|
|+ reference (contains information about the reference genome, its annotations etc.)
||- genome : Contains the genome sequence of the reference genome as one continuos string
||+ atlasfeature (These elements define the genome annotations that will be printed on the atlas)
||.+ feature (array)
|| |- type : Either CDS, rRNA, tRNA
|| |- begin : Start position of annotations
|| |- end : Stop position of annotation
|| |- dir : Strand (‘-‘ or ‘+’ )
|| .- label : Label of the annotation. This is mostly relevant when annotations
|| are marked with the optional ‘ann’ element, mentioned above.
||
||+ blastfeature : This is the annotation of the reference genome, used in the blast
| | search to construct the homology map
| .+ feature (array)
| |- proteins : The translation of the DNA region. Mandatory if lanes
| | are included using programs blastp or tblastn.
| |- orfs : The protein coding sequences. Optional (but prioritized) when lanes
| | are present using blastn or blastx. If ‘orfs’ are left out, the region will be extracted
| | from the genome sequence.
| |- begin : Start of the protein coding region
| |- end : Stop of the protein coding region
| .- dir : Translation direction (+,-)
|
|
.+ db (array) Each element corresponds to a genome/proteme to be included as a lane.
|- legend : Lane legend
|- program : Either ‘blastp’, ‘blastn’,’tblastn’, or ‘blastx’
|+ color
||+ from (Lower color of range)
|||- r : Red component
|||- g : Green component
||.- b : Blue component
||
||+ to’ (Upper color of range)
|||- r’ : Red component
|||- g’ : Green component
||.- b’ : Blue component
||
|.+ via’ (optional), mid part of the color range
| |- r : Red component
| |- g : Green component
| .- b : Blue component
|+ byrange (Ranges is fixed …)
||- bottom : Lower boundary of the color range
|.- top : Upper boundary of the color range
|
|+ byaverage (Range is dynamic …)
|.- stddev : Color range will go this many standard deviations around
| the average
|
| (there is an option of the next two elements ‘orfs’ and ‘proteins’ at this level)
|- orfs : (array) ORFs to be searched (required for blastn or tblastn)
.- proteins : (array) Protein sequences to be searched (required for blastp or blastx)

OUTPUT

|- jobid : The 32 byte identification string of the job
|- datetime : The last timepoint at which the status of the job has changed
|- status : Possible values are QUEUED, ACTIVE, FINISHED, WAITING, REJECTED,
| UNKNOWN JOBID or QUEUE DOWN
.- expire : Amount of hours since your job will expire

2. pollQueue

INPUT
.- jobid : The 32 byte identification string of the job

OUTPUT
| jobid : The 32 byte identification string of the job
| datetime : The last timepoint at which the status of the job has changed
| status : Possible values are QUEUED, ACTIVE, FINISHED, WAITING, REJECTED,
| UNKNOWN JOBID or QUEUE DOWN
. expire : Amount of hours since your job will expire

3. fetchAtlasResult

INTPUT
.- jobid : The 32 byte identification string of the job

OUTPUT
|+ pdf
||- comment : ‘Genewiz output’
||- encoding : ‘base54’
||- MIMEtype : ‘application/pdf’
|.- content : The encoded document
|
.+ ps
|- comment : ‘Genewiz output’
|- encoding : ‘base54’
|- MIMEtype : ‘application/ps’
.- content : The encoded document

pfhallin (9 months ago)This BLASTatlas service creates maps of genome homology of a list of sequences against a reference
genome using either blastp, blastn, tblastn, or blastx.

The resolution is per-residue or per nucleotide depending on the regime of the blast
search: For each annotation in the reference genome, the best hit in the database
genome is found using one of the above algorithms. Each matching or mismatching
residue/nucleotide of the best hit (based on BLAST score) is then mapped back to the
genome sequence, using the coordinates provided in the annotations.
In the following, ‘database’ refers to the sequences being searched whereas
‘reference’ refers to the genome on to which the search is mapped. A map will always
have the same length as the reference genome on to which the search is mapped.
For DNA blast, there will be a 1:1 map of alignment:genome, and values can attain either
0 (gap/mismatch) or 1 (match). For protein searches, there will be a 1:3 map of alignment:genome.
For a given amino acid match, all three bases in the genome that encodes this residue, will
attain the value ‘1’. For conservative mismatches all three positions will attain the
value of 0.5, whereas gaps and non-conservative matches will attain the value ‘0’.

The genome BLASTatlas-a GeneWiz extension for visualization ofwhole-genome homology.Mol Biosyst 4:363-71 (2008) PubMed

Login to add a description License(s): None Login to add license info Cost: No info yet Login to add cost info Usage Conditions: No info yet Login to add usage conditions info Contact Info: None Login to add contact info Publications: for this service. This can be a URI to the publication and/or a DOI. None Login to add publication info Citations: None Login to add a citation