hmmsearch - search a sequence database with a profile HMM HMMER 2.1.1 (Dec 1998) Copyright (C) 1992-1998 Washington University School of Medicine HMMER is freely distributed under the GNU General Public License (GPL). - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - HMM file: sequences.prf [sequences] Sequence database: /usr/people/latek/smalldb.fasta - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query HMM: sequences|| [HMM has been calibrated; E-values are empirical estimates] Scores for complete sequences (score includes all domains): Sequence Description Score E-value N -------- ----------- ----- ------- --- gi|4557671|ref|NP_000198.1| (NM_000207) proinsulin 135.3 5.4e-40 1 gi|266377|sp|P30410|INS_PANTR INSULIN PRECURSOR >gi| 134.0 1.3e-39 1 gi|208668|gb|AAA72172.1| (J02547) synthetic pre 132.6 3.6e-39 1 gi|266373|sp|P30406|INS_MACFA INSULIN PRECURSOR >gi| 132.2 4.9e-39 1 gi|266362|sp|P30407|INS_CERAE INSULIN PRECURSOR >gi| 131.6 7e-39 1 gi|58103|emb|CAA23424.1| (V00082) reading frame 129.0 4.3e-38 1 gi|223965|prf||1006230A insulin,pro- [Macaca m 127.1 1.7e-37 1 gi|1352467|sp|P01311|INS_RABIT INSULIN PRECURSOR >gi| 124.0 1.4e-36 1 gi|124575|sp|P01321|INS_CANFA INSULIN PRECURSOR >gi| 121.2 9.5e-36 1 gi|307072|gb|AAA59179.1| (L15440) insulin [Homo 119.8 2.7e-35 1 Parsed for domains: Sequence Domain seq-f seq-t hmm-f hmm-t score E-value -------- ------- ----- ----- ----- ----- ----- ------- gi|4557671|ref|NP_000198.1| 1/1 22 76 .. 1 55 [] 135.3 5.4e-40 gi|266377|sp|P30410|INS_PANTR 1/1 22 76 .. 1 55 [] 134.0 1.3e-39 gi|208668|gb|AAA72172.1| 1/1 26 80 .. 1 55 [] 132.6 3.6e-39 gi|266373|sp|P30406|INS_MACFA 1/1 22 76 .. 1 55 [] 132.2 4.9e-39 gi|266362|sp|P30407|INS_CERAE 1/1 22 76 .. 1 55 [] 131.6 7e-39 gi|58103|emb|CAA23424.1| 1/1 1 53 [. 1 55 [] 129.0 4.3e-38 gi|223965|prf||1006230A 1/1 1 52 [. 1 55 [] 127.1 1.7e-37 gi|1352467|sp|P01311|INS_RABIT 1/1 22 76 .. 1 55 [] 124.0 1.4e-36 gi|124575|sp|P01321|INS_CANFA 1/1 22 76 .. 1 55 [] 121.2 9.5e-36 gi|307072|gb|AAA59179.1| 1/1 22 73 .. 1 55 [] 119.8 2.7e-35 Alignments of top-scoring domains: gi|4557671|ref|NP_000198.1|: domain 1 of 1, from 22 to 76: score 135.3, E = 5.4e-40 *->aaaFvNqhLaesklvEALalkCGERGFFYTPKTRREAEdlqVgqvel aaaFvNqhL++s+lvEAL+l+CGERGFFYTPKTRREAEdlqVgqvel gi|4557671 22 AAAFVNQHLCGSHLVEALYLVCGERGFFYTPKTRREAEDLQVGQVEL 68 gGgpaaie<-* gGgp+a+ gi|4557671 69 GGGPGAGS 76 gi|266377|sp|P30410|INS_PANTR: domain 1 of 1, from 22 to 76: score 134.0, E = 1.3e-39 *->aaaFvNqhLaesklvEALalkCGERGFFYTPKTRREAEdlqVgqvel a aFvNqhL++s+lvEAL+l+CGERGFFYTPKTRREAEdlqVgqvel gi|266377| 22 ASAFVNQHLCGSHLVEALYLVCGERGFFYTPKTRREAEDLQVGQVEL 68 gGgpaaie<-* gGgp+a+ gi|266377| 69 GGGPGAGS 76 gi|208668|gb|AAA72172.1|: domain 1 of 1, from 26 to 80: score 132.6, E = 3.6e-39 *->aaaFvNqhLaesklvEALalkCGERGFFYTPKTRREAEdlqVgqvel + FvNqhL++s+lvEAL+l+CGERGFFYTPKTRREAEdlqVgqvel gi|208668| 26 FRMFVNQHLCGSHLVEALYLVCGERGFFYTPKTRREAEDLQVGQVEL 72 gGgpaaie<-* gGgp+a+ gi|208668| 73 GGGPGAGS 80 gi|266373|sp|P30406|INS_MACFA: domain 1 of 1, from 22 to 76: score 132.2, E = 4.9e-39 *->aaaFvNqhLaesklvEALalkCGERGFFYTPKTRREAEdlqVgqvel a aFvNqhL++s+lvEAL+l+CGERGFFYTPKTRREAEd+qVgqvel gi|266373| 22 APAFVNQHLCGSHLVEALYLVCGERGFFYTPKTRREAEDPQVGQVEL 68 gGgpaaie<-* gGgp+a+ gi|266373| 69 GGGPGAGS 76 gi|266362|sp|P30407|INS_CERAE: domain 1 of 1, from 22 to 76: score 131.6, E = 7e-39 *->aaaFvNqhLaesklvEALalkCGERGFFYTPKTRREAEdlqVgqvel + aFvNqhL++s+lvEAL+l+CGERGFFYTPKTRREAEd+qVgqvel gi|266362| 22 VPAFVNQHLCGSHLVEALYLVCGERGFFYTPKTRREAEDPQVGQVEL 68 gGgpaaie<-* gGgp+a+ gi|266362| 69 GGGPGAGS 76 gi|58103|emb|CAA23424.1|: domain 1 of 1, from 1 to 53: score 129.0, E = 4.3e-38 *->aaaFvNqhLaesklvEALalkCGERGFFYTPKTRREAEdlqVgqvel FvNqhL++s+lvEAL+l+CGERGFFYTPKTRREAEdlqVgqvel gi|58103|e 1 --MFVNQHLCGSHLVEALYLVCGERGFFYTPKTRREAEDLQVGQVEL 45 gGgpaaie<-* gGgp+a+ gi|58103|e 46 GGGPGAGS 53 gi|223965|prf||1006230A: domain 1 of 1, from 1 to 52: score 127.1, E = 1.7e-37 *->aaaFvNqhLaesklvEALalkCGERGFFYTPKTRREAEdlqVgqvel FvNqhL++s+lvEAL+l+CGERGFFYTPKTRREAEd+qVgqvel gi|223965| 1 ---FVNQHLCGSHLVEALYLVCGERGFFYTPKTRREAEDPQVGQVEL 44 gGgpaaie<-* gGgp+a+ gi|223965| 45 GGGPGAGS 52 gi|1352467|sp|P01311|INS_RABIT: domain 1 of 1, from 22 to 76: score 124.0, E = 1.4e-36 *->aaaFvNqhLaesklvEALalkCGERGFFYTPKTRREAEdlqVgqvel a aFvNqhL++s+lvEAL+l+CGERGFFYTPK+RRE+E+lqVgq el gi|1352467 22 AQAFVNQHLCGSHLVEALYLVCGERGFFYTPKSRREVEELQVGQAEL 68 gGgpaaie<-* gGgp+a+ gi|1352467 69 GGGPGAGG 76 gi|124575|sp|P01321|INS_CANFA: domain 1 of 1, from 22 to 76: score 121.2, E = 9.5e-36 *->aaaFvNqhLaesklvEALalkCGERGFFYTPKTRREAEdlqVgqvel + aFvNqhL++s+lvEAL+l+CGERGFFYTPK+RRE+EdlqV++vel gi|124575| 22 TRAFVNQHLCGSHLVEALYLVCGERGFFYTPKARREVEDLQVRDVEL 68 gGgpaaie<-* +G p+++ gi|124575| 69 AGAPGEGG 76 gi|307072|gb|AAA59179.1|: domain 1 of 1, from 22 to 73: score 119.8, E = 2.7e-35 *->aaaFvNqhLaesklvEALalkCGERGFFYTPKTRREAEdlqVgqvel aaaFvNqhL++s+lvEAL+l+CGERGFFYTPKTRREAE Vgqvel gi|307072| 22 AAAFVNQHLCGSHLVEALYLVCGERGFFYTPKTRREAE---VGQVEL 65 gGgpaaie<-* gGgp+a+ gi|307072| 66 GGGPGAGS 73 Histogram of all scores: score obs exp (one = represents 1 sequences) ----- --- --- -37 1 2|=* -36 1 2|=* -35 1 3|= * -34 0 3| * -33 9 3|==*====== -32 5 2|=*=== -31 2 2|=* > -30 11 -|=========== % Statistical details of theoretical EVD fit: mu = -33.9379 lambda = 0.3059 chi-sq statistic = 0.0000 P(chi-square) = 0 Total sequences searched: 30 Whole sequence top hits: tophits_s report: Total hits: 10 Satisfying E cutoff: 10 Total memory: 18K Domain top hits: tophits_s report: Total hits: 10 Satisfying E cutoff: 10 Total memory: 21K