# 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/SpoU_methylase.hmm.txt # target sequence database: proteomes/Anolis_carolinensis.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: SpoU_methylase [M=142] Accession: PF00588.18 Description: SpoU rRNA Methylase family 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2e-30 105.8 0.0 3.1e-30 105.2 0.0 1.2 1 tr|H9G6C9|H9G6C9_ANOCA Mitochondrial rRNA methyltransferase 3.8e-23 82.3 0.0 1.8e-13 50.9 0.0 2.2 2 tr|H9GMT3|H9GMT3_ANOCA Uncharacterized protein OS=Anolis car 5e-23 81.9 0.0 1.2e-22 80.7 0.0 1.6 1 tr|G1K946|G1K946_ANOCA TAR (HIV-1) RNA binding protein 1 OS= ------ inclusion threshold ------ 0.011 16.0 0.0 0.25 11.5 0.0 2.4 2 tr|H9G369|H9G369_ANOCA Uncharacterized protein OS=Anolis car 0.051 13.8 0.0 0.088 13.0 0.0 1.3 1 tr|H9G6C2|H9G6C2_ANOCA Short chain dehydrogenase/reductase f Domain annotation for each sequence (and alignments): >> tr|H9G6C9|H9G6C9_ANOCA Mitochondrial rRNA methyltransferase 1 OS=Anolis carolinensis OX=28377 GN=MRM1 PE=4 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 105.2 0.0 8.5e-34 3.1e-30 4 142 .] 159 308 .. 157 308 .. 0.90 Alignments for each domain: == domain 1 score: 105.2 bits; conditional E-value: 8.5e-34 EEEES-SSHHHHHHHHHHHHHTT-SEEEEESGSSCTCGSHHHHHHCTTHHHHHHEEEESSHHHHHHHHCCCTEEEEEESTT-SEGC. CS SpoU_methylase 4 vvlveiedpgNlGaiaRtakafgaselvlvkprsadllnpeakraaagalevlelkvvdsleealkelkeegfalvatslkgksls. 89 ++l++i+dp+NlGa++R+a +g++++v+++++s +l+p + +a++ga evl++ +d+l +lk +eg+++++t +++k + tr|H9G6C9|H9G6C9_ANOCA 159 LALEGIQDPMNLGAVLRSAHFLGVDGIVMSRKNSC-PLTPVVSKASSGAMEVLDVFSADDLPGFLKVKSKEGWEVLGTVAHAKAQDg 244 789**************************998888.******************************************997775443 PP .......CSHHHTTEEEEEEBTTTBS-HHHHCTSSEEEE-....---TTSSS--HHHHHHHHHH CS SpoU_methylase 90 .......eidlkekvalvfGnEdkGlsnelleladrlvkI....pmsgeveSLNvavAvailly 142 +++++++++lv+GnE++Gls+e+ +l++r+++I ++ +eSLNv+vA++ill+ tr|H9G6C9|H9G6C9_ANOCA 245 ipvvsclDFRWMKPTILVLGNEGHGLSPETKSLCRRMLTIapgrELEPGIESLNVSVAAGILLH 308 344456677999****************************4333368999***********985 PP >> tr|H9GMT3|H9GMT3_ANOCA Uncharacterized protein OS=Anolis carolinensis OX=28377 GN=mrm3 PE=4 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 50.9 0.0 4.8e-17 1.8e-13 2 63 .. 191 249 .. 190 268 .. 0.91 2 ! 29.2 0.0 2.4e-10 8.8e-07 95 142 .] 336 387 .. 313 387 .. 0.85 Alignments for each domain: == domain 1 score: 50.9 bits; conditional E-value: 4.8e-17 EEEEEES-SSHHHHHHHHHHHHHTT-SEEEEESGSSCTCGSHHHHHHCTTHHHHHHEEEESS CS SpoU_methylase 2 livvlveiedpgNlGaiaRtakafgaselvlvkprsadllnpeakraaagalevlelkvvds 63 l +++++i+dpgNlG+i+R+a+ +g+s+++lvk ++d ++p+++ra++ga++ l++ v++ tr|H9GMT3|H9GMT3_ANOCA 191 LSLICDNIRDPGNLGTILRSAAGAGCSRVLLVK-GCVDAWEPKVLRAGMGAHFRLPI--VSN 249 6799****************************9.**********************5..444 PP == domain 2 score: 29.2 bits; conditional E-value: 2.4e-10 TTEEEEEEBTTTBS-HHHHCTSS....EEEE----TTSSS--HHHHHHHHHH CS SpoU_methylase 95 ekvalvfGnEdkGlsnellelad....rlvkIpmsgeveSLNvavAvailly 142 +va+v+G E++G+s+e+++la+ r + Ip+++ v+SLN a+A++ill+ tr|H9GMT3|H9GMT3_ANOCA 336 LSVAVVIGGETHGVSEEAMRLAQdtggRRLLIPIAPGVDSLNAAMAASILLF 387 689*******************955556899******************985 PP >> tr|G1K946|G1K946_ANOCA TAR (HIV-1) RNA binding protein 1 OS=Anolis carolinensis OX=28377 GN=TARBP1 PE=4 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 80.7 0.0 3.2e-26 1.2e-22 2 127 .. 1426 1553 .. 1425 1556 .. 0.90 Alignments for each domain: == domain 1 score: 80.7 bits; conditional E-value: 3.2e-26 EEEEEES-SSHHHHHHHHHHHHHTT-SEEEEESGSSCTCGSHHHHHHCTTHHHHHHEEEES..SHHHHHHHHCCCTEEEEEESTT CS SpoU_methylase 2 livvlveiedpgNlGaiaRtakafgaselvlvkprsadllnpeakraaagalevlelkvvd..sleealkelkeegfalvatslk 84 livv +++p+NlG ++Rt+++fg+s+l++ + + + +++++ ++++a+++l l +v+ +l+++l++ k eg++++++++ tr|G1K946|G1K946_ANOCA 1426 LIVVASLLDKPTNLGGLCRTCEIFGVSALIVGS-IHY-IHDKQFQYLSVSAEQWLSLLEVKpfQLADYLEQKKIEGYTIIGVEQT 1508 7999999*******************9999966.554.899***************98886347899*****************6 PP -SEGC..CSHHHTTEEEEEEBTTTBS-HHHHCTSSEEEE----TT CS SpoU_methylase 85 gksls..eidlkekvalvfGnEdkGlsnelleladrlvkIpmsge 127 +ks++ e++++ek +l++GnE +G++ +l++++d +v+Ip++g tr|G1K946|G1K946_ANOCA 1509 AKSCDltEYSFPEKSLLLLGNEHEGIPANLIQQLDTCVEIPQQGV 1553 6655554***********************************985 PP >> tr|H9G369|H9G369_ANOCA Uncharacterized protein OS=Anolis carolinensis OX=28377 GN=spout1 PE=4 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.5 0.0 6.9e-05 0.25 16 41 .. 43 68 .. 42 102 .. 0.87 2 ? 1.4 0.0 0.089 3.3e+02 75 102 .. 227 257 .. 208 269 .. 0.83 Alignments for each domain: == domain 1 score: 11.5 bits; conditional E-value: 6.9e-05 HHHHHHHHHTT-SEEEEESGSSCTCG CS SpoU_methylase 16 GaiaRtakafgaselvlvkprsadll 41 G iaR++++f ++e+v+ ++r++d tr|H9G369|H9G369_ANOCA 43 GQIARSCAVFSVDEVVVFDERGEDAK 68 99*****************9998654 PP == domain 2 score: 1.4 bits; conditional E-value: 0.089 TEEE.EEESTT-SEGCCSHHH..TTEEEEEE CS SpoU_methylase 75 gfal.vatslkgkslseidlk..ekvalvfG 102 g++l v+ts++g+++++++l+ +++++vfG tr|H9G369|H9G369_ANOCA 227 GYDLsVGTSERGTPVDQATLPafRHALIVFG 257 5666589999******999999999*****9 PP >> tr|H9G6C2|H9G6C2_ANOCA Short chain dehydrogenase/reductase family 42E, member 1 OS=Anolis carolinensis OX=28377 GN=S # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.0 0.0 2.4e-05 0.088 17 87 .. 81 148 .. 71 162 .. 0.83 Alignments for each domain: == domain 1 score: 13.0 bits; conditional E-value: 2.4e-05 HHHHHHHHTT-SEEEEESGSSCTCGSHHHHHHCTTHHHHHHEEEESSHHHHHHHHCCCTEEEEEESTT-SE CS SpoU_methylase 17 aiaRtakafgaselvlvkprsadllnpeakraaagalevlelkvvdsleealkelkeegfalvatslkgks 87 +++R+a+a+g+++lv++++ ++ ++ +++r + +++ +l+l+ +++ + k+l e a++++++ g++ tr|H9G6C2|H9G6C2_ANOCA 81 NVLRAARAAGVRRLVYTSTYNV-VFGGQVVRGGDESWPYLPLRLHPDHYSRTKALAEM--AVLGADGAGRE 148 79****************9998.**********************9999888887765..99999884443 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (142 nodes) Target sequences: 18528 (9796552 residues searched) Passed MSV filter: 391 (0.0211032); expected 370.6 (0.02) Passed bias filter: 359 (0.0193761); expected 370.6 (0.02) Passed Vit filter: 30 (0.00161917); expected 18.5 (0.001) Passed Fwd filter: 5 (0.000269862); expected 0.2 (1e-05) Initial search space (Z): 18528 [actual number of targets] Domain search space (domZ): 5 [number of targets reported over threshold] # CPU time: 0.11u 0.01s 00:00:00.12 Elapsed: 00:00:00.05 # Mc/sec: 27822.21 // [ok]