# 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/Pan_troglodytes.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2.7e-71 240.5 0.9 4.7e-71 239.7 0.5 1.6 1 tr|H2QJX7|H2QJX7_PANTR tRNA (adenine(58)-N(1))-methyltr ------ inclusion threshold ------ 0.24 10.8 1.0 0.34 10.3 1.0 1.2 1 tr|A0A2I3RP28|A0A2I3RP28_PANTR Sodium/potassium-transporting AT Domain annotation for each sequence (and alignments): >> tr|H2QJX7|H2QJX7_PANTR tRNA (adenine(58)-N(1))-methyltransferase non-catalytic subunit TRM6 OS=Pan troglodytes OX=95 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 239.7 0.5 4.1e-75 4.7e-71 2 266 .. 20 254 .. 19 322 .. 0.84 Alignments for each domain: == domain 1 score: 239.7 bits; conditional E-value: 4.1e-75 Gcd10p 2 IkpnqhvllkLpsenlkivqvkpnttisLgKfgsfplnliigrpygltfeildkreeeeksrlrvvpaaeleaeslaeeeaeeeeee 88 I+ +++v+lk ++++k vqv++ ++++++K++ f+l+++ig+ yg+ fe+ + + + k+ +++e +ae++ tr|H2QJX7|H2QJX7_PANTR 20 IRDGDFVVLK-REDVFKAVQVQRRKKVTFEKQW-FYLDNVIGHSYGTAFEVTSGGSLQPKK---------KREEPTAETK------- 88 999*******.777*******************.********************9654433.........3333333332....... PP Gcd10p 89 eeardnreiiddgarQkLtkeeIeeLKkegasagkeiIaklleshtafdqKTaFSqeKYlkrKkkKYlkrftvlpldvsllleylle 175 e+++dnr+i+ddg++QkLt+++I++LK++g++ g+ei+++l+e++t+f++KT+F+q+KY+k+KkkKY +tv +++++ l ++++ tr|H2QJX7|H2QJX7_PANTR 89 EAGTDNRNIVDDGKSQKLTQDDIKALKDKGIK-GEEIVQQLIENSTTFRDKTEFAQDKYIKKKKKKYEAIITVVKPSTRILSIMYYA 174 2389****************************.****************************************************** PP Gcd10p 176 kkdaqkilelreetlglllslanvhaggryLvvDdtgGLlvaalaeRmgifessakegtitliheneqpnlsllkyfnydaaepehp 262 ++++ki ++r +tl+++l+l+n++ag++++v+++++GL+++a++eRmg+ g+i++++++ p ++ +f++ p+++ tr|H2QJX7|H2QJX7_PANTR 175 -REPGKINHMRYDTLAQMLTLGNIRAGNKMIVMETCAGLVLGAVMERMGG------FGSIIQLYPGGGPVRAATACFGF----PKSF 250 .9************************************************......***********************....6666 PP Gcd10p 263 lkkh 266 l+ tr|H2QJX7|H2QJX7_PANTR 251 LSGL 254 5543 PP >> tr|A0A2I3RP28|A0A2I3RP28_PANTR Sodium/potassium-transporting ATPase subunit beta OS=Pan troglodytes OX=9598 GN=ATP1B # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.3 1.0 3e-05 0.34 46 90 .. 15 71 .. 4 108 .. 0.70 Alignments for each domain: == domain 1 score: 10.3 bits; conditional E-value: 3e-05 Gcd10p 46 ygltfeildkreeeeks............rlrvvpaaeleaeslaeeeaeeeeeeee 90 y +++ d +e ++++ r+ vvp++e+e+e+ ++ee+eeee+eee tr|A0A2I3RP28|A0A2I3RP28_PANTR 15 YSYRYRLDDPDEANQNYladeeeeaeeeaRVTVVPKSEEEEEEEEKEEEEEEEKEEE 71 55555555555444444555555555556*******999999999888887777655 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (300 nodes) Target sequences: 23008 (11942255 residues searched) Passed MSV filter: 1363 (0.0592403); expected 460.2 (0.02) Passed bias filter: 470 (0.0204277); expected 460.2 (0.02) Passed Vit filter: 40 (0.00173853); expected 23.0 (0.001) Passed Fwd filter: 2 (8.69263e-05); expected 0.2 (1e-05) Initial search space (Z): 23008 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.29u 0.01s 00:00:00.30 Elapsed: 00:00:00.10 # Mc/sec: 35826.76 // [ok]