prophet

Function

Description

prophet scans one or more sequences with the supplied Gribskov or Henikoff profile (produced by prophecy). It writes the the highest-scoring (gapped) alignment for each sequence to file, which includes the score for the match. Gap opening coefficient and gap extension coefficient must be specified.

Algorithm

prophet does NOT use the 'simple frequency matrices' produced by prophecy. If you have a 'simple frequency matrix'you should use the program profit to scan sequences.

Usage

Command line arguments


Input file format

prophet reads a GRIBSKOV or a HENIKOFF profile produced by prophecy and searches one or more sequence USAs.

Output file format

Data files

None.

Notes

None.

References

None.

Warnings

prophet does NOT use the 'simple frequency matrices' produced by prophecy. If you have a 'simple frequency matrix'you should use the program profit to scan sequences.

Diagnostic Error Messages

If you see something like:

   EMBOSS An error in prophet.c at line 110:
Unrecognised profile/matrix file format

then you have probably given prophet a 'simple frequency matrix' and not a a GRIBSKOV or HENIKOFF profile produced by prophecy.

If you see something like:

Uncaught exception
Allocation Failed
raised at prophet.c:147

then you have run out of memory. This will occur if your sequence and/or your profile are large.

prophet allocates space based on the length of the sequence and the size of the profile (8 x seqlen x profilelen)

Exit status

It exits with a status of 0.

Known bugs

None.

Author(s)

History

Target users

Comments