# 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_344.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_344 [L=236] Description: # 379178 # 379885 # -1 # ID=1_344;partial=00;start_type=ATG;rbs_motif=GGAGG;rbs_spacer=5-10bp;gc_cont=0.326 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 9.7e-53 178.8 0.0 1.2e-52 178.5 0.0 1.1 1 Ribosomal_L1 Ribosomal protein L1p/L10e family Domain annotation for each model (and alignments): >> Ribosomal_L1 Ribosomal protein L1p/L10e family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 178.5 0.0 7.1e-57 1.2e-52 2 205 .] 32 221 .. 30 221 .. 0.98 Alignments for each domain: == domain 1 score: 178.5 bits; conditional E-value: 7.1e-57 HHHC-SSS--EEEEEEEESSS-TTSGGG--EEEEE-STSTSSCS-EEEEE-HHHHH...HHHHTT...-SEEETCHHHHHHT-H CS Ribosomal_L1 2 lkhakkkfdesveialnlttkkpkksdqrkpgkipLphplgkevkiclitkdpqakyakeakeagiriakviglkklkkkyksf 85 k++++kf+e+ve+al+l +++p+++dq+++g+++Lph++gk+vki+ it +++ + +a +ag a++ g ++ ++++ 23726_FusoPortal_Gene_344 32 QKTKTAKFTETVEVALRL-GVDPRHADQQIRGTVVLPHGTGKTVKILAITSGENIE---KALAAG---ADYAGAEEYINQI--- 105 577788************.*************************************...99****...*************... PP HHHHHTCCH-SEEEEEGGGHHHHHHHCHHHHGGGT-S-ECC..CT.BTS-HHHHHHHHCTTEEEEE-TTTSEEEEEEEETTSSH CS Ribosomal_L1 86 eakrkllkefdlfladdrimpllpklLGkvfgkrkkmPipvltkgtvkedvkkevekalkstvvvrldkgtclsvkvGklsfta 169 ++++ +fdl++a++++mp++++ LGk++g++++mP+p+ +gtv++d++++v +++k+++ +r+dk +++++++Gk++f+ 23726_FusoPortal_Gene_344 106 ---QQGWLDFDLVIATPDMMPKIGR-LGKILGTKGLMPNPK--SGTVTPDIAAAVSEFKKGKLAFRVDKLGSIHAPIGKVDFDL 183 ...9*********************.***************..9**************************************** PP HHHHHHHHHHHHHHHHTSTTT..SGSEEEEEEEETTS- CS Ribosomal_L1 170 eelaeNieavvkalvekvpkk..wkniksihlkttmsp 205 +++ eN++a+++++ + +p++ +++++++ ++ tm+p 23726_FusoPortal_Gene_344 184 DKIEENFKAFMDQIIRLKPASskGQYLRTVAVSLTMGP 221 *****************998888*************98 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (236 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 421 (0.0251915); expected 334.2 (0.02) Passed bias filter: 389 (0.0232767); expected 334.2 (0.02) Passed Vit filter: 27 (0.00161561); expected 16.7 (0.001) Passed Fwd filter: 1 (5.98372e-05); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.25u 0.24s 00:00:00.49 Elapsed: 00:00:00.33 # Mc/sec: 2078.97 // [ok]