# 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/Macaca_mulatta.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-72 244.9 1.0 1.6e-72 244.4 0.5 1.5 1 tr|F6RWA5|F6RWA5_MACMU tRNA (adenine(58)-N(1))-methyltr ------ inclusion threshold ------ 0.24 10.7 1.0 0.34 10.2 1.0 1.2 1 tr|A0A1D5QQR7|A0A1D5QQR7_MACMU Sodium/potassium-transporting AT Domain annotation for each sequence (and alignments): >> tr|F6RWA5|F6RWA5_MACMU tRNA (adenine(58)-N(1))-methyltransferase non-catalytic subunit TRM6 OS=Macaca mulatta OX=954 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 244.4 0.5 1.6e-76 1.6e-72 2 266 .. 20 254 .. 19 322 .. 0.85 Alignments for each domain: == domain 1 score: 244.4 bits; conditional E-value: 1.6e-76 Gcd10p 2 IkpnqhvllkLpsenlkivqvkpnttisLgKfgsfplnliigrpygltfeildkreeeeksrlrvvpaaeleaeslaeeeaeeeeee 88 I+ +++v+lk ++++k vqv++ ++++++K++ f+l+++ig+ yg+tfe+ + + + k+ +++e +ae++ tr|F6RWA5|F6RWA5_MACMU 20 IRDGDFVVLK-REDVFKAVQVQRRKKVTFEKQW-FYLDNVIGHSYGTTFEVTSGGSLQPKK---------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|F6RWA5|F6RWA5_MACMU 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|F6RWA5|F6RWA5_MACMU 175 -REPGKINHMRYDTLAQMLTLGNIRAGNKMIVMETCAGLVLGAMMERMGG------FGSIIQLYPGGGPVRAATACFGF----PKSF 250 .9************************************************......***********************....7777 PP Gcd10p 263 lkkh 266 l+ tr|F6RWA5|F6RWA5_MACMU 251 LNGL 254 6643 PP >> tr|A0A1D5QQR7|A0A1D5QQR7_MACMU Sodium/potassium-transporting ATPase subunit beta OS=Macaca mulatta OX=9544 GN=ATP1B4 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.2 1.0 3.2e-05 0.34 50 90 .. 31 71 .. 5 105 .. 0.69 Alignments for each domain: == domain 1 score: 10.2 bits; conditional E-value: 3.2e-05 Gcd10p 50 feildkreeeeksrlrvvpaaeleaeslaeeeaeeeeeeee 90 + +++e ee+ r+ vvp++e+e+e+ ++ee+eeee+eee tr|A0A1D5QQR7|A0A1D5QQR7_MACMU 31 YLADEEEEAEEEARVTVVPKSEEEEEEEEKEEEEEEEKEEE 71 3344455555555*******999999999888887777655 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (300 nodes) Target sequences: 21211 (10888574 residues searched) Passed MSV filter: 1260 (0.0594031); expected 424.2 (0.02) Passed bias filter: 433 (0.0204139); expected 424.2 (0.02) Passed Vit filter: 35 (0.00165009); expected 21.2 (0.001) Passed Fwd filter: 2 (9.42907e-05); expected 0.2 (1e-05) Initial search space (Z): 21211 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.27u 0.01s 00:00:00.28 Elapsed: 00:00:00.09 # Mc/sec: 36295.25 // [ok]