BLASTatlas
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
WSDL Changes
| Date | Description of Changes |
Test #1: BLASTatlas test case
Construction of a protein BLASTatlas using Candidatus Sulcia muelleri GWSS as reference and Candidatus Carsonella ruddii PV as query genome.