# hmmsearch :: search profile(s) against a sequence database # HMMER 3.1b2 (February 2015); http://hmmer.org/ # Copyright (C) 2015 Howard Hughes Medical Institute. # Freely distributed under the GNU General Public License (GPLv3). # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - # query HMM file: pfam_list/TrmO.hmm.txt # target sequence database: proteomes/Trichomonas_vaginalis.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: TrmO [M=119] Accession: PF01980.16 Description: tRNA-methyltransferase O Scores for complete sequences (score includes all domains): --- full sequence --- --- best 1 domain --- -#dom- E-value score bias E-value score bias exp N Sequence Description ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 4.8e-34 118.6 0.0 5.8e-34 118.3 0.0 1.0 1 tr|A2FP34|A2FP34_TRIVA Uncharacterized protein OS=Trichomona 1.3e-25 91.3 0.0 1.5e-25 91.1 0.0 1.1 1 tr|A2FC71|A2FC71_TRIVA Uncharacterized protein OS=Trichomona 2.1e-25 90.7 0.0 2.4e-25 90.5 0.0 1.1 1 tr|A2EBW9|A2EBW9_TRIVA Uncharacterized protein OS=Trichomona Domain annotation for each sequence (and alignments): >> tr|A2FP34|A2FP34_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_181300 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 118.3 0.0 3.5e-38 5.8e-34 2 119 .] 21 133 .. 18 133 .. 0.92 Alignments for each domain: == domain 1 score: 118.3 bits; conditional E-value: 3.5e-38 S-SSGGG--.-EEEEEE-HHH..HGGGTTGGG-SEEEEEEE-TTS...-.SS--E-TT-.SS-E.EGGGS--S-SSS-EEEEEEEEE CS TrmO 2 iPrqpglaeeaeaeielepeynseealegleefshlwllfvfhenaekklkakvrpprlggnervGvfatRsphRpnpiglsvvkle 88 +Pr ++ ++ e+++++ pey +++l gle+fsh++++ ++h++ +++ ++vrp++ ++ +GvfatRsp Rpnpig+++vk++ tr|A2FP34|A2FP34_TRIVA 21 VPRCAKDCT-EESTVTFYPEY--SDCLVGLENFSHVFVFAWMHQS-DRD-LQLVRPHFGAAKKDIGVFATRSPARPNPIGMTLVKVN 102 577777666.799********..********************99.344.699999999999************************* PP EEETTEEEEES----TT-EEEEEEE--HHHH CS TrmO 89 kvegnvlevsgvDlldgtpvlDiKPYvpyaD 119 k+eg +l+v+g+Dlldgtp++DiK Y +++D tr|A2FP34|A2FP34_TRIVA 103 KIEGITLYVTGLDLLDGTPIIDIKKYDTTLD 133 **************************88777 PP >> tr|A2FC71|A2FC71_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_335440 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 91.1 0.0 9.2e-30 1.5e-25 2 119 .] 25 137 .. 24 137 .. 0.89 Alignments for each domain: == domain 1 score: 91.1 bits; conditional E-value: 9.2e-30 S-SSGGG--.-EEEEEE-HHH..HGGGTTGGG-SEEEEEEE-TTS...-.SS--E-TT-.SS-E.EGGGS--S-SSS-EEEEEEEEE CS TrmO 2 iPrqpglaeeaeaeielepeynseealegleefshlwllfvfhenaekklkakvrpprlggnervGvfatRsphRpnpiglsvvkle 88 +P+ ++ e e+ +e++pe+ +++l+ lee + + l++++++++++ +++p+r n +Gvf tRsp Rpnpi+ s+v++ tr|A2FC71|A2FC71_TRIVA 25 TPKCGNGDE-KEGIVEIFPEF--TDCLKCLEEEKFIHLICWLDKAQRD--VLSCHPRRDPANPITGVFGTRSPSRPNPISCSLVEVL 106 566666555.89*********..********************99333..466778888889999********************** PP EEETTEEEEES----TT-EEEEEEE--HHHH CS TrmO 89 kvegnvlevsgvDlldgtpvlDiKPYvpyaD 119 +++ n ++v+g+Dl+dgtp+lDiKP++ +D tr|A2FC71|A2FC71_TRIVA 107 EISANNIKVKGLDLFDGTPLLDIKPFYFPFD 137 *************************986665 PP >> tr|A2EBW9|A2EBW9_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_259180 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 90.5 0.0 1.4e-29 2.4e-25 13 118 .. 31 132 .. 20 133 .. 0.85 Alignments for each domain: == domain 1 score: 90.5 bits; conditional E-value: 1.4e-29 EEEEEE-HHH..HGGGTTGGG-SEEEEEEE-TTS...-.SS--E-TT-.SS-E.EGGGS--S-SSS-EEEEEEEEEEEETTEEEEES CS TrmO 13 eaeielepeynseealegleefshlwllfvfhenaekklkakvrpprlggnervGvfatRsphRpnpiglsvvklekvegnvlevsg 99 + ie+ p++ +e++ gle+++++w++ +h + +++ ++v+p + ++ + Gvfat sp Rpnpig s ++++k++g +++v+g tr|A2EBW9|A2EBW9_TRIVA 31 GGVIEILPQF--SECILGLEKVKYIWVFSFLHLA-DRS-IQQVQPVHRRHLHPRGVFATCSPVRPNPIGKSFMEIQKIDGLKIYVKG 113 5899******..9*******************88.344.455665444444455********************************* PP ----TT-EEEEEEE--HHH CS TrmO 100 vDlldgtpvlDiKPYvpya 118 +Dl+dgtp+lDiK Y + + tr|A2EBW9|A2EBW9_TRIVA 114 LDLVDGTPILDIKLYDEGS 132 *************997665 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (119 nodes) Target sequences: 50190 (17055200 residues searched) Passed MSV filter: 517 (0.0103009); expected 1003.8 (0.02) Passed bias filter: 403 (0.00802949); expected 1003.8 (0.02) Passed Vit filter: 23 (0.000458259); expected 50.2 (0.001) Passed Fwd filter: 3 (5.97729e-05); expected 0.5 (1e-05) Initial search space (Z): 50190 [actual number of targets] Domain search space (domZ): 3 [number of targets reported over threshold] # CPU time: 0.17u 0.01s 00:00:00.18 Elapsed: 00:00:00.10 # Mc/sec: 20295.69 // [ok]