# 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/Amphimedon_queenslandica.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2e-39 136.8 7.5 2.3e-36 126.7 3.1 2.2 2 tr|A0A1X7UZN4|A0A1X7UZN4_AMPQE Uncharacterized protein OS=Amphi Domain annotation for each sequence (and alignments): >> tr|A0A1X7UZN4|A0A1X7UZN4_AMPQE Uncharacterized protein OS=Amphimedon queenslandica OX=400682 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 126.7 3.1 5.4e-41 2.3e-36 2 182 .. 6 173 .. 5 176 .. 0.92 2 ! 12.4 0.1 3.6e-06 0.16 206 225 .. 176 195 .. 174 253 .. 0.85 Alignments for each domain: == domain 1 score: 126.7 bits; conditional E-value: 5.4e-41 Gcd10p 2 IkpnqhvllkLpsenlkivqvkpnttisLgKfgsfplnliigrpygltfeildkreeeeksrlrvvpaaeleaeslaee 80 I++n++vll ++ +k +++++ +i+++K f +++ig pyg+ fe+ + + l +v+++e++ e+ + + tr|A0A1X7UZN4|A0A1X7UZN4_AMPQE 6 IEENKVVLLH-VKDIYKPIRLVKRRKIKFDKIS-FVPDAAIGLPYGSEFEVARGN-------LILVKEEERRGEEGTIS 75 99********.999******************9.*****************9999.......66887788777776666 PP Gcd10p 81 eaeeeeeeeeardnreiiddgarQkLtkeeIeeLKkegasagkeiIaklleshtafdqKTaFSqeKYlkrKkkKYlkrf 159 + +e e+a+dnr i+ddg +QkL++eeI+++Kk+g+s g+ei + +++ +++f+++T+FSqeKY+k+K k++ + + tr|A0A1X7UZN4|A0A1X7UZN4_AMPQE 76 QLTE--VVESAKDNRLIVDDGLSQKLSHEEIHQMKKDGMS-GQEILESVVQGSSSFKDRTKFSQEKYVKKKYKRHSHVL 151 5444..4457******************************.************************************** PP Gcd10p 160 tvlpldvsllleyllekkdaqki 182 ++ +++++l+++y+++ k ++ki tr|A0A1X7UZN4|A0A1X7UZN4_AMPQE 152 SIIKPTTRLIAQYYYA-KGPSKI 173 ****************.899998 PP == domain 2 score: 12.4 bits; conditional E-value: 3.6e-06 Gcd10p 206 LvvDdtgGLlvaalaeRmgi 225 +v+++++GLl++a++eR+g tr|A0A1X7UZN4|A0A1X7UZN4_AMPQE 176 VVMETCQGLLTGAVLERLGD 195 699***************97 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (300 nodes) Target sequences: 43435 (14716620 residues searched) Passed MSV filter: 2600 (0.0598596); expected 868.7 (0.02) Passed bias filter: 950 (0.0218718); expected 868.7 (0.02) Passed Vit filter: 87 (0.00200299); expected 43.4 (0.001) Passed Fwd filter: 2 (4.60458e-05); expected 0.4 (1e-05) Initial search space (Z): 43435 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.37u 0.01s 00:00:00.38 Elapsed: 00:00:00.12 # Mc/sec: 36791.55 // [ok]