1 hmmpfam - search one or more sequences against HMM database
3 Copyright (C) 1992-2003 HHMI/Washington University School of Medicine
4 Freely distributed under the GNU General Public License (GPL)
5 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
6 HMM file: testInput.hmm
7 Sequence file: testInput.fasta
8 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
14 Scores for sequence family classification (score includes all domains):
15 Model Description Score E-value N
16 -------- ----------- ----- ------- ---
20 Model Domain seq-f seq-t hmm-f hmm-t score E-value
21 -------- ------- ----- ----- ----- ----- ----- -------
22 TEST 1/1 8 97 .] 1 95 [] 184.7 2.5e-56
24 Alignments of top-scoring domains:
25 TEST: domain 1 of 1, from 8 to 97: score 184.7, E = 2.5e-56
26 *->svfqqqqssksttgstvtAiAiAigYRYRYRAvtWnsGsLssGvnDn
27 sv+qqqq+ + +vtAiAiAigYRYRYRAv Wn GsLs G nDn
28 Test 8 SVYQQQQGGSA----MVTAIAIAIGYRYRYRAVVWNKGSLSTGTNDN 50
30 DnDqqsdgLYtiYYsvtvpssslpsqtviHHHaHkasstkiiikiePr<-
31 DnDq +d LYtiYYsvtv +ss+p q+v+HHHaH+asstkiiiki P
32 Test 51 DNDQAAD-LYTIYYSVTVSASSWPGQSVTHHHAHPASSTKIIIKIAPS 97
41 Query sequence: lcl|gi|340783625|Plus1
43 Description: megaplasmid, complete sequence [UNKNOWN]
45 Scores for sequence family classification (score includes all domains):
46 Model Description Score E-value N
47 -------- ----------- ----- ------- ---
48 IS66_ORF3.uniq -30.0 1.7e-15 1
51 Model Domain seq-f seq-t hmm-f hmm-t score E-value
52 -------- ------- ----- ----- ----- ----- ----- -------
53 IS66_ORF3.uniq 1/1 60792 61184 .. 1 558 [] -30.0 1.7e-15
55 Alignments of top-scoring domains:
56 IS66_ORF3.uniq: domain 1 of 1, from 60792 to 61184: score -30.0, E = 1.7e-15
57 *->mSSPLDLSLFPNLmtevvkmnaalpddlslLkAalqkqQavvea..l
58 L S FP + + aa+p+ l + l + v + +
59 lcl|gi|340 60792 SCTGLWDSEFPPVFPVNGSQRAATPF-LPVGPDELSSPPSSVLXscY 60837
61 dakianleelIEkLka.qldkLRRmkFGkrSEKkRhkLeldiaQiekale
63 lcl|gi|340 60838 DFPLCQSLSLCFRLGFpQALALRVRQRAPRGCH----------------- 60870
65 dLEddLnetqaeiaeaekkieaPAssppkeaPiaskPesPRqerkRkPLP
66 + +a+ + A+ +p++ + + LP
67 lcl|gi|340 60871 ------AHRRARGLFFTDTPHPVAVPRPQR----------DLTGSLVTLP 60904
69 eeLPReerRlePestvCPcggGqLkriGEDvsEqLDlvpaaFeVIqTvRp
70 + +P e P + + P + + D vp+ +
71 lcl|gi|340 60905 APMPCSET---PVESQHPGH-----------TGRSDAVPTPNTMKTS--- 60937
73 KyACrqCdtiVQAPaPakpIErGiptaGLLArvlVSKyaEHlPLYRQsEI
75 lcl|gi|340 - -------------------------------------------------- -
77 yaRqGVeiaRstLadWVgrtgarLaPLvdALaeyVLkeGklHADeTPVqV
79 lcl|gi|340 60938 ------AIMISGLIHGVSARCLRF-------------------------- 60955
81 LaPGnkKTKTGyLWAYvRDDRnagsSlppaVvfayspdRkGEHaqvhLae
82 k TG R D Slp+ + + R G e
83 lcl|gi|340 60956 -----KRSVTG-----PRQD-----SLPADGL---RLCREG-------VE 60980
85 ysGkLqaDAyaGYnalyesgRVGGDpikEAgCwAHaRRKifDlhvrnpSe
86 sG + + + ++ wAH RR ++D+ pS+
87 lcl|gi|340 60981 PSGSQXKVSVTSILLSRTXP---------VASWAHVRRHFYDVNQGAPSP 61021
89 lteEALerIaaLYeIEaeIRGspaEdRlavRqelSvPlvkslseWLtaqr
90 ++ AL rI +LYeIEaeI p+E+ Rq+++ Pl+ s+ WL r
91 lcl|gi|340 61022 VAQTALLRIHELYEIEAEIKDDPPEQKMLARQQRAAPLLESFWTWLNDTR 61071
93 KtlsrgseLAKALaYlLnrWdALtryldDGqVpIDNNivENAiRrvAlGR
94 + +s+ AKA Y+L+rW ALt+yl+ G+ IDNN vE A+R vA+GR
95 lcl|gi|340 61072 AQVAPKSAIAKAIGYALKRWKALTLYLQEGHLGIDNNPVERALRGVAIGR 61121
97 kNyLFaGSDrGGeraAliySLIgTCKmNgVePyaYLrdVltrLadwplnk
98 kN+LFaG D GGeraA++y LI TCK+NgVeP aYL+dVl++L++wp +
99 lcl|gi|340 61122 KNFLFAGNDAGGERAAAFYTLIETCKLNGVEPFAYLCDVLEKLPTWPNKR 61171
103 lcl|gi|340 61172 -LHELLPWNWKNLT 61184