# 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/Gcd10p.hmm.txt # target sequence database: proteomes/Monosiga_brevicollis.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: Gcd10p [M=300] Accession: PF04189.12 Description: Gcd10p family 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-41 142.0 0.0 1.5e-41 141.6 0.0 1.3 1 tr|A9UTS9|A9UTS9_MONBE Predicted protein OS=Monosiga brevico Domain annotation for each sequence (and alignments): >> tr|A9UTS9|A9UTS9_MONBE Predicted protein OS=Monosiga brevicollis OX=81824 GN=23559 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 141.6 0.0 1.6e-45 1.5e-41 26 236 .. 24 210 .. 18 331 .. 0.90 Alignments for each domain: == domain 1 score: 141.6 bits; conditional E-value: 1.6e-45 Gcd10p 26 ttisLgKfgsfplnliigrpygltfeildkreeeeksrlrvvpaaeleaeslaeeeaeeeeeeeeardnreiiddgarQkLtkeeIe 112 ++++gK + f+++++ig yg +++ ++r + vp+++++ e+ a + dnr ii +a+Q L+ + I+ tr|A9UTS9|A9UTS9_MONBE 24 RHVRMGKIR-FNSKHLIGTAYGEFYRVDRDR------LIPYVPEEDNDDEDAAV-----------EGDNRGIIATKANQGLDDQAIQ 92 6899*****.****************99998......34556554444444433...........469******************* PP Gcd10p 113 eLKkegasagkeiIaklleshtafdqKTaFSqeKYlkrKkkKYlkrftvlpldvsllleyllekkdaqkilelreetlglllslanv 199 +LK eg++ g++++++l+e++++f+qK aF +eKY+krKk+ +++ +++l++++++l + +++ + ++ki +r +tl+ ll anv tr|A9UTS9|A9UTS9_MONBE 93 NLKSEGVT-GEALVQTLVENSKTFKQKHAFTREKYIKRKKEIHVNAVQILRPTIRNLHDAFYQ-HFPDKIASIRVDTLASLLMRANV 177 ********.******************************************************.8********************** PP Gcd10p 200 haggryLvvDdtgGLlvaalaeRmgifessakegtit 236 ++++r++v+++++GL++aa+a ++g +++ ti+ tr|A9UTS9|A9UTS9_MONBE 178 QPHSRVVVFETCSGLMTAAVARQLGA----ESKATIV 210 *************************9....5555555 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (300 nodes) Target sequences: 9188 (5511963 residues searched) Passed MSV filter: 542 (0.05899); expected 183.8 (0.02) Passed bias filter: 241 (0.0262299); expected 183.8 (0.02) Passed Vit filter: 24 (0.0026121); expected 9.2 (0.001) Passed Fwd filter: 1 (0.000108838); expected 0.1 (1e-05) Initial search space (Z): 9188 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.14u 0.00s 00:00:00.14 Elapsed: 00:00:00.05 # Mc/sec: 33071.78 // [ok]