# 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_1580.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1580 [L=456] Description: # 1694125 # 1695492 # -1 # ID=1_1580;partial=00;start_type=ATG;rbs_motif=GGAGG;rbs_spacer=5-10bp;gc_cont=0.264 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.6e-15 56.5 22.0 2.1e-09 37.7 2.8 2.7 2 Methylase_S Type I restriction modification DNA specificity Domain annotation for each model (and alignments): >> Methylase_S Type I restriction modification DNA specificity domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 37.7 2.8 1.3e-13 2.1e-09 1 167 [. 19 209 .. 19 210 .. 0.67 2 ! 25.0 10.2 1e-09 1.7e-05 6 156 .. 228 394 .. 223 406 .. 0.75 Alignments for each domain: == domain 1 score: 37.7 bits; conditional E-value: 1.3e-13 ETT-EEEEHHHHEEE.EE-----TT.....GGTTT--EEEE-HHHCCTEE-HHHHCCTT---B-.TT-EEEE-S.SSTT-EEE CS Methylase_S 1 pedwevvkLgdvfei.kkgktpkkk.....elkengkypyiragilnnkvvkdglgykkkyilp.kgdilissr.atiGkvay 75 p+ we++k+g +f+i k++++p k +g+ y +++ + + + +++ ++y+l gdil++s + +G+v + 23726_FusoPortal_Gene_1580 19 PSHWETKKIGKIFDIrKEKNSPVKTkevlsLSSMYGVSLYSERKEKGG---NKPKENLEAYNLCyPGDILVNSMnIVAGSVGI 98 789************867777766545554444455555555554444...57777777777666677777766666666666 PP ESS-B.EE-TTEEEEEESTTT--HHHHHHHHHHCHHHHHHH................CCSSTTT---HHHHCC-EEEE--HHH CS Methylase_S 76 rdkky.ylnqsvivlipknkkkineflylylllklkklksl................svgstqpkltkenlkkikIplPplee 141 +++ + + +++l+ ++kk +++++ yl+ +++ + sl ++ + + +++ ++lk++++p Pp ee 23726_FusoPortal_Gene_1580 99 SNY-FgAISPVYYSLQNLSEKKYSKYYLEYLFRNYNFQRSLvglgkgiqmsetedgrLFTVRM-RISWDTLKSQEFPTPPIEE 179 666.22444444444444444456788888887777444445555555556678885555555.******************* PP HHHHHHHHH....HHHHHHHHHHHHHHHHH CS Methylase_S 142 QkkIakiLs....qldntielekkikneie 167 Q +Ia +L+ +d++i +ek++++e+e 23726_FusoPortal_Gene_1580 180 QIQIANYLDwkinEIDRLILIEKEQIKELE 209 *******98555568888889998888888 PP == domain 2 score: 25.0 bits; conditional E-value: 1e-09 EEEHHHHEEEEE-----TT..GGTTT--EEEE-HHHCCT......EE-HHHHCCTT---B-TT-EEEE-S....SSTT-.EEE CS Methylase_S 6 vvkLgdvfeikkgktpkkk..elkengkypyiragilnn......kvvkdglgykkkyilpkgdilissr....atiGk.vay 75 +++L+ ++++ kg ++++ + ++g + y + +++n + ++d y + + ++++ +l++ + + iGk vay 23726_FusoPortal_Gene_1580 228 FISLSKIGTFFKGGGFSREnlSDSDYGAILYGDIYTKYNyffeecISKIDENAYFNSKCIDGNVVLFTGSgetkEDIGKnVAY 310 56666777774555555544778888888555555666666667***********************9889999999965777 PP ESS.-BEE-TTEEEEEESTTT--.HHHHHHHHH.HCHHHHHHH.CCSSTTT---HHHHCC-EEEE--HHHHHHHHHHHHHHHH CS Methylase_S 76 rdk.kyylnqsvivlipknkkki.neflylyll.lklkklksl.svgstqpkltkenlkkikIplPpleeQkkIakiLsqldn 154 + k+ l+ ++i l+p+ k+ +f++++ + k++k+ s g + ++t+ +k+ikIp+ + eeQk I+k ++ ++ 23726_FusoPortal_Gene_1580 311 VGTkKIALGGDIIALKPN--KNFsPKFIAYFSNtSNIKAFKHMkSTGDIIVHITLGAIKSIKIPFISIEEQKDIVKKIDEYIL 391 7663689**********5..4442577777776666666666645666689****************************9987 PP .HH CS Methylase_S 155 .ti 156 + 23726_FusoPortal_Gene_1580 392 nLK 394 333 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (456 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 872 (0.0521781); expected 334.2 (0.02) Passed bias filter: 387 (0.023157); expected 334.2 (0.02) Passed Vit filter: 27 (0.00161561); expected 16.7 (0.001) Passed Fwd filter: 2 (0.000119674); 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.30u 0.24s 00:00:00.54 Elapsed: 00:00:00.33 # Mc/sec: 4016.99 // [ok]