# hmmsearch :: search profile(s) against a sequence 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 HMM file: pfam_list/SpoU_methylase.hmm.txt # target sequence database: proteomes/Drosophila_melanogaster.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: SpoU_methylase [M=142] Accession: PF00588.18 Description: SpoU rRNA Methylase family Scores for complete sequences (score includes all domains): --- full sequence --- --- best 1 domain --- -#dom- E-value score bias E-value score bias exp N Sequence Description ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 3.4e-25 88.5 0.0 8e-25 87.3 0.0 1.7 1 tr|Q9VD22|Q9VD22_DROME LP02352p OS=Drosophila melanogaster O 5.1e-21 75.0 0.0 8.4e-21 74.3 0.0 1.4 1 sp|Q9VW14|MRM3_DROME rRNA methyltransferase 3, mitochondri Domain annotation for each sequence (and alignments): >> tr|Q9VD22|Q9VD22_DROME LP02352p OS=Drosophila melanogaster OX=7227 GN=Dmel\CG18596 PE=1 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 87.3 0.0 1.2e-28 8e-25 2 141 .. 986 1127 .. 985 1128 .. 0.92 Alignments for each domain: == domain 1 score: 87.3 bits; conditional E-value: 1.2e-28 EEEEEES-SSHHHHHHHHHHHHHTT-SEEEEESGSSCTCGSHHHHHHCTTHHHHHHEEEE..SSHHHHHHHHCCCTEEEEEEST. CS SpoU_methylase 2 livvlveiedpgNlGaiaRtakafgaselvlvkprsadllnpeakraaagalevlelkvv..dsleealkelkeegfalvatsl. 83 livv i++ NlG +aRt++++g+++l+l + + ++++ ++++ a+++l++ +v +sl+ +l e + eg+++v++++ tr|Q9VD22|Q9VD22_DROME 986 LIVVASLIDKLPNLGGLARTCEVLGVNTLILGLKSQ--AEKSDFTNLSMTAEKTLNILEVnpESLAGFLLEKQMEGYKIVGAEQt 1068 79999*************************965444..69999*************77776789******************996 PP .T-SEGCCSHHHTTEEEEEEBTTTBS-HHHHCTSSEEEE----TTSSS--HHHHHHHHH CS SpoU_methylase 84 .kgkslseidlkekvalvfGnEdkGlsnelleladrlvkIpmsgeveSLNvavAvaill 141 +++++ ++++++k +l++G E++G++ +l+ +d+ v+Ip++g v+SLNv+vA ++++ tr|Q9VD22|Q9VD22_DROME 1069 aHSTNFVDFKFPKKSILLLGHEKHGIPANLIGFLDYAVEIPQYGLVRSLNVHVAGSLFI 1127 6777777************************************************9987 PP >> sp|Q9VW14|MRM3_DROME rRNA methyltransferase 3, mitochondrial OS=Drosophila melanogaster OX=7227 GN=CG14100 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 74.3 0.0 1.2e-24 8.4e-21 2 141 .. 236 392 .. 235 393 .. 0.80 Alignments for each domain: == domain 1 score: 74.3 bits; conditional E-value: 1.2e-24 EEEEEES-SSHHHHHHHHHHHHHTT-SEEEEESGSSCTCGSHHHHHHCTTHHHHHHEEEESSHHH...HHHHHCCCTEEEEEEST...T CS SpoU_methylase 2 livvlveiedpgNlGaiaRtakafgaselvlvkprsadllnpeakraaagalevlelkvvdslee...alkelkeegfalvatsl...k 84 ++vv+++i++p+NlG+i+Rt++a+ +s++v+++ + d+++++a+r ++g ++ ++++ +++e ++ +++++ +++ k sp|Q9VW14|MRM3_DROME 236 ITVVCDNIREPNNLGSIIRTCAALPCSQVVVTH-GCCDPWESKALRGGCGGQFRVPIRDDVTWDElalTIPPEAADDCHVFIAETnqrK 323 78*****************************88.8889******************987777777211444444446888777664333 PP -SEGCCSHHH......TTEEEEEEBTTTBS-HHHHCTSS......EEEE----TTSSS--HHHHHHHHH CS SpoU_methylase 85 gkslseidlk......ekvalvfGnEdkGlsnellelad......rlvkIpmsgeveSLNvavAvaill 141 ++ ++id++ + +l++G E++G+s+e+ + ++ ++++Ip + ++SLNva A +ll sp|Q9VW14|MRM3_DROME 324 RENNQTIDYAdikglgAHNLLIIGGESHGVSEEAYRFLNlvggkgKCIYIPLAAGIDSLNVASALTLLL 392 33333666667779999****************98777632243469*****************99987 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (142 nodes) Target sequences: 13786 (7394964 residues searched) Passed MSV filter: 311 (0.0225591); expected 275.7 (0.02) Passed bias filter: 255 (0.018497); expected 275.7 (0.02) Passed Vit filter: 20 (0.00145075); expected 13.8 (0.001) Passed Fwd filter: 2 (0.000145075); expected 0.1 (1e-05) Initial search space (Z): 13786 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.08u 0.00s 00:00:00.08 Elapsed: 00:00:00.04 # Mc/sec: 26252.12 // [ok]