# 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/TruB_N.hmm.txt # target sequence database: proteomes/Trichomonas_vaginalis.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: TruB_N [M=149] Accession: PF01509.17 Description: TruB family pseudouridylate synthase (N terminal domain) 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-22 80.4 3.0 2.3e-13 52.2 0.0 4.0 3 tr|A2F410|A2F410_TRIVA Pseudouridine syntase related protein ------ inclusion threshold ------ 0.42 12.4 0.1 0.42 12.4 0.1 3.7 4 tr|A2DJE2|A2DJE2_TRIVA Uncharacterized protein OS=Trichomona 1.2 10.9 2.3 9.7 8.0 0.0 2.2 2 tr|A2FHQ9|A2FHQ9_TRIVA Uncharacterized protein OS=Trichomona Domain annotation for each sequence (and alignments): >> tr|A2F410|A2F410_TRIVA Pseudouridine syntase related protein, putative OS=Trichomonas vaginalis OX=5722 GN=TVAG_4071 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 52.2 0.0 1.4e-17 2.3e-13 1 53 [. 77 129 .. 77 132 .. 0.96 2 ! 22.5 0.1 1.9e-08 0.00032 72 149 .] 130 192 .. 127 192 .. 0.90 3 ? 1.3 0.1 0.068 1.1e+03 102 149 .] 296 343 .. 291 343 .. 0.90 Alignments for each domain: == domain 1 score: 52.2 bits; conditional E-value: 1.4e-17 HHHTT-S-EEESS---TT-EEEEEEEEGGGGGGHHHCTTS-EEEEEEEEETEE CS TruB_N 1 krllkakkvGhtGtLDplatGvLvvavgeaTklleylleadKeYeatirlgae 53 kr+l++ k Gh+GtLDp++tG+L++++++aT + ++ ++a KeY+a++rl ++ tr|A2F410|A2F410_TRIVA 77 KRILGVAKTGHSGTLDPKVTGCLIICIENATCIGKSQQNAGKEYVAVLRLFET 129 89**********************************************99765 PP == domain 2 score: 22.5 bits; conditional E-value: 1.9e-08 --HHHHHHHHHHC-EEE---------------------------------EEEEEEEEEEEEEE..EEEEEEEESTTTTH CS TruB_N 72 lteekleevlkkftGeieqvppmySAvkvnGkrlyelaregeeverkkrkvtiyelelleveep..evelevecskGtYi 149 + +eklee++k +tG+++q pp++SAvk + r ++y+ el+e++++ + ++v+c++GtYi tr|A2F410|A2F410_TRIVA 130 VPREKLEEAIKFLTGKLYQCPPAISAVK-----------------KQLRVRQVYSNELIEYDDEqkLAVMRVACEAGTYI 192 6789************************.................67899999*******9877668999*********8 PP == domain 3 score: 1.3 bits; conditional E-value: 0.068 --------------------EEEEEEEEEEEEEEEEEEEEEESTTTTH CS TruB_N 102 GkrlyelaregeeverkkrkvtiyelelleveepevelevecskGtYi 149 G ++ ++ +ge++ ++ ++t +++++ + + +v + +GtY+ tr|A2F410|A2F410_TRIVA 296 GDQVVVITAKGEAIALAIATMTAQQMSVNQKDVAATLKRVIMDRGTYV 343 667777889999****************9999999999*********7 PP >> tr|A2DJE2|A2DJE2_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_136600 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -1.5 0.1 0.49 8.1e+03 73 108 .. 66 101 .. 45 117 .. 0.71 2 ? -1.2 2.4 0.39 6.5e+03 30 107 .. 180 257 .. 173 268 .. 0.69 3 ? -2.0 0.1 0.71 1.2e+04 66 92 .. 423 448 .. 413 453 .. 0.72 4 ? 12.4 0.1 2.5e-05 0.42 76 136 .. 493 553 .. 461 562 .. 0.84 Alignments for each domain: == domain 1 score: -1.5 bits; conditional E-value: 0.49 -HHHHHHHHHHC-EEE-------------------- CS TruB_N 73 teekleevlkkftGeieqvppmySAvkvnGkrlyel 108 ++ kl++ ++f+++i+ + + S + + ++ly + tr|A2DJE2|A2DJE2_TRIVA 66 NKIKLDNKTERFQSKIKSIEIKLSHLTEEHQNLYFY 101 445567777888888888877777777777777655 PP == domain 2 score: -1.2 bits; conditional E-value: 0.39 GGGGHHHCTTS-EEEEEEEEETEEETTSBTTS-EEE-.--TT--HHHHHHHHHHC-EEE------------------- CS TruB_N 30 aTklleylleadKeYeatirlgaetdteDaegeileeekveklteekleevlkkftGeieqvppmySAvkvnGkrlye 107 Tkl++ ++ K +++++ +t+++ ++ +i +e + +++++++ +k++ e++ + m+ + k + +r y+ tr|A2DJE2|A2DJE2_TRIVA 180 DTKLMNEKIKELKDENTQLQKELRTKSQQKDIDIRTELADLSKKKKEVSNEMKDLIEEYKLQIKMFENEKNDRERRYQ 257 578888888888877777777777777777777766633444466667777777777777777777666666666665 PP == domain 3 score: -2.0 bits; conditional E-value: 0.71 -.--TT--HHHHHHHHHHC-EEE---- CS TruB_N 66 eekveklteekleevlkkftGeieqvp 92 e +ek++++++e++ ++f+ ++++ tr|A2DJE2|A2DJE2_TRIVA 423 E-MTEKIKKQRVENINNAFKAPQSKQK 448 3.3788999999999999988776655 PP == domain 4 score: 12.4 bits; conditional E-value: 2.5e-05 HHHHHHHHC-EEE---------------------------------EEEEEEEEEEEEEEE CS TruB_N 76 kleevlkkftGeieqvppmySAvkvnGkrlyelaregeeverkkrkvtiyelelleveepe 136 + ++ l+kf e++ ++ A+k +G + y+l+ +++e++ k+ + ++ + e++e+++++ tr|A2DJE2|A2DJE2_TRIVA 493 EAKDGLDKFLQELDAEREKVNAMKSEGSEYYNLVEKNKELKEKSEHFKVMNEEIIEMRKRK 553 45666777888888888889999*********************************99875 PP >> tr|A2FHQ9|A2FHQ9_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_005740 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 8.0 0.0 0.00058 9.7 71 115 .. 24 70 .. 9 92 .. 0.73 2 ? 1.7 1.1 0.051 8.5e+02 27 80 .. 107 161 .. 106 171 .. 0.75 Alignments for each domain: == domain 1 score: 8.0 bits; conditional E-value: 0.00058 T--HHHHHHHH..HHC-EEE--------------------------- CS TruB_N 71 klteekleevl..kkftGeieqvppmySAvkvnGkrlyelaregeev 115 +++ +++++++ ++ft +++ + + S+ +nGk+ly+ +r+ + + tr|A2FHQ9|A2FHQ9_TRIVA 24 TFDIRDIKTIMkySRFTADQYVTLLKQSSSTINGKELYTCTRKANVT 70 55566666666225789999999999***************996433 PP == domain 2 score: 1.7 bits; conditional E-value: 0.051 EGGGGGGHHHCTTS-EEEEEEEEETEEETTSBTTS-EEE-..--TT--HHHHHHH CS TruB_N 27 vgeaTklleylleadKeYeatirlgaetdteDaegeilee.ekveklteekleev 80 +++ T +e+++++ K ++ + ++ e +e++++e++++ ++ ++ +++++ + tr|A2FHQ9|A2FHQ9_TRIVA 107 INDSTNKTEIFQDKCKVFQHENEIENENVSENSTKETTTNvTNENHNHSQDILTQ 161 677889999999999**********999999999999877544444444444444 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (149 nodes) Target sequences: 50190 (17055200 residues searched) Passed MSV filter: 1824 (0.0363419); expected 1003.8 (0.02) Passed bias filter: 1031 (0.0205419); expected 1003.8 (0.02) Passed Vit filter: 61 (0.00121538); 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.25u 0.01s 00:00:00.26 Elapsed: 00:00:00.11 # Mc/sec: 23102.04 // [ok]