# 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/Ciona_savignyi.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 4.8e-62 209.1 0.1 1.7e-61 207.4 0.1 1.8 1 tr|H2YBX2|H2YBX2_CIOSA Uncharacterized protein OS=Ciona savi 0.0037 15.8 0.0 0.0053 15.3 0.0 1.2 1 tr|H2Z0U5|H2Z0U5_CIOSA tRNA (adenine(58)-N(1))-methyltransfe Domain annotation for each sequence (and alignments): >> tr|H2YBX2|H2YBX2_CIOSA Uncharacterized protein OS=Ciona savignyi OX=51511 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 207.4 0.1 2.9e-65 1.7e-61 2 289 .. 9 265 .. 8 283 .. 0.89 Alignments for each domain: == domain 1 score: 207.4 bits; conditional E-value: 2.9e-65 Gcd10p 2 IkpnqhvllkLpsenlkivqvkpnttisLgKfgsfplnliigrpygltfeildkreeeeksrlrvvpaaeleaeslaeeeaeeeeee 88 I+ n++v+l+ ++ lk v++ pn+ i ++K + f+ +++ g+++g++f+ + + ++ +sl+++++++ tr|H2YBX2|H2YBX2_CIOSA 9 ISLNKTVILD-KNGILKAVKLIPNKMIVFEKLR-FKPDKAEGMNFGTIFKGGYIT---------------TDVDSLDTHADDA---- 74 889*******.888*******************.***************976666...............4567777776554.... PP Gcd10p 89 eeardnr.eiid.dgarQkLtkeeIeeLKkegasagkeiIaklleshtafdqKTaFSqeKYlkrKkkKYlkrftvlpldvsllleyl 173 +++dnr ++id +g++Q++++e+I e+K+eg++ g+ei+++l+e++++f +KTaFSq+KYlk+KkkKYl++++v++++++l+++ + tr|H2YBX2|H2YBX2_CIOSA 75 -QGKDNRlSLIDtHGKSQSIKRENIIEMKDEGFT-GEEIVQRLVENSASFADKTAFSQNKYLKKKKKKYLSYVRVHKPNARLITQLY 159 .378999789997788******************.**************************************************** PP Gcd10p 174 lekkdaqkilelreetlglllslanvhaggryLvvDdtgGLlvaalaeRmgifessakegtitliheneq.pnlsllkyfnydaaep 259 e ++++ki+++r +t++++l+ an+++g ++ vv++++G++++a++eR+g+ +g+++ +h++e+ pnl ++ +fn+ p tr|H2YBX2|H2YBX2_CIOSA 160 QE-REPAKICNIRIDTISQILTAANIRHGANVAVVETCSGVVLGAVLERLGG------NGCVVNLHNGEApPNLFAVSHFNF----P 235 **.9************************************************......************9***********....9 PP Gcd10p 260 ehplkkhlktlswlqllepeedetyeeepe 289 +++++++l+++ + +l ++++ e+e++ tr|H2YBX2|H2YBX2_CIOSA 236 DSYWDDNLHSVPLHKLGPLLAKRSFEDEFP 265 *******99998877655444444444444 PP >> tr|H2Z0U5|H2Z0U5_CIOSA tRNA (adenine(58)-N(1))-methyltransferase catalytic subunit TRMT61A OS=Ciona savignyi OX=5151 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 15.3 0.0 9.2e-07 0.0053 2 53 .. 10 60 .. 9 64 .. 0.92 Alignments for each domain: == domain 1 score: 15.3 bits; conditional E-value: 9.2e-07 Gcd10p 2 IkpnqhvllkLpsenlkivqvkpnttisLgKfgsfplnliigrpygltfeil 53 ++++++v++ + ++++ ++vk+++ is Kfg+++ +++ig++yg+ f tr|H2Z0U5|H2Z0U5_CIOSA 10 VEAGDVVIIFMGHDSMFPINVKKGE-ISQTKFGAIKHDQVIGHEYGTKFTCS 60 78999***************99877.9**********************865 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (300 nodes) Target sequences: 11592 (4231397 residues searched) Passed MSV filter: 694 (0.0598689); expected 231.8 (0.02) Passed bias filter: 214 (0.018461); expected 231.8 (0.02) Passed Vit filter: 16 (0.00138026); expected 11.6 (0.001) Passed Fwd filter: 2 (0.000172533); expected 0.1 (1e-05) Initial search space (Z): 11592 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.11u 0.00s 00:00:00.11 Elapsed: 00:00:00.03 # Mc/sec: 42313.97 // [ok]