# 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_1954.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1954 [L=359] Description: # 2102121 # 2103197 # -1 # ID=1_1954;partial=00;start_type=ATG;rbs_motif=AGGAG;rbs_spacer=5-10bp;gc_cont=0.323 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 3.4e-54 182.0 0.0 4.9e-54 181.4 0.0 1.2 1 CPSase_sm_chain Carbamoyl-phosphate synthase small chain, CP 6.8e-46 156.3 0.0 9.7e-46 155.8 0.0 1.2 1 GATase Glutamine amidotransferase class-I 1.7e-05 24.6 0.4 0.049 13.3 0.1 2.4 2 Peptidase_C26 Peptidase C26 ------ inclusion threshold ------ 0.038 14.3 0.0 0.092 13.0 0.0 1.7 1 Radical_SAM Radical SAM superfamily Domain annotation for each model (and alignments): >> CPSase_sm_chain Carbamoyl-phosphate synthase small chain, CPSase domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 181.4 0.0 1.2e-57 4.9e-54 3 127 .] 5 131 .. 3 131 .. 0.97 Alignments for each domain: == domain 1 score: 181.4 bits; conditional E-value: 1.2e-57 CPSase_sm_chain 3 vLvLedGtvfeGksfga..eavgEvvfntsmtGYqeiltDpsYagqivvltypliGnyGvneedvesekikvaglvvrelsee 83 L+LedGtv++G++fga e+vgEvvfntsmtGYqeil+DpsY+gqiv+ltypliGnyG+n++d+es k+ ++g+vv+e++++ 23726_FusoPortal_Gene_1954 5 QLILEDGTVYKGYAFGAdvENVGEVVFNTSMTGYQEILSDPSYNGQIVTLTYPLIGNYGINRDDFESMKPCIKGMVVKEVCTT 87 69**************54488************************************************************** PP CPSase_sm_chain 84 psnwraeesLeewlkeegipgisgvDTRaLtkklrekgamkgvi 127 psn+r+e++L+e lk+ gipgi+g+DTRaLt+klr+kg++kg++ 23726_FusoPortal_Gene_1954 88 PSNFRSEKTLDEALKDFGIPGIYGIDTRALTRKLRSKGVVKGCL 131 ****************************************9975 PP >> GATase Glutamine amidotransferase class-I # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 155.8 0.0 2.3e-49 9.7e-46 2 189 .. 175 351 .. 174 352 .. 0.95 Alignments for each domain: == domain 1 score: 155.8 bits; conditional E-value: 2.3e-49 EEESSTSTHHHHHHHHHHTT-EEEEEETTSHHHHHHHTHHSTSEEEEEEESSCTT.TTTHHHHHHHHCTSTS-EEEETHHHHH CS GATase 2 liDnydsvtenlvhalvelgvkvevvpidasaeeilellknpdgilvsgGpgspg.vegkieaikearenkiPiLGICLGmQl 83 l+D + ++ +v++lv++g++v vvp++++aee l + +pdg+++ +Gpg+p+ ++++ie+ik +++ k+ i+GIC+G+Ql 23726_FusoPortal_Gene_1954 175 LVDLG--MKIGIVRELVSRGCDVIVVPYNTTAEEVL--RLEPDGVMLTNGPGDPEdAKESIEMIKGIIN-KVTIFGICMGHQL 252 67777..88899***********************5..79*****************************.************* PP HHHHTTTEEEEEECEEEEEEEEEEEETTTCGGTTTTSEEEEEEEEEEEEESGCHHHTTEEEEEEEETTEEEEEEEEETTTEEE CS GATase 84 laiefgrkvvkekekehhggtmrlgddetallyglknevierhrHryevnpkelepaglevtakseddeelmeiiekknpfvg 166 + ++ g+k++k k ++h+gg++++++ +l+g ++ i++++H+y+v+ ++l +++le t+ +d + +++++kk+p++ 23726_FusoPortal_Gene_1954 253 VSLACGAKTYKLK-FGHRGGNHPVKN----ILTG--RVDITSQNHGYAVDIDSLNDTDLELTHIAINDRSCEGVRHKKYPVFT 328 *************.************....5555..6889***********9999999************************* PP ESSBTTSTTSTTHHHHHHHHHHH CS GATase 167 vQFHPEststpegpellfnflik 189 vQFHPE+ +p+++++lf+++ik 23726_FusoPortal_Gene_1954 329 VQFHPEAAAGPHDTSYLFDEFIK 351 *********************97 PP >> Peptidase_C26 Peptidase C26 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 13.3 0.1 1.2e-05 0.049 94 128 .. 230 263 .. 216 276 .. 0.80 2 ! 9.1 0.0 0.00022 0.92 174 216 .] 290 334 .. 271 334 .. 0.72 Alignments for each domain: == domain 1 score: 13.3 bits; conditional E-value: 1.2e-05 Peptidase_C26 94 elaLiraalergvPilgiCrGlqllnvalGGtlyq 128 ++ +i+ ++ ++v i+giC G+ql+ a G ++y+ 23726_FusoPortal_Gene_1954 230 SIEMIKGII-NKVTIFGICMGHQLVSLACGAKTYK 263 455565555.5788****************99885 PP == domain 2 score: 9.1 bits; conditional E-value: 0.00022 Peptidase_C26 174 Hhqaike...laeelrveAvae.DgviEaieseeakafvlgvQwHPE 216 H a++ +++l+ + +a D E +++++ ++v+ vQ+HPE 23726_FusoPortal_Gene_1954 290 HGYAVDIdslNDTDLELTHIAInDRSCEGVRHKK--YPVFTVQFHPE 334 5555544333456677777776699999999988..78********9 PP >> Radical_SAM Radical SAM superfamily # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.0 0.0 2.2e-05 0.092 75 150 .. 168 247 .. 109 252 .. 0.78 Alignments for each domain: == domain 1 score: 13.0 bits; conditional E-value: 2.2e-05 CCTEEEEEEECSSCH.HHHHHHHHHTT-SEEEEEEBSSSHH......HHHHHHTCTSHHHHHHHHHHHHHHCTSET.EEEEEE CS xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx......xxxxxxxxxxxxxxxxxxxxxxxxxxxxx.xxxxxx RF Radical_SAM 75 vpgiritvetngtlpdeeklerllelgvdrvslgleaadde......kvlklmnrghtfeevlealeklreagikr.vvdriv 150 +g+r++++ g++ ++++l+ +g+d++ ++ +++ +e + l+n+++++e+ +e++e+++ i++ ++++i 23726_FusoPortal_Gene_1954 168 GRGRRVVLVDLGMK--IGIVRELVSRGCDVIVVPYNTTAEEvlrlepDGVMLTNGPGDPEDAKESIEMIKGI-INKvTIFGIC 247 56777777777777..5899*********************999999888899999**************99.7754444444 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (359 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 370 (0.0221398); expected 334.2 (0.02) Passed bias filter: 345 (0.0206438); expected 334.2 (0.02) Passed Vit filter: 30 (0.00179512); expected 16.7 (0.001) Passed Fwd filter: 4 (0.000239349); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 4 [number of targets reported over threshold] # CPU time: 0.26u 0.23s 00:00:00.49 Elapsed: 00:00:00.31 # Mc/sec: 3366.53 // [ok]