# 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_1844.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1844 [L=342] Description: # 1974417 # 1975442 # -1 # ID=1_1844;partial=00;start_type=ATG;rbs_motif=AGGAGG;rbs_spacer=5-10bp;gc_cont=0.298 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2e-98 329.3 0.1 2.3e-98 329.1 0.1 1.0 1 Peptidase_M22 Glycoprotease family 1.3e-07 31.5 0.2 4.4e-05 23.2 0.0 2.2 2 Carbam_trans_N Carbamoyltransferase N-terminus ------ inclusion threshold ------ 0.012 15.0 0.1 1 8.7 0.0 2.2 2 BcrAD_BadFG BadF/BadG/BcrA/BcrD ATPase family Domain annotation for each model (and alignments): >> Peptidase_M22 Glycoprotease family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 329.1 0.1 4.1e-102 2.3e-98 1 270 [. 23 305 .. 23 306 .. 0.97 Alignments for each domain: == domain 1 score: 329.1 bits; conditional E-value: 4.1e-102 XCEXEEEEEEEEE....SSSSHHHHHHHHHHHHHHHHHHHHHHHTSTGGGEEEEEEEEESS-HHHHHHHHHHHHHHHHHCT-- CS Peptidase_M22 1 evlanvvasqkelhakyGGvvPelasraHserllplieealaeaglsledldaIavtkgPGlftgLrvGvsfAkglalalnkP 83 e+l+n+++sq+e+h++yGGvvPe+asr+H +++ ++ee l+ea+++l+d+d Iavt PGl+++L+vGvsfAkgl++a n+P 23726_FusoPortal_Gene_1844 23 EILSNNISSQIEIHKEYGGVVPEIASRQHIKNIATVLEESLEEAKITLDDVDYIAVTYAPGLIGALLVGVSFAKGLSYAKNIP 105 69********************************************************************************* PP EEEEEHHHHHHHCCHCTSECSSC.EEEEESSSTEEEEEECT.TEEEEEEEESSSEHHHHHHHHHHHHT..SSCCHHHHHHHCT CS Peptidase_M22 84 lvgVnhleahalaarleekeefP.vvllvsGghtqvylvke.gryeilgetlddaagealdkvarllglpypggpkieklake 164 +++V+h+++h++a++le++ e P ++l+vsGght++++++e +++ ++getldda+ge++dkvar+lgl+ypggp+i+k++ + 23726_FusoPortal_Gene_1844 106 IIPVHHIKGHMYANFLEHDVELPcISLVVSGGHTNIIYIDEnHNFINIGETLDDAVGESCDKVARVLGLGYPGGPVIDKMYYK 188 *****************7779************************************************************** PP C..S.......EET.TEEE.CHHHHHHHHHHST.T.......HHHHHHHHHHHHHHHHHHHHHHHHHHCTBSEEEEESGGGGX CS Peptidase_M22 165 g..kkiefprplkg.ldfsfsGlktavlkelek.e.......vedlaasfqeavvdhlvekteralklpkakelvilGGvaan 236 g + +++++p+++ +dfsfsG+kta+++ ++ + +edlaasf +vvd l++kt +a+ ++++k+++++GGvaan 23726_FusoPortal_Gene_1844 189 GdrDFLKITKPKVSrFDFSFSGIKTAIINFDNNmKmknqeykKEDLAASFLGTVVDILCDKTLDAAVEKNVKTIMLAGGVAAN 271 *85559999999999**************8855545678999***************************************** PP HHHHHCHCCHHHHTTSEEEXXXHHHHSSXCHHHH CS Peptidase_M22 237 krLrealeeeaeelglkliiPplkyctdnaamia 270 ++Lr++l+e a+e g+k+i+P++k+ctdnaamia 23726_FusoPortal_Gene_1844 272 SLLRSQLTEKAAEKGIKVIYPSMKLCTDNAAMIA 305 ********99999********************8 PP >> Carbam_trans_N Carbamoyltransferase N-terminus # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 6.5 0.1 0.00097 5.4 2 74 .. 3 85 .. 2 126 .. 0.76 2 ! 23.2 0.0 7.8e-09 4.4e-05 244 317 .. 212 283 .. 151 295 .. 0.78 Alignments for each domain: == domain 1 score: 6.5 bits; conditional E-value: 0.00097 Carbam_trans_N 2 ilGikltlflHDgavalvddgklvaaveeE.rltrikrdqk.......lpklkaidyale....eagvelsdiDqivfdeidgei 74 ilGi + ++a+v+dgk++ + + ++ +k +++ +++k+i+ +le ea + l d+D+i++ +g i 23726_FusoPortal_Gene_1844 3 ILGIESS--CDETSIAVVKDGKEILSNNISsQIEIHKEYGGvvpeiasRQHIKNIATVLEesleEAKITLDDVDYIAVTYAPGLI 85 7888866..78899*********998765414444455555888888778888898887511116789999*****988777665 PP == domain 2 score: 23.2 bits; conditional E-value: 7.8e-09 Carbam_trans_N 244 kkkvhdffeasalrlkearekdedlaAsiqvflErlvvevlnivyeklkirnLclaGGvalNiklNsallesgl 317 k ++ +f ++++++++e++++ dlaAs+ + +++ + + + + +++ + laGGva+N l s+l e++ 23726_FusoPortal_Gene_1844 212 KTAIINFDNNMKMKNQEYKKE--DLAASFLGTVVDILCDKTLDAAVEKNVKTIMLAGGVAANSLLRSQLTEKAA 283 344556666667777777777..***************9988888889*********************99865 PP >> BcrAD_BadFG BadF/BadG/BcrA/BcrD ATPase family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 8.7 0.0 0.00018 1 9 70 .. 12 81 .. 4 153 .. 0.80 2 ? 4.4 0.0 0.0037 20 217 251 .. 257 289 .. 220 297 .. 0.64 Alignments for each domain: == domain 1 score: 8.7 bits; conditional E-value: 0.00018 EEEEEEEETTTEEEEEEEE--........TTS-HHH.........HHHHHHHHHHT--GGGECEE..EEE CS BcrAD_BadFG 9 ttkavlmdekgevlgrgitns........anlesvrfeaaekvleealtealerageeeadieetylvgt 70 t +++++ +e+l + i+ + + ++ k++ ++l+e le+a+ + +d+++++++ + 23726_FusoPortal_Gene_1844 12 ETSIAVVKDGKEILSNNISSQieihkeygGVVPEIASRQHIKNIATVLEESLEEAKITLDDVDYIAVTYA 81 57778888888888888888888887777666667778899**********************9987754 PP == domain 2 score: 4.4 bits; conditional E-value: 0.0037 SS...EEEESGGGG-HHHHHHH.....HHHH-S-EE- CS BcrAD_BadFG 217 tde..vtltGGvakneaavkelrklvkenygevkini 251 ++ ++l GGva n+ lr + e+ e+ i++ 23726_FusoPortal_Gene_1844 257 KNVktIMLAGGVAANSL----LRSQLTEKAAEKGIKV 289 22234**********44....5555666666665555 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (342 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 556 (0.0332695); expected 334.2 (0.02) Passed bias filter: 400 (0.0239349); expected 334.2 (0.02) Passed Vit filter: 38 (0.00227382); expected 16.7 (0.001) Passed Fwd filter: 3 (0.000179512); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 3 [number of targets reported over threshold] # CPU time: 0.25u 0.22s 00:00:00.47 Elapsed: 00:00:00.31 # Mc/sec: 3207.11 // [ok]