A brief introduction to SSEARCH
SSEARCH performs a rigorous Smith-Waterman alignment between a protein sequence and a protein database, or with DNA sequence to another
DNA database (not available in NPS@).
So, for example, you can use SSEARCH to find similar sequences in a database built
from an ACNUC query or
after a PATTINPROT search.
SSEARCH parameters are not currently available for the user.
By default the number of description and the number of alignment are set to 500.
The expected threshold is set to 10.0.
The comparison matrix is BLOSUM50.
The gap opening penalty is -12 and gap extending penalty is -2.
To retrieve these informations see the line "Information and statistics : [SSEARCH]" in
NPS@ SSEARCH result file.
NPS@ SSEARCH output example
The NPS@ SSEARCH output is divided into three parts.
A graphic indicating the similarity found along the query sequence. It's computed with SSEARCH alignment. This graphic can
be see in MPSA/ANTHEPROT to determine interesting boundaries in the query sequence.
A form to select subject sequences. You can select sequences by indicated the maximal E threshold to do it (by default
NPS@ sets it to 1e-6). You can also select subject sequences matching a particular region of the
It's the SSEARCH description block in a 'HTMLized' form.
You can see :
A checkbox to select/unselect subject sequence.
The NPSA link allows you to apply
NPS@ methods on the corresponding sequence after it is extracted from the query database.
The database link to retrieve the database entry.
The alignment link (on expected value) to see SSEARCH alignment between the query sequence and the current subject sequence.
The number of displayed sequences.
The number of selected sequences.
You have :
An extract form to extract selected sequence and make a database. You can extract full or partial sequences.
The full extract is made from the database.
The partial extract, is made for subject sequence that match the query one in the indicated region. The extracted sequence
is retrieved from the SSEARCH alignment. You can set the minimal length of the extracted sequence.