# 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/Trypanosoma_brucei.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.4e-28 98.9 0.9 2.7e-16 58.6 0.2 3.3 3 tr|Q385R2|Q385R2_TRYB2 Uncharacterized protein OS=Trypanosom Domain annotation for each sequence (and alignments): >> tr|Q385R2|Q385R2_TRYB2 Uncharacterized protein OS=Trypanosoma brucei brucei (strain 927/4 GUTat10.1) OX=185431 GN=Tb # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 36.1 0.0 2.3e-13 1.9e-09 3 60 .. 252 309 .. 250 357 .. 0.89 2 ! 58.6 0.2 3.2e-20 2.7e-16 57 178 .. 474 601 .. 444 622 .. 0.82 3 ? -0.1 0.0 0.025 2.1e+02 219 250 .. 814 844 .. 797 867 .. 0.80 Alignments for each domain: == domain 1 score: 36.1 bits; conditional E-value: 2.3e-13 Rsm22 3 esiaYalarlpatyavlrrvLdelaervpdflPkslldvgegpgtalwaaselweele 60 ++i Y l r +a ya ++rv el ++ p+f+Pk++ld+g+g+gta+++a+e+++ tr|Q385R2|Q385R2_TRYB2 252 DAIGYFLFRGAAMYAGMHRVFFELSKQLPHFVPKTMLDFGAGTGTAILVAKEVYDPGS 309 7899**************************************************8544 PP == domain 2 score: 58.6 bits; conditional E-value: 3.2e-20 Rsm22 57 eeleevlvidaseelleiskelaadapelkeaalrrsvirealeleeadLviisyvLleled.esreklvknlWakaskilvivEeG 142 + l+ev+ ++ s + ei+ +++d + ++++r ++ e e + dLv+ +y L e++ e+r +v++lW+ ++ +lv vE tr|Q385R2|Q385R2_TRYB2 474 RPLQEVTAVEPSPGMMEIGTMVLHDDVPN--VTWKRYLLPED-EAIQHDLVVAAYSLSEIATsENRRRIVQQLWKMTKGVLVFVEFA 557 56899**************9988765443..45667777554.55567*************99************************ PP Rsm22 143 tpaGfrrvleaReal........kaagfhivAPCPHeaaCPlva 178 + f+ ++eaR+ + + ivAPCPHe CPl + tr|Q385R2|Q385R2_TRYB2 558 NLNNFNILMEARDWIleekdvglWDWQPTIVAPCPHEHRCPLRH 601 *************9977655554455679************975 PP == domain 3 score: -0.1 bits; conditional E-value: 0.025 Rsm22 219 prasaaaRvvrptkvrskkvlidlCsedetlq 250 +a++ Rvvr ++ +kv d C+ ++ l tr|Q385R2|Q385R2_TRYB2 814 RAAPKYLRVVRDPRC-HGKVQADFCTPEGDLV 844 4556688******97.589*******999875 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (275 nodes) Target sequences: 8561 (4332363 residues searched) Passed MSV filter: 257 (0.0300199); expected 171.2 (0.02) Passed bias filter: 182 (0.0212592); expected 171.2 (0.02) Passed Vit filter: 22 (0.00256979); expected 8.6 (0.001) Passed Fwd filter: 1 (0.000116809); expected 0.1 (1e-05) Initial search space (Z): 8561 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.10u 0.00s 00:00:00.10 Elapsed: 00:00:00.06 # Mc/sec: 19856.66 // [ok]