# 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/Dictyostelium_discoideum.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1e-63 214.7 11.8 1e-63 214.7 11.8 1.8 2 sp|Q54UB1|TRM6_DICDI tRNA (adenine(58)-N(1))-methyltransfera Domain annotation for each sequence (and alignments): >> sp|Q54UB1|TRM6_DICDI tRNA (adenine(58)-N(1))-methyltransferase non-catalytic subunit trm6 OS=Dictyostelium discoideu # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 214.7 11.8 8.2e-68 1e-63 1 285 [. 34 320 .. 34 336 .. 0.81 2 ? -2.6 0.2 0.12 1.6e+03 106 106 .. 479 479 .. 403 519 .. 0.58 Alignments for each domain: == domain 1 score: 214.7 bits; conditional E-value: 8.2e-68 Gcd10p 1 lIkpnqhvllkLpse.nlkivqvkpnttisLgKfgsfplnliigrpygltfeildkreeeeksrlrvvpaaeleaeslaeeeaeeeeee 88 +Ik+++hv+l+++++ ++ +++vk ++++++gK++ + +n+iig y++ f++ +++++ l+ + ++e +++ ++ e + sp|Q54UB1|TRM6_DICDI 34 IIKEGDHVILDINNGeKFSVIKVKLGSKVKIGKKQ-ILINSIIGESYYSSFQVSNEKNT-----LERITQKEI-----DDRLNNLVELN 111 59************99*******************.******************99954.....444533333.....33333333334 PP Gcd10p 89 eeardnreiiddgarQkLtkeeIeeLKkegasagkeiIaklleshtafdqKTaFSqeKYlkrKkkKYlkrftvlpldvsllleyllekk 177 ++ +dnr++ ++++QkLt+e+I+e+K++g++ ++ iI++++e++++f++KT+FSq KYlk+K kKY + +++ +++ ++l+e +++ k sp|Q54UB1|TRM6_DICDI 112 QNDADNRNLDQNNTAQKLTQEDINEMKQKGTD-SNTIIKTIVENSESFKTKTSFSQIKYLKKKIKKYSTIVKIIKPTLKSLTEAYYK-K 198 4579****************************.******************************************************.9 PP Gcd10p 178 daqkilelreetlglllslanvhaggryLvvDdtgGLlvaalaeRmgifessakegtitliheneqpnlsllkyfnydaaepehplkkh 266 d ki+ lr +++g+ll+l+n++a++++Lvv++++GL+++++aeRm++ +gti + ++ p+ls++++f+++++ +++ + sp|Q54UB1|TRM6_DICDI 199 DSRKICGLRFDSFGQLLTLGNIRANSQVLVVETCMGLVTGSIAERMNG------QGTILSAYIGKGPSLSIVNNFGFNTDVLNTIYPFN 281 ************************************************......***********************443333333334 PP Gcd10p 267 lktlswlqlle....................peedetye 285 l++ s+l+ e +e++++++ sp|Q54UB1|TRM6_DICDI 282 LNITSVLNKGEdisklpptassgiydkqvkeKEKEKEKD 320 777777733334555555444443333221133333333 PP == domain 2 score: -2.6 bits; conditional E-value: 0.12 Gcd10p 106 L 106 L sp|Q54UB1|TRM6_DICDI 479 L 479 0 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (300 nodes) Target sequences: 12739 (6856946 residues searched) Passed MSV filter: 4129 (0.324123); expected 254.8 (0.02) Passed bias filter: 249 (0.0195463); expected 254.8 (0.02) Passed Vit filter: 38 (0.00298297); expected 12.7 (0.001) Passed Fwd filter: 1 (7.84991e-05); expected 0.1 (1e-05) Initial search space (Z): 12739 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.30u 0.01s 00:00:00.31 Elapsed: 00:00:00.09 # Mc/sec: 22856.49 // [ok]