# 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/Drosophila_melanogaster.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2.4e-102 341.5 0.0 3.3e-102 341.1 0.0 1.2 1 tr|Q9VL60|Q9VL60_DROME RE01590p OS=Drosophila melanogaster O ------ inclusion threshold ------ 0.04 12.8 0.0 0.05 12.5 0.0 1.2 1 sp|O97067|PTH2_DROME Probable peptidyl-tRNA hydrolase 2 OS Domain annotation for each sequence (and alignments): >> tr|Q9VL60|Q9VL60_DROME RE01590p OS=Drosophila melanogaster OX=7227 GN=anon-WO0118547.208 PE=1 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 341.1 0.0 4.8e-106 3.3e-102 2 275 .] 157 441 .. 156 441 .. 0.96 Alignments for each domain: == domain 1 score: 341.1 bits; conditional E-value: 4.8e-106 Rsm22 2 eesiaYalarlpatyavlrrvLdelaervpdflPkslldvgegpgtalwaaselwe.eleevlvidaseelleiskelaadapelke 87 +es++Ya+ar +++yavl+rvL+ela+r+++f+P+s++d+g+g+gt++w+aselw+ +++e++++d+s+e++eis+ +++d +e+k+ tr|Q9VL60|Q9VL60_DROME 157 YESLVYAIARGAQEYAVLKRVLTELAQRDEQFRPRSFFDFGSGIGTGMWVASELWReHIFEYYNVDRSREMNEISELILRDGHENKQ 243 789*****************************************************9****************************** PP Rsm22 88 aalrrsvirealeleea..dLviisyvLleled.esreklvknlWakaskilvivEeGtpaGfrrvleaReal..kaagfhivAPCP 169 +alr++++r+ l e+ dLviis++L+ele+ e+re++++nlW+k++ +lvivEeGt++G+++v+eaR++l +++++h+vAPCP tr|Q9VL60|Q9VL60_DROME 244 IALRNVFYRQFLPAIETnyDLVIISHTLFELENkEQREDVLRNLWRKCDGYLVIVEEGTRRGSELVNEARQFLleQEQEGHTVAPCP 330 ********76665555589**************9***************************************999*********** PP Rsm22 170 HeaaCPlvate...dwChFsqrva..Rlslhrlaksasvpvedekfsy..vakerqprasa.aaRvvrptkvrskkvlidlCsedet 248 H+++CP+++++ ++C+F++++a Rl++++la+ +++sy ++k++q+++s+ ++R+vrpt vrsk+++++lC+e++t tr|Q9VL60|Q9VL60_DROME 331 HDLVCPRLRDLadrTPCNFPISYAplRLGSDFLADR-----HTSLYSYaiLKKGSQSDSSSsWPRIVRPTLVRSKHAICRLCTEKGT 412 ********9988889********************7.....5799***99999999887656************************* PP Rsm22 249 lqrlvvtK.rkG.eaykaArkaeWGDrlp 275 lq++++tK ++G +ay++A++++WGDrlp tr|Q9VL60|Q9VL60_DROME 413 LQEVIFTKsKHGkSAYSCAKVSRWGDRLP 441 ********99**99************986 PP >> sp|O97067|PTH2_DROME Probable peptidyl-tRNA hydrolase 2 OS=Drosophila melanogaster OX=7227 GN=CG1307 PE=2 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.5 0.0 7.2e-06 0.05 144 175 .. 69 100 .. 66 167 .. 0.81 Alignments for each domain: == domain 1 score: 12.5 bits; conditional E-value: 7.2e-06 Rsm22 144 paGfrrvleaRealkaagfhivAPCPHeaaCP 175 + f++vl++R+ lk +++i+A C H a+ sp|O97067|PTH2_DROME 69 NDNFKMVLVVRNDLKMGKGKIAAQCGHGAVGA 100 6779***********************99865 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (275 nodes) Target sequences: 13786 (7394964 residues searched) Passed MSV filter: 361 (0.026186); expected 275.7 (0.02) Passed bias filter: 244 (0.0176991); expected 275.7 (0.02) Passed Vit filter: 18 (0.00130567); expected 13.8 (0.001) Passed Fwd filter: 2 (0.000145075); expected 0.1 (1e-05) Initial search space (Z): 13786 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.15u 0.01s 00:00:00.16 Elapsed: 00:00:00.06 # Mc/sec: 33893.58 // [ok]