Course home > Smith Waterman worked example

We are going to do the same search that we performed with BLAST and Fasta, with the Felis catus cannabinoid receptor 1, to see how similar the results are to those that with got with the previous searches.

Paste your sequence into the main text-entry box in the Smith-Waterman home page at EBI. We shall select the Swissprot protein database against which to perform our search. We can leave the default settings for the moment. Be sure to change the database that we are searching against to Swissprot. Also change the Program type from Normal to Affine to allow the use of affine gap penalties (i.e., gap extension penalties are a fraction of the penalty to open a gap). Click 'Run'.

Show SMITH-WATERMAN QUERY FORM image

Our results will appear in the same window shortly. As with FASTA, the sequence is not filtered for low complexity regions.

The main results page contains a summary table, similar to the FASTA result page. Note the affine gap penalties.

Show SMITH-WATERMAN SUMMARY image

The 'MPsrch Result' button takes you to a page which details the parameters of this run (again, similarly to the FASTA results page). From the text at the top of this file, we can again notice that the Swissprot database that is used in the Smith-Waterman searches is different from both the one used by NCBI and the one used by the FASTA search.

Show SMITH-WATERMAN DETAILS image

However, the results that we get are very similar to those that we get doing the FASTA search. The layout is slightly different. And, as in BLAST, only the portions of the sequence that are actually involved in the match are shown.

Show SMITH-WATERMAN ALIGNMENT image
News
Jul, 2009; ChIPseeqer, a comprehensive framework for analysis of ChIP-seq data developed in the Elemento lab, is now available for download. [More]
Apr, 2009; The BDVal program developed by the Campagne laboratory for MAQC-II is now available from http://bdval.org. The software supports the development and evaluation of predictive biomarker models from high-throughput data. The web site offers binary and source distributions. [More]
Jan, 2009; Twease now supports searching MEDLINE articles by Author, Journal, and Publication Year. Examples for performing these searches can be found in the updated Twease tutorial. [More]

[News Archives] [Mailing List]


Events
Dec 11th; 4:00pm-5:00pm: Institute for Computational Biomedicine Research in Progress Seminar Series - Fabien Campagne; ICB Conference Room - Y.1301
Jan 15th; 4:00pm-5:00pm: Institute for Computational Biomedicine Research in Progress Seminar Series - Lei Shi; ICB Conference Room - Y.1301
Feb 12th; 4:00pm-5:00pm: Institute for Computational Biomedicine Research in Progress Seminar Series - Christopher E. Mason; ICB Conference Room - Y.1301
Mar 12th; 4:00pm-5:00pm: Institute for Computational Biomedicine Research in Progress Seminar Series - Olivier Elemento; ICB Conference Room - Y.1301
Apr 9th; 4:00pm-5:00pm: Institute for Computational Biomedicine Research in Progress Seminar Series - Emre Aksay; ICB Conference Room - Y.1301
May 14th; 4:00pm-5:00pm: Institute for Computational Biomedicine Research in Progress Seminar Series - Jonathan D. Victor; ICB Conference Room - Y.1301
Jun 11th; 4:00pm-5:00pm: Institute for Computational Biomedicine Research in Progress Seminar Series - Harel Weinstein; ICB Conference Room - Y.1301
Jul 9th; 4:00pm-5:00pm: Institute for Computational Biomedicine Research in Progress Seminar Series - Duane Hassane; ICB Conference Room - Y.1301