# hmmscan :: search sequence(s) against a profile 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 sequence file: AA/23726_FusoPortal_Gene_121.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_121 [L=248] Description: # 140035 # 140778 # 1 # ID=1_121;partial=00;start_type=ATG;rbs_motif=AGGAGG;rbs_spacer=5-10bp;gc_cont=0.304 Scores for complete sequence (score includes all domains): --- full sequence --- --- best 1 domain --- -#dom- E-value score bias E-value score bias exp N Model Description ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2.1e-96 321.2 0.3 2.4e-96 321.0 0.3 1.0 1 Ribosomal_S2 Ribosomal protein S2 ------ inclusion threshold ------ 0.015 15.4 0.1 0.038 14.1 0.0 1.7 2 Egg_lysin Egg lysin (Sperm-lysin) 0.047 13.6 3.0 14 5.6 0.1 3.1 3 DUF4363 Domain of unknown function (DUF4363) 0.065 12.7 4.3 0.36 10.3 4.2 1.9 2 DUF4795 Domain of unknown function (DUF4795) 0.097 12.4 3.6 0.29 10.8 3.6 1.7 1 TACC Transforming acidic coiled-coil-containing prot Domain annotation for each model (and alignments): >> Ribosomal_S2 Ribosomal protein S2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 321.0 0.3 7.3e-100 2.4e-96 1 215 [] 9 225 .. 9 225 .. 0.99 Alignments for each domain: == domain 1 score: 321.0 bits; conditional E-value: 7.3e-100 Ribosomal_S2 1 LLeagvhfGhqtrrwnpkmkpyIygerngihiidlektlaalrraykvvkevaakggkilfvgtkkqakeavkeaakrangyyv 84 LLeagvhfGhq++rwnpkm++yI++erngih+idl+k l+++++ay+ ++++a++ggk+lfvgtkkqa+ea+ke+a+r+++yyv 23726_FusoPortal_Gene_121 9 LLEAGVHFGHQAKRWNPKMAKYIFTERNGIHVIDLHKSLKKIEEAYEEMRKIAEDGGKVLFVGTKKQAQEAIKEQAERSGMYYV 92 8*********************************************************************************** PP Ribosomal_S2 85 nerWlgGtLTNfktikksikklkeleeleedgeld..ltkkealklkkekekleknlggikkmkklPdlvvvldpkkekiavkE 166 n+rWlgG+LTNf+tikk+i+++kele+l+++g ld +tkkea++++ke +kl+knl+gi++m+k+Pd++ v+d+k e++ vkE 23726_FusoPortal_Gene_121 93 NSRWLGGMLTNFSTIKKRIERMKELEKLDAEGILDtdYTKKEAAEFRKELSKLSKNLSGIRDMEKVPDAIYVVDVKMEELPVKE 176 ********************************99778*********************************************** PP Ribosomal_S2 167 akklgipviaivDtnvdpdlvdypIpaNddsirsiklilgvladaieeg 215 a+ lgipv+a++Dtnvdpdl++ypIpaNdd+irs+kli++v+a+ai eg 23726_FusoPortal_Gene_121 177 AHLLGIPVFAMIDTNVDPDLITYPIPANDDAIRSVKLITSVIANAIVEG 225 **********************************************997 PP >> Egg_lysin Egg lysin (Sperm-lysin) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -3.8 0.0 4.2 1.4e+04 103 113 .. 22 32 .. 19 38 .. 0.72 2 ? 14.1 0.0 1.1e-05 0.038 33 68 .. 82 116 .. 69 124 .. 0.82 Alignments for each domain: == domain 1 score: -3.8 bits; conditional E-value: 4.2 ---HHHHHHHC CS Egg_lysin 103 kiraymeklla 113 ++++ m+k++ 23726_FusoPortal_Gene_121 22 RWNPKMAKYIF 32 68888888765 PP == domain 2 score: 14.1 bits; conditional E-value: 1.1e-05 -HHHHHHHHHHHHHHHHHCHHHHHHHHHHHHHCCTS CS Egg_lysin 33 skvqkktlrfvnrrrmqtlwqnYaewvskkirkLgr 68 ++ ++ +++vn+r++ + n++++ +k+i+++++ 23726_FusoPortal_Gene_121 82 EQAERSGMYYVNSRWLGGMLTNFSTI-KKRIERMKE 116 6778889*******************.*****9973 PP >> DUF4363 Domain of unknown function (DUF4363) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 5.3 0.1 0.0054 18 59 107 .. 33 81 .. 16 87 .. 0.80 2 ? 1.7 0.0 0.071 2.4e+02 96 120 .] 101 125 .. 90 125 .. 0.83 3 ? 5.6 0.1 0.0043 14 88 115 .. 131 158 .. 124 159 .. 0.85 Alignments for each domain: == domain 1 score: 5.3 bits; conditional E-value: 0.0054 DUF4363 59 kkkkklsifiehdeldeietslaklkayietedkeealaelaelkflle 107 +++ +++ h l++ie+++ ++++ +e + k ++++ +++++ ++ 23726_FusoPortal_Gene_121 33 TERNGIHVIDLHKSLKKIEEAYEEMRKIAEDGGKVLFVGTKKQAQEAIK 81 56777888899*******************9999999988888777665 PP == domain 2 score: 1.7 bits; conditional E-value: 0.071 DUF4363 96 laelaelkfllehikekeklsleNI 120 l + +++k+ +e +ke ekl e I 23726_FusoPortal_Gene_121 101 LTNFSTIKKRIERMKELEKLDAEGI 125 6778899999999999999987765 PP == domain 3 score: 5.6 bits; conditional E-value: 0.0043 DUF4363 88 etedkeealaelaelkfllehikekekl 115 +++ e+ el++l + l i+++ek+ 23726_FusoPortal_Gene_121 131 TKKEAAEFRKELSKLSKNLSGIRDMEKV 158 56888999****************9985 PP >> DUF4795 Domain of unknown function (DUF4795) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 1.9 0.4 0.041 1.4e+02 16 59 .. 40 84 .. 6 94 .. 0.73 2 ? 10.3 4.2 0.00011 0.36 111 185 .. 57 158 .. 42 178 .. 0.58 Alignments for each domain: == domain 1 score: 1.9 bits; conditional E-value: 0.041 DUF4795 16 ieklqeelekleekvkslekrledkekklea.lkkqieeLekvka 59 + +l ++l+k+ee+ +++ k ed k l++ kkq++e k+ a 23726_FusoPortal_Gene_121 40 VIDLHKSLKKIEEAYEEMRKIAEDGGKVLFVgTKKQAQEAIKEQA 84 667777777777777777777777666665414445555555544 PP == domain 2 score: 10.3 bits; conditional E-value: 0.00011 DUF4795 111 lkkaldtLqeelkekldrlelenLkeqle..........................erlkelqkelqeamekedk.edeaaglrk 167 + k+ ++ + l+ ++ e +keq e er+kel+k e + +d ++eaa++rk 23726_FusoPortal_Gene_121 57 MRKIAEDGGKVLFVGTKKQAQEAIKEQAErsgmyyvnsrwlggmltnfstikkriERMKELEKLDAEGILDTDYtKKEAAEFRK 140 44444444444444444444444444444444444444444444444444444445999999766666655555799******* PP DUF4795 168 kLlerfdcISCdrpvemk 185 +L++ +++S r++e 23726_FusoPortal_Gene_121 141 ELSKLSKNLSGIRDMEKV 158 ***999999977777655 PP >> TACC Transforming acidic coiled-coil-containing protein (TACC) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.8 3.6 8.6e-05 0.29 22 104 .. 38 117 .. 26 148 .. 0.76 Alignments for each domain: == domain 1 score: 10.8 bits; conditional E-value: 8.6e-05 TACC 22 eeiqelkkkveelkkenqemkkivaefektiaqlieekqkqkelakkeiqkvleekeqlladlnsveksfsdlfkryeklkev 104 ++ +l+k+ ++ ++ ++em+ki ++ k++ ++++ +++ ++++e + + ++ l + + fs + kr+e++ke 23726_FusoPortal_Gene_121 38 IHVIDLHKSLKKIEEAYEEMRKIAEDGGKVLFVGTKKQAQEAIKEQAERSGMY---YVNSRWLGGMLTNFSTIKKRIERMKEL 117 46889**********************99877666666555555555555554...3445557777899************94 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (248 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 1125 (0.0673169); expected 334.2 (0.02) Passed bias filter: 692 (0.0414074); expected 334.2 (0.02) Passed Vit filter: 76 (0.00454763); expected 16.7 (0.001) Passed Fwd filter: 5 (0.000299186); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 5 [number of targets reported over threshold] # CPU time: 0.27u 0.23s 00:00:00.50 Elapsed: 00:00:00.31 # Mc/sec: 2325.63 // [ok]