# 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/Monodelphis_domestica.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1.1e-61 208.8 0.0 1.5e-61 208.4 0.0 1.2 1 tr|F6PFR6|F6PFR6_MONDO Methyltransferase like 17 OS=Monodelp Domain annotation for each sequence (and alignments): >> tr|F6PFR6|F6PFR6_MONDO Methyltransferase like 17 OS=Monodelphis domestica OX=13616 GN=METTL17 PE=4 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 208.4 0.0 7e-66 1.5e-61 3 275 .] 175 462 .. 173 462 .. 0.92 Alignments for each domain: == domain 1 score: 208.4 bits; conditional E-value: 7e-66 Rsm22 3 esiaYalarlpatyavlrrvLdelaervpdflPkslldvgegpgtalwaaselwe.eleevlvidaseelleiskelaadapelkea 88 s++Y++arl + ya++ r+ +e+++rvpdf+P++l+d+g+g+g ++waa w +l+e+l +d s+++l++++ l++ e + tr|F6PFR6|F6PFR6_MONDO 175 LSLVYMAARLDGGYAAVSRAFHEIQARVPDFKPQTLMDFGSGTGSVTWAAYSSWGqSLREYLCVDSSAAMLALAERLLKGGLESGKP 261 689***************************************************999****************************** PP Rsm22 89 alrrsvirealeleea...dLviisyvLleled.esreklvknlWakaskilvivEeGtpaGfrrvleaReal.........kaagf 162 + ++ r+ l + ++ d+v+ ++ L el ++r+++++ lW+k+ + lv+vE+Gt++G++++++aR+ + + + tr|F6PFR6|F6PFR6_MONDO 262 HIPGVFFRQFLPVSPKvqfDVVVSAFSLSELPSmADRTEVIQTLWRKTGHFLVLVENGTKSGHHLLMDARDLVlrgqekfplDPRPG 348 *******99999999899***************9*************************************9889877777667889 PP Rsm22 163 hivAPCPHeaaCPlvate..dwChFsqrvaRlslhrlaksasvpvedekfsyvakerqprasa..aaRvvrptkvrskkvlidlCse 245 h+ APCPHe+ CP + ++C Fsq ++ ++ +ekfs v +r+ ++ +aR+ +p +r ++v ++lC tr|F6PFR6|F6PFR6_MONDO 349 HVFAPCPHELPCPKLSASrpQPCSFSQAYSTIPFSWN-----HSHKEEKFSLVILSRGRPGEYhqWARLTQPVLKRPRHVHCHLCFP 430 *************9988889*******9986655443.....35689*****5555555444444********************** PP Rsm22 246 detlqrlvvtK.rkG.eaykaArkaeWGDrlp 275 d++lq++vvt r+G ++y++Ar ++WGDrlp tr|F6PFR6|F6PFR6_MONDO 431 DGHLQHVVVTPrRHGrDLYRCARASSWGDRLP 462 **********989**999***********986 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (275 nodes) Target sequences: 21272 (11266248 residues searched) Passed MSV filter: 511 (0.0240222); expected 425.4 (0.02) Passed bias filter: 383 (0.0180049); expected 425.4 (0.02) Passed Vit filter: 37 (0.00173938); expected 21.3 (0.001) Passed Fwd filter: 1 (4.70102e-05); expected 0.2 (1e-05) Initial search space (Z): 21272 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.22u 0.01s 00:00:00.23 Elapsed: 00:00:00.09 # Mc/sec: 34424.65 // [ok]