# 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/Natronomonas_moolapensis.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1e-28 97.6 0.0 1.4e-28 97.2 0.0 1.2 1 tr|M1Y1N6|M1Y1N6_NATM8 tRNA (Cytidine/uridine-2'-O-)-methylt Domain annotation for each sequence (and alignments): >> tr|M1Y1N6|M1Y1N6_NATM8 tRNA (Cytidine/uridine-2'-O-)-methyltransferase OS=Natronomonas moolapensis (strain DSM 18674 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 97.2 0.0 5.1e-32 1.4e-28 2 142 .] 4 156 .. 3 156 .. 0.79 Alignments for each domain: == domain 1 score: 97.2 bits; conditional E-value: 5.1e-32 EEEEEES-SSHHHHHHHHHHHHHTT-SEEEEESGSSCTCGSHHHHHHCTTHHHHHHEEEESSHHHHHHHHCCCTEEEEEEST..T-S CS SpoU_methylase 2 livvlveiedpgNlGaiaRtakafgaselvlvkprsadllnpeakraaagalevlelkvvdsleealkelkeegfalvatsl..kgk 86 ++v++v+ e+pgN+G+iaR++k+fg+sel+lv+p + d++ + a++ ++vl+ ++ +e + +++ + a++at++ k tr|M1Y1N6|M1Y1N6_NATM8 4 IVVAVVGAETPGNVGSIARSMKNFGLSELYLVDPPELDPDGEAYGFAGHAREDVLPNATEVGFEYLTERFYTV--ACTATTNqdARK 88 7899*********************************999988777888889999944444555555555554..445555443333 PP EGC.....CSHHH.......TTEEEEEEBTTTBS-HHHHCTSSEEEE----TTSSS--HHHHHHHHHH CS SpoU_methylase 87 sls.....eidlk.......ekvalvfGnEdkGlsnelleladrlvkIpmsgeveSLNvavAvailly 142 +++ ++dl+ v++vfG+E+ Gl+n++l+ +dr+++Ip sge++ LN+a+A+ ++ly tr|M1Y1N6|M1Y1N6_NATM8 89 HVRypfvePADLPrelegvdADVCVVFGRERVGLTNDELAVLDRICSIPASGEYPVLNLAQAATVVLY 156 34343333223337777888**********************************************98 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (142 nodes) Target sequences: 2723 (809592 residues searched) Passed MSV filter: 67 (0.0246052); expected 54.5 (0.02) Passed bias filter: 62 (0.022769); expected 54.5 (0.02) Passed Vit filter: 10 (0.00367242); expected 2.7 (0.001) Passed Fwd filter: 1 (0.000367242); expected 0.0 (1e-05) Initial search space (Z): 2723 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.02u 0.00s 00:00:00.02 Elapsed: 00:00:00.01 # Mc/sec: 11496.21 // [ok]