# 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/TRM13.hmm.txt # target sequence database: proteomes/Trypanosoma_brucei.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: TRM13 [M=262] Accession: PF05206.13 Description: Methyltransferase TRM13 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2.1e-37 128.0 5.7 1.4e-27 95.8 3.2 3.3 2 tr|Q57Y07|Q57Y07_TRYB2 Uncharacterized protein OS=Trypanosom 5.1e-16 57.9 0.2 8.1e-16 57.3 0.2 1.3 1 tr|Q385P4|Q385P4_TRYB2 Uncharacterized protein OS=Trypanosom Domain annotation for each sequence (and alignments): >> tr|Q57Y07|Q57Y07_TRYB2 Uncharacterized protein OS=Trypanosoma brucei brucei (strain 927/4 GUTat10.1) OX=185431 GN=Tb # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 95.8 3.2 3.2e-31 1.4e-27 21 179 .. 245 468 .. 232 474 .. 0.89 2 ! 30.3 0.1 3.2e-11 1.4e-07 214 261 .. 471 519 .. 469 520 .. 0.92 Alignments for each domain: == domain 1 score: 95.8 bits; conditional E-value: 3.2e-31 TRM13 21 ayvEfGaGkgelsryvnqaleed............................kkenskyvLidrasnrlkaDrklrkdeselaikRlk 79 ++E+GaGk+ l+ +++al + + ++ v +d + r kaD ++r+ l + Rl+ tr|Q57Y07|Q57Y07_TRYB2 245 GFLELGAGKAGLAVSLQDALLSGsirychshpgfqerkegrtavnealsadNGGFPRIVVVDMDNFRRKADACVRNSR--LPLVRLR 329 689***************99888999**********************9988999*********************99..699**** PP TRM13 80 idIkdLnlsaleelekkkkv............................vavsKHLCGaATDLtLrcllnsekask.kakleglliAl 137 i+IkdL+l + + + + vav+KHLCGa TD++L c+++ + +++ a+++++++A+ tr|Q57Y07|Q57Y07_TRYB2 330 INIKDLDLAKALCGA--STCkrsreetsstsdqctnaikqsvsveeqwVAVGKHLCGACTDFALSCVTAPNLNTEgFASVCAVVFAT 414 ********9987633..2223345577999**************************************998887756899******* PP TRM13 138 CChhvCswkeyvn............key.leelgitkeefqiltkmvSWavsgkr 179 CCh++C+ k+ + +++ ++++ef+++t+m+SWavsg+ tr|Q57Y07|Q57Y07_TRYB2 415 CCHQRCELKHINAllspsecrrgriA-IpGTNHTFSEREFAAITSMSSWAVSGEA 468 *********9888999****988642.2357899******************987 PP == domain 2 score: 30.3 bits; conditional E-value: 3.2e-11 TRM13 214 eereeiGlkakrlidegRllalkekgfe.aelvkYvekevslEnvlLla 261 +er+ G+ +kr+id+ Rl++lk+ gf + l +Y e++v+ En++++a tr|Q57Y07|Q57Y07_TRYB2 471 AERRLTGMCCKRVIDAFRLEYLKQSGFRrTYLCQYAESSVTGENICIVA 519 589999*******************9963789**************987 PP >> tr|Q385P4|Q385P4_TRYB2 Uncharacterized protein OS=Trypanosoma brucei brucei (strain 927/4 GUTat10.1) OX=185431 GN=Tb # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 57.3 0.2 1.9e-19 8.1e-16 73 261 .. 297 484 .. 225 485 .. 0.80 Alignments for each domain: == domain 1 score: 57.3 bits; conditional E-value: 1.9e-19 TRM13 73 laikRlkidIkdLnlsaleelekkkkvvavsKHLCGaATDLtLrcllnsekaskkakleglliAlCChhvCswkeyvnkeyleel.g 158 ++R+ ++d + + + + +++ +v+KHLCG+ D Lr l++s+k + l++ CC ++ ey n +l++ + tr|Q385P4|Q385P4_TRYB2 297 YPLHRISKRLQDVDWGNTHI--DFDSCAVVAKHLCGTSVDSCLRHLMASNKL-----PRVLVVVPCCFNKGRYAEYCNPGFLSHVaD 376 55666666666666666666..6788999*******************9997.....589*********************887658 PP TRM13 159 itke.efqiltkmvSWavsgkrkse........aeeedekeekekeekkeeeeseeekksklsseereeiGlkakrlidegRllalk 236 ++ e +++ t+++ W s ++ ++++ + +++ ekk++++ + e++ l + i + +i++gR+l+l+ tr|Q385P4|Q385P4_TRYB2 377 VEGElSWEQCTRLTDWNKSCYQQAGnplscsscCSCCNGSYDRRAGEKKVSRKRKMEHF--LPCM--DDIATLVEVIINYGRVLWLR 459 99999***********9776664446999998866677777777778877777777776..7666..589***************** PP TRM13 237 ekgfeaelvkYvekevslEnvlLla 261 e+g+++++v+Yv++ v++ n +++a tr|Q385P4|Q385P4_TRYB2 460 EMGYQTSVVEYVPRCVTPKNRAIVA 484 *********************9987 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (262 nodes) Target sequences: 8561 (4332363 residues searched) Passed MSV filter: 498 (0.0581708); expected 171.2 (0.02) Passed bias filter: 265 (0.0309543); expected 171.2 (0.02) Passed Vit filter: 17 (0.00198575); expected 8.6 (0.001) Passed Fwd filter: 2 (0.000233618); expected 0.1 (1e-05) Initial search space (Z): 8561 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.10u 0.00s 00:00:00.10 Elapsed: 00:00:00.04 # Mc/sec: 28376.98 // [ok]