# 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/Rsm22.hmm.txt # target sequence database: proteomes/Strongylocentrotus_purpuratus.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: Rsm22 [M=275] Accession: PF09243.9 Description: Mitochondrial small ribosomal subunit Rsm22 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 3e-63 214.3 1.2 1.6e-30 107.0 0.0 2.1 2 tr|W4Z116|W4Z116_STRPU Uncharacterized protein OS=Strongyloc Domain annotation for each sequence (and alignments): >> tr|W4Z116|W4Z116_STRPU Uncharacterized protein OS=Strongylocentrotus purpuratus OX=7668 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 107.0 0.0 5.5e-35 1.6e-30 3 158 .. 176 336 .. 174 346 .. 0.92 2 ! 105.8 0.8 1.3e-34 3.8e-30 156 275 .] 363 486 .. 353 486 .. 0.89 Alignments for each domain: == domain 1 score: 107.0 bits; conditional E-value: 5.5e-35 Rsm22 3 esiaYalarlpatyavlrrvLdelaervpdflPkslldvgegpgtalwaaselwe.eleevlvidaseelleiskelaadapelkea 88 +s+aY+lar+p +yavl+rvL+e+++ +p+++P ++ld+g+g gt++waa+ +w ++e+++++ +++++++s+ l++ p +e tr|W4Z116|W4Z116_STRPU 176 KSLAYMLARMPPNYAVLYRVLNEIKQGDPEYRPSTILDFGSGLGTTTWAAHSVWGdAVKEYYMVEVAKAMHQLSERLLRVEPGKEEL 262 799****************************************************999*************************9999 PP Rsm22 89 alrrsvirealeleea...dLviisyvLleled.esreklvknlWakaskilvivEeGtpaGfrrvleaRealk 158 ++ r + + + + +y Llel ++r k v+ lW+ ++++lv++E+G + + ++eaR+a+ tr|W4Z116|W4Z116_STRPU 263 LIPNVYFRFFMPVSLRarySMTVSAYSLLELPSfQDRVKSVRSLWRLTDDYLVLIENGSYESYLALMEARDAIL 336 8888888665555443567888889*********************************************9873 PP == domain 2 score: 105.8 bits; conditional E-value: 1.3e-34 Rsm22 156 alkaagfhivAPCPHeaaCPlvate..dwChFsqrvaRlslhrlaksasvpvedekfsy..vakerqprasaaaRvvrptkvrskkv 238 ++ +g+h+ APCPH+ CP+++++ +C+F q+++ l++ + +s v e+f y ++k++ ++++++Rvv+p ++s+++ tr|W4Z116|W4Z116_STRPU 363 QIDRQGGHVFAPCPHDQPCPRLQDGsgIPCNFEQSYHPLTQFIGK--RSSMVSTERFTYmvLKKGSRTNKHSWPRVVKPVLKKSRHI 447 56789******************99889*********99765544..4566777788887799999999999*************** PP Rsm22 239 lidlCsedetlqrlvvtK.rkG.eaykaArkaeWGDrlp 275 +++lC +d++l++ v+tK ++ ++y++A+ ++WGD lp tr|W4Z116|W4Z116_STRPU 448 ICRLCCDDGSLREAVFTKaKHSkDMYRCAKYTNWGDMLP 486 ******************9998899***********985 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (275 nodes) Target sequences: 28594 (14266044 residues searched) Passed MSV filter: 781 (0.0273134); expected 571.9 (0.02) Passed bias filter: 570 (0.0199343); expected 571.9 (0.02) Passed Vit filter: 47 (0.0016437); expected 28.6 (0.001) Passed Fwd filter: 1 (3.49724e-05); expected 0.3 (1e-05) Initial search space (Z): 28594 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.28u 0.01s 00:00:00.29 Elapsed: 00:00:00.10 # Mc/sec: 39231.62 // [ok]