prophet.txt

emboss的linux版本的源代码

0.46 0.00 -0.89 2.13 1.78 -1.45 0.86 0.63 -0.33 0.00 0.67 -0.72 -0.50 0.93 0.00
 0.17 0.97 0.11 0.33 0.33 0.00 -0.33 -1.67 0.00 -0.85 0.00 4.55 4.55
-0.28 0.00 -0.20 -1.30 -0.80 1.86 -0.70 -0.17 0.67 0.00 -0.72 1.52 0.57 -0.60 0
.00 -0.70 -0.99 -0.52 -0.21 -0.33 0.00 0.19 1.45 0.00 1.66 0.00 4.55 4.55
0.79 0.00 0.39 0.59 0.52 -0.69 1.39 -0.24 -0.02 0.00 0.15 -0.46 -0.23 0.47 0.00
 0.52 0.02 -0.26 0.86 1.64 0.00 0.29 -1.13 0.00 -0.70 0.00 4.55 4.55
0.84 0.00 0.07 0.03 0.22 -0.34 0.36 0.17 -0.03 0.00 -0.05 0.01 -0.06 -0.01 0.00
 1.35 0.18 0.02 0.43 0.34 0.00 0.19 -0.78 0.00 -0.50 0.00 4.55 4.55
1.12 0.00 -0.39 1.32 1.48 -1.09 0.88 0.39 -0.24 0.00 0.32 -0.47 -0.30 0.62 0.00
 0.64 0.75 -0.09 0.47 0.45 0.00 -0.08 -1.63 0.00 -0.80 0.00 4.55 4.55
0.59 0.00 -0.23 -0.13 -0.06 0.22 0.13 -0.28 0.98 0.00 0.21 0.96 0.90 -0.16 0.00
 0.16 0.01 -0.21 -0.02 0.29 0.00 1.32 -0.60 0.00 -0.24 0.00 4.55 4.55
0.10 0.00 -0.78 0.58 0.62 -0.35 -0.04 1.19 -0.17 0.00 0.31 0.22 0.17 0.63 0.00
0.23 1.43 0.43 -0.25 -0.12 0.00 -0.10 -0.31 0.00 -0.22 0.00 4.55 4.55
1.60 0.00 0.24 0.38 0.40 -0.54 1.03 -0.27 0.07 0.00 0.06 0.05 0.31 0.26 0.00 0.
55 0.23 -0.29 0.68 0.53 0.00 0.41 -1.10 0.00 -0.55 0.00 4.55 4.55
1.41 0.00 0.80 0.40 0.40 -0.64 1.00 -0.27 -0.10 0.00 0.20 -0.47 -0.30 0.43 0.00
 0.74 0.03 -0.10 1.77 0.57 0.00 0.03 -0.24 0.00 -0.60 0.00 4.55 4.55
-0.29 0.00 -0.52 -0.86 -0.64 1.61 -0.82 -0.18 0.74 0.00 -0.25 1.68 1.44 -0.47 0
.00 -0.72 -0.37 -0.13 -0.51 -0.26 0.00 0.61 0.72 0.00 0.86 0.00 4.55 4.55
0.35 0.00 -0.50 1.14 0.87 -0.84 0.47 0.50 0.16 0.00 0.41 -0.19 -0.11 0.89 0.00
0.16 0.94 0.09 0.17 0.40 0.00 0.08 -1.08 0.00 -0.54 0.00 4.55 4.55
0.09 0.00 -0.26 0.19 0.16 -0.68 -0.11 0.05 0.03 0.00 1.63 -0.49 0.15 0.36 0.00
0.12 0.27 0.85 0.41 0.33 0.00 -0.03 -0.31 0.00 -0.42 0.00 4.55 4.55
0.55 0.00 0.24 -0.49 -0.29 0.76 -0.02 -0.31 1.20 0.00 -0.43 1.11 0.67 -0.36 0.0
0 -0.07 -0.47 -0.60 -0.03 0.19 0.00 1.14 -0.28 0.00 0.48 0.00 4.55 4.55
-0.07 0.00 -0.67 -0.71 -0.49 1.61 -0.56 -0.30 1.36 0.00 -0.50 2.05 1.57 -0.58 0
.00 -0.41 -0.28 -0.66 -0.54 -0.07 0.00 1.50 0.38 0.00 0.51 0.00 4.55 4.55
1.16 0.00 -0.16 0.64 0.79 -0.38 0.64 0.01 0.12 0.00 0.11 0.11 0.16 0.34 0.00 0.
38 0.35 -0.28 0.63 0.40 0.00 0.24 -0.79 0.00 -0.44 0.00 4.55 4.55
0.36 0.00 -0.06 0.32 0.42 -0.24 0.73 -0.04 -0.03 0.00 0.24 -0.07 0.16 0.28 0.00
 0.32 0.21 0.24 0.90 0.28 0.00 0.18 0.13 0.00 -0.60 0.00 4.55 4.55
0.12 0.00 0.21 -0.52 -0.36 0.75 0.09 -0.39 1.27 0.00 -0.36 1.21 0.79 -0.40 0.00
 -0.03 -0.49 -0.46 0.20 0.16 0.00 1.50 -0.16 0.00 0.26 0.00 4.55 4.55
1.38 0.00 0.45 0.35 0.35 -0.74 0.71 -0.02 -0.15 0.00 0.36 -0.41 -0.08 0.35 0.00
 0.72 0.25 0.28 1.10 0.46 0.00 0.02 -0.11 0.00 -0.65 0.00 4.55 4.55
0.19 0.00 0.10 0.27 0.22 -0.07 0.10 0.24 0.15 0.00 0.55 -0.11 0.00 0.78 0.00 -0
.09 0.03 0.03 0.20 0.82 0.00 0.03 -0.12 0.00 0.08 0.00 4.55 4.55
0.85 0.00 0.09 0.49 0.33 -0.41 0.57 -0.10 0.63 0.00 -0.03 0.30 0.25 0.09 0.00 0
.30 0.15 -0.31 0.17 0.36 0.00 0.96 -1.20 0.00 -0.35 0.00 0.92 0.92
-0.09 0.00 -0.63 -0.54 -0.36 1.36 -0.59 -0.32 1.43 0.00 -0.36 1.73 1.36 -0.50 0
.00 -0.36 -0.23 -0.50 -0.41 0.00 0.00 1.24 0.21 0.00 0.32 0.00 0.92 0.92
1.31 0.00 0.23 0.36 0.36 -0.60 0.63 -0.11 0.07 0.00 0.31 -0.17 0.03 0.30 0.00 0
.53 0.16 -0.16 0.47 1.05 0.00 0.22 -0.87 0.00 -0.45 0.00 0.92 0.92
0.65 0.00 0.26 0.40 0.40 -0.44 0.39 0.51 -0.20 0.00 0.35 -0.35 -0.24 0.48 0.00
0.46 0.32 0.27 0.81 0.26 0.00 -0.14 -0.14 0.00 -0.25 0.00 0.92 0.92
0.04 0.00 -0.82 0.51 0.58 -0.82 -0.05 0.18 -0.27 0.00 1.89 -0.41 0.22 0.56 0.00
 0.14 0.58 0.98 0.27 0.27 0.00 -0.27 -0.03 0.00 -0.81 0.00 0.92 0.92
-0.37 0.00 1.14 -0.67 -0.64 1.67 -0.78 0.33 0.21 0.00 -0.71 0.54 0.09 -0.17 0.0
0 -1.01 -0.74 -0.71 -0.53 -0.37 0.00 -0.04 1.31 0.00 1.83 0.00 0.92 0.92
-0.16 0.00 -0.43 0.35 0.34 -0.55 -0.20 0.97 -0.38 0.00 0.98 -0.43 0.00 0.60 0.0
0 0.26 0.63 1.21 0.07 0.00 0.00 -0.38 0.58 0.00 -0.38 0.00 0.92 0.92
0.14 0.00 -0.27 0.54 0.64 -0.34 0.25 0.19 -0.09 0.00 0.14 -0.17 -0.12 0.25 0.00
 0.05 0.28 0.00 0.09 0.09 0.00 -0.09 -0.52 0.00 -0.24 0.00 0.92 0.92
0.18 0.00 -0.22 -0.12 -0.06 0.32 -0.08 -0.08 0.26 0.00 -0.10 0.48 0.39 -0.10 0.
00 -0.03 0.00 -0.17 -0.07 0.02 0.00 0.29 0.05 0.00 0.06 0.00 0.92 0.92
0.41 0.00 0.11 0.24 0.21 -0.27 0.58 -0.08 -0.10 0.00 -0.03 -0.18 -0.10 0.16 0.0
0 0.17 0.09 -0.14 0.25 0.19 0.00 0.09 -0.44 0.00 -0.24 0.00 0.92 0.92
-0.10 0.00 0.33 -0.16 -0.16 0.43 -0.20 0.10 0.03 0.00 -0.20 0.10 -0.03 -0.03 0.
00 -0.26 -0.20 -0.20 -0.13 -0.10 0.00 -0.03 0.36 0.00 0.49 0.00 0.92 0.92
0.07 0.00 -0.20 0.20 0.20 -0.26 0.07 0.20 -0.10 0.00 0.13 -0.03 0.00 0.13 0.00
0.10 0.49 0.13 -0.03 -0.03 0.00 -0.07 -0.16 0.00 -0.20 0.00 0.92 0.92
0.20 0.00 0.07 0.20 0.16 -0.20 0.49 -0.07 -0.10 0.00 -0.03 -0.16 -0.10 0.13 0.0
0 0.10 0.07 -0.10 0.20 0.13 0.00 0.07 -0.33 0.00 -0.20 0.00 0.92 0.92

Output file format

  Output files for usage example

  File: hd_fugru.prophet

########################################
# Program: prophet
# Rundate: Sat 15 Jul 2006 12:00:00
# Commandline: prophet
#    -sequence "tsw:h*"
#    -infile ../../data/globins.gribskov
# Align_format: simple
# Report_file: hd_fugru.prophet
########################################

#=======================================
#
# Aligned_sequences: 2
# 1: globins
# 2: HD_FUGRU
# Matrix: EBLOSUM62
#
# Length: 167
# Identity:      20/167 (12.0%)
# Similarity:    46/167 (27.5%)
# Gaps:           6/167 ( 3.6%)
# Score: 74.34
#
#
#=======================================

globins            3 VDTGSV-VSLS-EEELSAVDKAWVKANSVAEVGGHALERGLFASEPMTLE     50
                     :....| ..|| .|..:.........:.:......:....|...|....:
HD_FUGRU        2289 LQQPCVWNKLSTPEYNTHTCSLIYCLHHIILAVAVSPGDQLLHPERKKTK   2338

globins           51 FFDTFKYLSTFDLSKGSADVKAHGKKVLDALGDAVAHLDDLEGT-LAALS     99
                     ...........|....:..::....:::..|.:.:..:..|... ..|..
HD_FUGRU        2339 ALRHSDDEDQVDSVHDNHTLEWQACEIMAELVEGLQSVLSLGHHRNTAFP   2388

globins          100 DLHAHKLKKGVDPVNFKLLSHCLLVVLASHLPGDFTPEVQASMDKFLASV    149
                     ......|:..:..::...|.:....|........::|:........|..:
HD_FUGRU        2389 AFLTPTLRNIIISLSRLPLVNSHTRVPPLVWKLGWSPQPGGEFGTTLPEI   2438

globins          150 ATV-LASK--YRELGYQ    163
                     ... |..|  :||..|:
HD_FUGRU        2439 PVDFLQEKDVFREFLYR   2455


#=======================================
#
# Aligned_sequences: 2
# 1: globins
# 2: HIRA_FUGRU
# Matrix: EBLOSUM62
#


  [Part of this file has been deleted for brevity]

#
#=======================================

globins           15 ELSAVDKAWVKAN-SVAEVGGHALERGLFASEPMTLEFFDTFKYLSTFDL     63
                     .||..||..|||. . .:||.||.|.|..|.|.|.|.|..|..|...|||
HBA_HUMAN          1 VLSPADKTNVKAAWG-GKVGAHAGEYGAEALERMFLSFPTTKTYFPHFDL     49

globins           64 SKGSADVKAHGKKVLDALGDAVAHLDDLEGTLAALSDLHAHKLK-KGVDP    112
                     |.|||.||.|||||.|||.:||||:||:...|:||||||||||: .  ||
HBA_HUMAN         50 SHGSAQVKGHGKKVADALTNAVAHVDDMPNALSALSDLHAHKLRVD--DP     97

globins          113 VNFKLLSHCLLVVLASHLPGDFTPEVQASMDKFLASVATVLASKYR    158
                     ||||||||||||.||:|||.:|||.|.||:|||||||:|||.||||
HBA_HUMAN         98 VNFKLLSHCLLVTLAAHLPAEFTPAVHASLDKFLASVSTVLTSKYR    143


#=======================================
#
# Aligned_sequences: 2
# 1: globins
# 2: HBB_HUMAN
# Matrix: EBLOSUM62
#
# Length: 153
# Identity:      77/153 (50.3%)
# Similarity:    91/153 (59.5%)
# Gaps:           7/153 ( 4.6%)
# Score: 178.29
#
#
#=======================================

globins            9 VSLSEEELSAVDKAWVKAN-SVAEVGGHALE-RGLFASEPMTLEFFDTFK     56
                     |.|:.||.|||...|.|.| . .||||.||. . |....|.|..||::|.
HBB_HUMAN          1 VHLTPEEKSAVTALWGKVNVD-DEVGGEALGRL-LLVVYPWTQRFFESFG     48

globins           57 YLSTFDLSKGSADVKAHGKKVLDALGDAVAHLDDLEGTLAALSDLHAHKL    106
                     .|||.|...|:..|||||||||.|..|.:||||:|:||.|.||:||..||
HBB_HUMAN         49 DLSTPDAVMGNPKVKAHGKKVLGAFSDGLAHLDNLKGTFATLSELHCDKL     98

globins          107 K-KGVDPVNFKLLSHCLLVVLASHLPGDFTPEVQASMDKFLASVATVLAS    155
                     . .  ||.||:||.:.|:.|||.|...:|||.|||:..|.:|.||..||.
HBB_HUMAN         99 HVD--DPENFRLLGNVLVCVLAHHFGKEFTPPVQAAYQKVVAGVANALAH    146

globins          156 KYR    158
                     ||.
HBB_HUMAN        147 KYH    149


#---------------------------------------
#---------------------------------------

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.

See also

   Program name                     Description
   profit       Scan a sequence or database with a matrix or profile
   prophecy     Creates matrices/profiles from multiple alignments

Author(s)

   Alan Bleasby (ajb