[ Team LiB ] Previous Section Next Section

profit

profit takes a simple frequency matrix produced by prophecy and uses it to find matches in the input sequence(s) you are searching.

Here is a sample session with profit:

My aligned set of sequences:

% more m.seq
>one
DEVGGEALGRLLVVYPWTQR
>two
DEVGREALGRLLVVYPWTQR
>three
DEVGGEALGRILVVYPWTQR
>four
DEVGGEAAGRVLVVYPWTQR

Make a simple frequency matrix using prophecy:

% prophecy
Creates matrices/profiles from multiple alignments
Input sequence set: m.seq
Profile type
         F : Frequency
         G : Gribskov
         H : Henikoff
Select type [F]:
Enter a name for the profile [mymatrix]:
Enter threshold reporting percentage [75]:
Output file [outfile.prophecy]:

Search using profit:

% profit
Scan a sequence or database with a matrix or profile
Profile or matrix file: outfile.prophecy
Input sequence(s): sw:*
Output file [outfile.profit]:

Mandatory qualifiers:

[-infile] (infile)

Profile or matrix file.

[-sequence] (seqall)

Sequence database USA.

[-outfile] (outfile)

Output filename.

    [ Team LiB ] Previous Section Next Section