# 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/Ciona_savignyi.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.8e-53 180.3 0.0 4.6e-53 179.7 0.0 1.3 1 tr|H2Z5Q2|H2Z5Q2_CIOSA Uncharacterized protein OS=Ciona savi ------ inclusion threshold ------ 0.053 12.2 0.0 0.081 11.6 0.0 1.2 1 tr|H2YTR8|H2YTR8_CIOSA Methyltransferase-like protein OS=Cio 0.13 10.9 4.1 0.11 11.2 3.0 1.3 1 tr|H2YUM2|H2YUM2_CIOSA Uncharacterized protein OS=Ciona savi Domain annotation for each sequence (and alignments): >> tr|H2Z5Q2|H2Z5Q2_CIOSA Uncharacterized protein OS=Ciona savignyi OX=51511 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 179.7 0.0 1.2e-56 4.6e-53 3 273 .. 111 401 .. 109 403 .. 0.87 Alignments for each domain: == domain 1 score: 179.7 bits; conditional E-value: 1.2e-56 Rsm22 3 esiaYalarlpatyavlrrvLdelaervpdflPkslldvgegpgtalwaaselwe.eleevlvidaseelleiskelaadapelkea 88 +s a++lar p+ ya rvL +l++ + f+Pks+ld+g+g+gt++waa++lw+ +l++++ +d se+++++s++l a + tr|H2Z5Q2|H2Z5Q2_CIOSA 111 NSHAFLLARSPGCYASSFRVLSNLKHEDFSFTPKSMLDFGSGTGTSMWAANKLWKeSLKQYVCVDSSEHMNRLSHTLHTTADTESDL 197 789****************************************************99*******************98877665554 PP Rsm22 89 alrrsvireal.eleea.dLviisyvLleled.esreklvknlWakaskilvivEeGtpaGfrrvleaReal............... 157 +ir+ l ++ + dLvi sy+L e ++r +l+k lW+k+s++lvi+E+G +G++ ++eaR+ l tr|H2Z5Q2|H2Z5Q2_CIOSA 198 -PPGIFIRKFLpSFKMTyDLVISSYTLSEQPSlKERMILLKLLWEKTSDYLVIIEHGNFNGYQIIMEARNLLvsgiegikgqktwkn 283 .67777755541455449************999**********************************************99988866 PP Rsm22 158 kaagfhivAPCPHeaaCPlvatedwChFsqrvaRlslhrlaksasvpvedekfsy..vakerqprasa.....aaRvvrptkvrskk 237 +++ +vAPCPH+ CP+ +++ C Fsq ++ s + e fsy ++k++ ++++ ++Rvv ++ ++ tr|H2Z5Q2|H2Z5Q2_CIOSA 284 WPDNATVVAPCPHDGPCPMWGSKTSCIFSQTYNVPSFAKMGTA------SESFSYvvLKKNSS-VENTqrthsWPRVVGGPNEQKSC 363 6789****************************98887766654......58899974455443.333246666************** PP Rsm22 238 vlidlCsedetlqrlvvtK.rkGeayk.aArkaeWGDr 273 v + +C++d+ l++l +tK +G++ k + +k + GD tr|H2Z5Q2|H2Z5Q2_CIOSA 364 VDCSVCTTDGVLENLHITKkIHGTNLKsMVKKCHRGDY 401 *******************66**666526788888886 PP >> tr|H2YTR8|H2YTR8_CIOSA Methyltransferase-like protein OS=Ciona savignyi OX=51511 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.6 0.0 2.1e-05 0.081 22 130 .. 66 176 .. 61 190 .. 0.80 Alignments for each domain: == domain 1 score: 11.6 bits; conditional E-value: 2.1e-05 Rsm22 22 vLdelaervpdflPkslldvgegpgtalwaaselweeleevlvidas.eelleiske..laadapelkeaalrrsvirealeleead 105 +e++er++ + ll+ g g+g a++ + +++l+ + d s ++++++ ++ +++ ++ + ++ + + ++ +++d tr|H2YTR8|H2YTR8_CIOSA 66 EFNEIEERIKSNEHTVLLEAGCGVGNAIFPLMDKYSNLFV-YACDFSpRAVQHVKDNpgFSKKNCNAFQCDISTQSLCDNIQENSVD 151 5678899999999999999999999999999999999985.556666255667766533566666666666777778899999999* PP Rsm22 106 LviisyvLleledesreklvknlWa 130 Lv++ +vL +++++ e + +n+ + tr|H2YTR8|H2YTR8_CIOSA 152 LVTLIFVLSAIHPNKFEFVFENIHK 176 *****************99999855 PP >> tr|H2YUM2|H2YUM2_CIOSA Uncharacterized protein OS=Ciona savignyi OX=51511 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.2 3.0 2.8e-05 0.11 166 207 .. 26 67 .. 15 84 .. 0.86 Alignments for each domain: == domain 1 score: 11.2 bits; conditional E-value: 2.8e-05 Rsm22 166 APCPHeaaCPlvatedwChFsqrvaRlslhrlaksasvpved 207 APC e++CP + ++C F + + l+ + ++ s + e tr|H2YUM2|H2YUM2_CIOSA 26 APCRDESICPQLIASTPCTFLKLSEQLRCPVSCNRCSSASES 67 9*******************9999999998888876766665 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (275 nodes) Target sequences: 11592 (4231397 residues searched) Passed MSV filter: 305 (0.0263112); expected 231.8 (0.02) Passed bias filter: 214 (0.018461); expected 231.8 (0.02) Passed Vit filter: 13 (0.00112146); expected 11.6 (0.001) Passed Fwd filter: 3 (0.000258799); expected 0.1 (1e-05) Initial search space (Z): 11592 [actual number of targets] Domain search space (domZ): 3 [number of targets reported over threshold] # CPU time: 0.09u 0.00s 00:00:00.09 Elapsed: 00:00:00.04 # Mc/sec: 29090.85 // [ok]