# 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/MT-A70.hmm.txt # target sequence database: proteomes/Trichoplax_adhaerens.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: MT-A70 [M=172] Accession: PF05063.13 Description: MT-A70 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 3e-16 59.2 0.0 1.2e-13 50.7 0.0 2.4 2 tr|B3S3T4|B3S3T4_TRIAD Uncharacterized protein OS=Trichoplax Domain annotation for each sequence (and alignments): >> tr|B3S3T4|B3S3T4_TRIAD Uncharacterized protein OS=Trichoplax adhaerens OX=10228 GN=TRIADDRAFT_58838 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 7.6 0.0 0.00019 2.2 22 53 .. 173 204 .. 163 209 .. 0.85 2 ! 50.7 0.0 1.1e-17 1.2e-13 80 172 .] 200 299 .. 198 299 .. 0.81 Alignments for each domain: == domain 1 score: 7.6 bits; conditional E-value: 0.00019 SS--B-HHHHHTS-GGGCECCEEEEEEEEETT CS MT-A70 22 pyktlneeelkaLpipelldkrslvflWvtnr 53 +y+ l++++l++Lpi ell+ ++lv+l +++ tr|B3S3T4|B3S3T4_TRIAD 173 KYHWLTTDQLLQLPITELLARHGLVVLTTGGE 204 578899*******************9976665 PP == domain 2 score: 50.7 bits; conditional E-value: 1.1e-17 -TTSCBS-BS-TTSSSBB-EEEEEEEEES---TTT.TTB-TTS........SBSEEEE.E-CTTTSS-THHHHHHHHHC.TCS-EEE CS MT-A70 80 nklgepiatlesgavlkkskEhllvgvkgnvkrstdkiinknv........dtdliis.perehsrKPpelfellehll.esarkLe 156 ++ ge +++l+s +kk++E ll+gvk+ ++ + ++ + ++i s p+ +hs+KPp l + ++ e+ ++Le tr|B3S3T4|B3S3T4_TRIAD 200 TTGGEFVYDLDSS--HKKPYEILLIGVKSGLSEDLPSACSSELpesiaafpKHNVICSiPSLVHSHKPP-LDGIFARFIrEQPNCLE 283 5567788887776..9*************998877654444434455778778888889**********.89999999956666*** PP ET--G--TTEEEEETT CS MT-A70 157 lFaRslrpgwltvGne 172 +FaRsl p+w+++Gne tr|B3S3T4|B3S3T4_TRIAD 284 MFARSLIPHWTSWGNE 299 **************97 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (172 nodes) Target sequences: 11518 (5222075 residues searched) Passed MSV filter: 376 (0.0326446); expected 230.4 (0.02) Passed bias filter: 252 (0.0218788); expected 230.4 (0.02) Passed Vit filter: 17 (0.00147595); expected 11.5 (0.001) Passed Fwd filter: 2 (0.000173641); expected 0.1 (1e-05) Initial search space (Z): 11518 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.09u 0.00s 00:00:00.09 Elapsed: 00:00:00.04 # Mc/sec: 22454.92 // [ok]