# 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/Mus_musculus.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 3.3e-61 207.3 0.0 4.2e-61 206.9 0.0 1.1 1 sp|Q3U2U7|MET17_MOUSE Methyltransferase-like protein 17, mit ------ inclusion threshold ------ 9.6 5.7 6.3 4.5 6.8 1.1 2.2 2 sp|Q9WUS4|CXA10_MOUSE Gap junction alpha-10 protein OS=Mus m Domain annotation for each sequence (and alignments): >> sp|Q3U2U7|MET17_MOUSE Methyltransferase-like protein 17, mitochondrial OS=Mus musculus OX=10090 GN=Mettl17 PE=2 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 206.9 0.0 3.8e-65 4.2e-61 4 275 .] 156 442 .. 153 442 .. 0.93 Alignments for each domain: == domain 1 score: 206.9 bits; conditional E-value: 3.8e-65 Rsm22 4 siaYalarlpatyavlrrvLdelaervpdflPkslldvgegpgtalwaaselwe.eleevlvidaseelleiskelaadapelkeaal 90 s++Y++arl + +a++ r+ +e+++r+p+f+P++l+d+g+g+g + waa+ +w +l+e++ +d s+++l ++++l++ e + sp|Q3U2U7|MET17_MOUSE 156 SLIYMAARLDGGFAAVFRAFHEIQARIPEFRPQTLMDFGSGTGSVAWAAHRTWGqSLREYVCVDSSAAMLGLAEKLLKGGSESGKPCI 243 899***************************************************99******************************** PP Rsm22 91 rrsvirealeleea...dLviisyvLleled.esreklvknlWakaskilvivEeGtpaGfrrvleaReal.........kaagfhiv 165 ++ r+ l + ++ d+v+ +y L el ++r ++++nlW+k+s+ lv+vE+Gt+aG+r++++aR+ + + + sp|Q3U2U7|MET17_MOUSE 244 PGVFFRQFLPVSPKvqfDVVVSAYALSELPSrADRIEVIQNLWRKTSHFLVLVENGTKAGHRLLMDARNLVlgekeksplDLRPSFVF 331 *****99999999899*************999*************************************98898776666667789** PP Rsm22 166 APCPHeaaCPlvate..dwChFsqrvaRlslhrlaksasvpvedekfsy..vakerqprasaaaRvvrptkvrskkvlidlCsedetl 249 APCPHe+ CP ++ C Fsq ++ + k +e fs +a+++ ++a+ ++R+ +p +r ++v ++lC d+++ sp|Q3U2U7|MET17_MOUSE 332 APCPHELPCPQLNASksLACSFSQAYHPIPFNWNKKP-----KEEIFSMviLARGSPKEANRWPRITQPVLKRPRHVHCHLCCPDGHM 414 **********988876789********9998887775.....68999997667777788888************************** PP Rsm22 250 qrlvvtK.rkG.eaykaArkaeWGDrlp 275 q+ vvt r+G ++y++Ar+++WGD+lp sp|Q3U2U7|MET17_MOUSE 415 QHAVVTArRHGrDLYRCARVSSWGDLLP 442 ******989**999***********985 PP >> sp|Q9WUS4|CXA10_MOUSE Gap junction alpha-10 protein OS=Mus musculus OX=10090 GN=Gja10 PE=1 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 1.0 0.3 0.023 2.6e+02 164 181 .. 190 207 .. 182 221 .. 0.83 2 ? 6.8 1.1 0.00041 4.5 160 184 .. 330 355 .. 321 377 .. 0.86 Alignments for each domain: == domain 1 score: 1.0 bits; conditional E-value: 0.023 Rsm22 164 ivAPCPHeaaCPlvated 181 + APCP + +C + ++++ sp|Q9WUS4|CXA10_MOUSE 190 TQAPCPNSVDCFVSRPTE 207 57*********9888765 PP == domain 2 score: 6.8 bits; conditional E-value: 0.00041 Rsm22 160 agfhivAPCPHeaaCPlvate.dwCh 184 a+ h+ PCPH+ +++++ ++Ch sp|Q9WUS4|CXA10_MOUSE 330 AQLHVHSPCPHDRSARIQHPGqQPCH 355 6789************9999889**9 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (275 nodes) Target sequences: 22289 (11757923 residues searched) Passed MSV filter: 580 (0.0260218); expected 445.8 (0.02) Passed bias filter: 432 (0.0193818); expected 445.8 (0.02) Passed Vit filter: 35 (0.00157028); expected 22.3 (0.001) Passed Fwd filter: 3 (0.000134596); expected 0.2 (1e-05) Initial search space (Z): 22289 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.24u 0.01s 00:00:00.25 Elapsed: 00:00:00.09 # Mc/sec: 35926.99 // [ok]