# 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/Drosophila_melanogaster.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2.8e-22 79.3 0.2 2.5e-19 69.7 0.2 3.4 1 sp|O44081|DKC1_DROME H/ACA ribonucleoprotein complex subun 2.6e-09 37.2 0.0 4.7e-09 36.4 0.0 1.3 1 tr|Q29QN4|Q29QN4_DROME IP08469p OS=Drosophila melanogaster O ------ inclusion threshold ------ 0.033 14.2 0.9 1 9.4 0.0 2.3 2 tr|Q9V3A8|Q9V3A8_DROME Ergic53, isoform A OS=Drosophila mela 0.13 12.3 0.0 0.16 11.9 0.0 1.1 1 tr|A1ZAJ6|A1ZAJ6_DROME GEO02592p1 OS=Drosophila melanogaster Domain annotation for each sequence (and alignments): >> sp|O44081|DKC1_DROME H/ACA ribonucleoprotein complex subunit 4 OS=Drosophila melanogaster OX=7227 GN=Nop60B PE=1 SV= # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 69.7 0.2 7.4e-23 2.5e-19 1 149 [] 108 224 .. 108 224 .. 0.88 Alignments for each domain: == domain 1 score: 69.7 bits; conditional E-value: 7.4e-23 HHHTT-S-EEESS---TT-EEEEEEEEGGGGGGHHHCTTS-EEEEEEEEETEEETTSBTTS-EEE-.--TT--HHHHHHHHHHC-EEE- CS TruB_N 1 krllkakkvGhtGtLDplatGvLvvavgeaTklleylleadKeYeatirlgaetdteDaegeileeekveklteekleevlkkftGeie 89 k++lk++k Gh+GtLDp++tG+L+v++++aT+l+++ ++a KeY+a ++l+ + ++ k+++ l+k++G++ sp|O44081|DKC1_DROME 108 KKILKVEKTGHSGTLDPKVTGCLIVCIDRATRLVKSQQSAGKEYVAIFKLHGAVESVA-----------------KVRQGLEKLRGALF 179 89************************************************98877653.................68999********* PP --------------------------------EEEEEEEEEEEEEE..EEEEEEEESTTTTH CS TruB_N 90 qvppmySAvkvnGkrlyelaregeeverkkrkvtiyelelleveep..evelevecskGtYi 149 q+pp +SAvk r+ r t+y+ +ll+++e + v+c++G+Yi sp|O44081|DKC1_DROME 180 QRPPLISAVK-----------------RQLRVRTVYDSKLLDYDETrnMGVFWVSCEAGSYI 224 **********.................77788888888888887665566778888888887 PP >> tr|Q29QN4|Q29QN4_DROME IP08469p OS=Drosophila melanogaster OX=7227 GN=Dmel\CG7849 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 36.4 0.0 1.4e-12 4.7e-09 18 147 .. 109 240 .. 103 241 .. 0.93 Alignments for each domain: == domain 1 score: 36.4 bits; conditional E-value: 1.4e-12 T-EEEEEEEEGGGGGGHHHCTTS-..EEEEEEEEETEEETTSBTTS-EEE-.--TT--HHHHHHHHHHC-EEE-------------- CS TruB_N 18 latGvLvvavgeaTklleyllead..KeYeatirlgaetdteDaegeileeekveklteekleevlkkftGeieqvppmySAvkvnG 102 +GvL+ +++++ ++++ +++ + Y+ +lg++t+++ + +++ +++ ++++ +++++++++++ + ++ + v+ + tr|Q29QN4|Q29QN4_DROME 109 HTSGVLLFGINRGAGQSSSIRKNRpvRVYHLIGHLGIATENNLPDSRVIMRSNHRHVSADRISSLAASMQASHQRKMFELCGVDLQT 195 679*********99999999877668***********************9988********************************** PP -------------------EEEEEEEEEEEEEEEEEEEEEESTTT CS TruB_N 103 krlyelaregeeverkkrkvtiyelelleveepevelevecskGt 147 +++yela +g + +y ++l+++++p+++le++ + t tr|Q29QN4|Q29QN4_DROME 196 QEAYELACRGLLRPADDSQPVVYGIKLIHFDRPHFTLELHAINET 240 *************************************99876655 PP >> tr|Q9V3A8|Q9V3A8_DROME Ergic53, isoform A OS=Drosophila melanogaster OX=7227 GN=ergic53 PE=1 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 9.4 0.0 0.00029 1 77 114 .. 137 174 .. 112 178 .. 0.78 2 ? 0.3 2.8 0.18 6.1e+02 35 84 .. 273 322 .. 263 332 .. 0.82 Alignments for each domain: == domain 1 score: 9.4 bits; conditional E-value: 0.00029 HHHHHHHC-EEE-------------------------- CS TruB_N 77 leevlkkftGeieqvppmySAvkvnGkrlyelaregee 114 l ++++f + + p +SAv +G++ly+ a g++ tr|Q9V3A8|Q9V3A8_DROME 137 LAIMFDSFDNDNKHNNPYISAVLNDGTKLYDHAEDGTT 174 5667899**************************99975 PP == domain 2 score: 0.3 bits; conditional E-value: 0.18 HHCTTS-EEEEEEEEETEEETTSBTTS-EEE-.--TT--HHHHHHHHHHC CS TruB_N 35 eylleadKeYeatirlgaetdteDaegeileeekveklteekleevlkkf 84 e l+++ KeY+ +++ +++ ++D+ e +ee++e++ e++ ++ l+++ tr|Q9V3A8|Q9V3A8_DROME 273 EKLTQEYKEYQDKLEKQKQEYKKDHPDEHKDEEDWEEFYESENQRELRQI 322 5688999*******************999999889999888877777665 PP >> tr|A1ZAJ6|A1ZAJ6_DROME GEO02592p1 OS=Drosophila melanogaster OX=7227 GN=Dmel\CG15919 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.9 0.0 4.7e-05 0.16 9 45 .. 43 78 .. 37 88 .. 0.86 Alignments for each domain: == domain 1 score: 11.9 bits; conditional E-value: 4.7e-05 EEESS---TT-EEEEEEEEGGGGGGHHHCTTS-EEEE CS TruB_N 9 vGhtGtLDplatGvLvvavgeaTklleylleadKeYe 45 v ht tLDp+a G v + + ++y + K Ye tr|A1ZAJ6|A1ZAJ6_DROME 43 VSHTDTLDPQADGFWVNKTTWNARWVKYWRAK-KIYE 78 88************999999999999998755.7776 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (149 nodes) Target sequences: 13786 (7394964 residues searched) Passed MSV filter: 478 (0.0346729); expected 275.7 (0.02) Passed bias filter: 280 (0.0203105); expected 275.7 (0.02) Passed Vit filter: 24 (0.0017409); expected 13.8 (0.001) Passed Fwd filter: 5 (0.000362687); expected 0.1 (1e-05) Initial search space (Z): 13786 [actual number of targets] Domain search space (domZ): 4 [number of targets reported over threshold] # CPU time: 0.11u 0.01s 00:00:00.12 Elapsed: 00:00:00.07 # Mc/sec: 15740.71 // [ok]