# 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/Anolis_carolinensis.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 7.5e-72 242.0 0.1 7.5e-72 242.0 0.1 1.7 2 tr|G1K8L7|G1K8L7_ANOCA tRNA (adenine(58)-N(1))-methyltransfe Domain annotation for each sequence (and alignments): >> tr|G1K8L7|G1K8L7_ANOCA tRNA (adenine(58)-N(1))-methyltransferase non-catalytic subunit TRM6 OS=Anolis carolinensis O # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 242.0 0.1 4e-76 7.5e-72 1 264 [. 8 243 .. 8 313 .. 0.90 2 ? -2.7 1.1 0.14 2.6e+03 139 156 .. 329 346 .. 269 353 .. 0.66 Alignments for each domain: == domain 1 score: 242.0 bits; conditional E-value: 4e-76 Gcd10p 1 lIkpnqhvllkLpsenlkivqvkpnttisLgKfgsfplnliigrpygltfeildkreeeeksrlrvvpaaeleaeslaeeeaeeeee 87 +I+++++ +lk +e++k v v++ +i ++K++ f+l+++ig+ yg+tfe+ + + k+ + e++a+e++ e tr|G1K8L7|G1K8L7_ANOCA 8 VIREGDWAVLK-RNEVFKAVPVARRRKIIFEKQW-FYLDNAIGQIYGTTFEVTGGGNLQPKK----------RVEETATEAK----E 78 699********.888*******************.*******************99765554..........2333333333....2 PP Gcd10p 88 eeeardnreiiddgarQkLtkeeIeeLKkegasagkeiIaklleshtafdqKTaFSqeKYlkrKkkKYlkrftvlpldvsllleyll 174 +++ +dnr+i+ddg++QkLt+++I++LK++g++ g+ei+++l+e++t+f++KT+F+q+KY+k+KkkKY +t+ +++++ l ++++ tr|G1K8L7|G1K8L7_ANOCA 79 AGTDTDNRNIVDDGKSQKLTHDDIKALKDKGIK-GQEIVQQLIENSTTFRDKTEFAQDKYIKKKKKKYEAIITIVKPSTRILSTMYY 164 34479****************************.***************************************************** PP Gcd10p 175 ekkdaqkilelreetlglllslanvhaggryLvvDdtgGLlvaalaeRmgifessakegtitliheneqpnlsllkyfnydaaepeh 261 + ++++ki lr +tl+++l+++nvhag++++v+++++GL+++a++eRmg+ g+i+ ++++ p ++ ++f++ p++ tr|G1K8L7|G1K8L7_ANOCA 165 A-REPGKINYLRYDTLAQMLTWGNVHAGNKMIVMETCAGLVLGAVMERMGG------YGSIINMYPGGGPVRAATNCFGF----PDS 240 *.9************************************************......***********************....555 PP Gcd10p 262 plk 264 + + tr|G1K8L7|G1K8L7_ANOCA 241 FFQ 243 444 PP == domain 2 score: -2.7 bits; conditional E-value: 0.14 Gcd10p 139 KTaFSqeKYlkrKkkKYl 156 Ta+ +eK K+kK++ tr|G1K8L7|G1K8L7_ANOCA 329 GTAHDKEKKQEEKRKKFV 346 355555555555555555 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (300 nodes) Target sequences: 18528 (9796552 residues searched) Passed MSV filter: 1267 (0.068383); expected 370.6 (0.02) Passed bias filter: 409 (0.0220747); expected 370.6 (0.02) Passed Vit filter: 28 (0.00151123); expected 18.5 (0.001) Passed Fwd filter: 1 (5.39724e-05); expected 0.2 (1e-05) Initial search space (Z): 18528 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.24u 0.01s 00:00:00.25 Elapsed: 00:00:00.08 # Mc/sec: 36737.07 // [ok]