RPS-BLAST 2.2.26 [Sep-21-2011]

Database: CDD.v3.10 
           44,354 sequences; 10,937,602 total letters

Searching..................................................done

Query= psy5757
         (139 letters)



>gnl|CDD|216320 pfam01135, PCMT, Protein-L-isoaspartate(D-aspartate)
           O-methyltransferase (PCMT). 
          Length = 210

 Score = 58.5 bits (142), Expect = 2e-11
 Identities = 25/63 (39%), Positives = 35/63 (55%), Gaps = 10/63 (15%)

Query: 22  IIQPH----------LNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLVNSSI 71
           I  PH          L    +VLEIGSGSGYLT   + ++   G V+ IEH+P+LV  + 
Sbjct: 56  ISAPHMHAMMLELLELKPGMRVLEIGSGSGYLTACFARMVGEVGLVVSIEHIPELVEIAR 115

Query: 72  QNI 74
           +N+
Sbjct: 116 RNL 118


>gnl|CDD|232816 TIGR00080, pimt, protein-L-isoaspartate(D-aspartate)
           O-methyltransferase.  This is an all-kingdom (but not
           all species) full-length ortholog enzyme for repairing
           aging proteins. Among the prokaryotes, the gene name is
           pcm. Among eukaryotes, pimt [Protein fate, Protein
           modification and repair].
          Length = 215

 Score = 52.5 bits (126), Expect = 3e-09
 Identities = 21/63 (33%), Positives = 35/63 (55%), Gaps = 1/63 (1%)

Query: 13  ISAILTYLSIIQ-PHLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLVNSSI 71
           ISA      + +   L    KVLEIG+GSGY   +++E++   G V+ IE +P+L   + 
Sbjct: 60  ISAPHMVAMMTELLELKPGMKVLEIGTGSGYQAAVLAEIVGRDGLVVSIERIPELAEKAE 119

Query: 72  QNI 74
           + +
Sbjct: 120 RRL 122


>gnl|CDD|184409 PRK13942, PRK13942, protein-L-isoaspartate O-methyltransferase;
           Provisional.
          Length = 212

 Score = 51.9 bits (125), Expect = 5e-09
 Identities = 19/49 (38%), Positives = 33/49 (67%)

Query: 27  LNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLVNSSIQNIL 75
           L E  KVLEIG+GSGY   +++E++  +G+V+ IE +P+L   + + + 
Sbjct: 74  LKEGMKVLEIGTGSGYHAAVVAEIVGKSGKVVTIERIPELAEKAKKTLK 122


>gnl|CDD|225316 COG2518, Pcm, Protein-L-isoaspartate carboxylmethyltransferase
           [Posttranslational modification, protein turnover,
           chaperones].
          Length = 209

 Score = 49.6 bits (119), Expect = 4e-08
 Identities = 17/48 (35%), Positives = 29/48 (60%), Gaps = 3/48 (6%)

Query: 27  LNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLVNSSIQNI 74
           L    +VLEIG+GSGY   +++ L    G+V+ IE + +L   + +N+
Sbjct: 70  LKPGDRVLEIGTGSGYQAAVLARL---VGRVVSIERIEELAEQARRNL 114


>gnl|CDD|222415 pfam13847, Methyltransf_31, Methyltransferase domain.  This
          family appears to be have methyltransferase activity.
          Length = 151

 Score = 46.3 bits (110), Expect = 4e-07
 Identities = 15/47 (31%), Positives = 29/47 (61%)

Query: 27 LNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLVNSSIQN 73
          L    KVL++G G+GYLT +++E +    +V+GI+   + +  + +N
Sbjct: 1  LKSGIKVLDLGCGTGYLTFILAEKLGPGAEVVGIDISEEAIEKAKEN 47


>gnl|CDD|225317 COG2519, GCD14, tRNA(1-methyladenosine) methyltransferase and
           related methyltransferases [Translation, ribosomal
           structure and biogenesis].
          Length = 256

 Score = 39.2 bits (92), Expect = 2e-04
 Identities = 16/49 (32%), Positives = 23/49 (46%), Gaps = 5/49 (10%)

Query: 13  ISAILTYLSIIQPHLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIE 61
              I+  L I     +  S+VLE G+GSG LT  ++  +   G V   E
Sbjct: 83  AGYIVARLGI-----SPGSRVLEAGTGSGALTAYLARAVGPEGHVTTYE 126


>gnl|CDD|140001 PRK13944, PRK13944, protein-L-isoaspartate O-methyltransferase;
           Provisional.
          Length = 205

 Score = 38.3 bits (89), Expect = 5e-04
 Identities = 17/45 (37%), Positives = 27/45 (60%)

Query: 32  KVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLVNSSIQNILH 76
           K+LE+G+GSGY   + +E +   G+V  +E V +L   + QNI  
Sbjct: 75  KILEVGTGSGYQAAVCAEAIERRGKVYTVEIVKELAIYAAQNIER 119


>gnl|CDD|237568 PRK13943, PRK13943, protein-L-isoaspartate O-methyltransferase;
           Provisional.
          Length = 322

 Score = 37.9 bits (88), Expect = 7e-04
 Identities = 20/65 (30%), Positives = 34/65 (52%), Gaps = 5/65 (7%)

Query: 26  HLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLVNSSIQNILHSNARLLTDG 85
            L++  +VLEIG G+GY   ++S ++   G V+ +E+     +  I  I   N R L   
Sbjct: 77  GLDKGMRVLEIGGGTGYNAAVMSRVVGEKGLVVSVEY-----SRKICEIAKRNVRRLGIE 131

Query: 86  HIKFV 90
           ++ FV
Sbjct: 132 NVIFV 136


>gnl|CDD|178974 PRK00312, pcm, protein-L-isoaspartate O-methyltransferase;
           Reviewed.
          Length = 212

 Score = 37.5 bits (88), Expect = 8e-04
 Identities = 16/62 (25%), Positives = 29/62 (46%), Gaps = 4/62 (6%)

Query: 13  ISAILTYLSIIQP-HLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLVNSSI 71
           IS       + +   L    +VLEIG+GSGY   +++ L+    +V  +E +  L   + 
Sbjct: 61  ISQPYMVARMTELLELKPGDRVLEIGTGSGYQAAVLAHLV---RRVFSVERIKTLQWEAK 117

Query: 72  QN 73
           + 
Sbjct: 118 RR 119


>gnl|CDD|173391 PTZ00098, PTZ00098, phosphoethanolamine N-methyltransferase;
           Provisional.
          Length = 263

 Score = 37.3 bits (86), Expect = 0.001
 Identities = 31/87 (35%), Positives = 44/87 (50%), Gaps = 11/87 (12%)

Query: 11  GGISAILTYLSIIQPHLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLVNSS 70
           GGI A    LS I+  LNENSKVL+IGSG G     I+E   +   V G++   ++VN  
Sbjct: 36  GGIEATTKILSDIE--LNENSKVLDIGSGLGGGCKYINEKYGA--HVHGVDICEKMVNI- 90

Query: 71  IQNILHSNARLLTDGHIKFVALGMIKR 97
                 +  R      I+F A  ++K+
Sbjct: 91  ------AKLRNSDKNKIEFEANDILKK 111


>gnl|CDD|201699 pfam01269, Fibrillarin, Fibrillarin. 
          Length = 229

 Score = 36.8 bits (86), Expect = 0.002
 Identities = 20/71 (28%), Positives = 34/71 (47%), Gaps = 14/71 (19%)

Query: 5   KIGAAI-GGISAILTYLSIIQPHLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHV 63
           K+ AAI  G+  I          +   SKVL +G+ SG   + +S+++   G V  +E  
Sbjct: 57  KLAAAILKGLDNI---------PIKPGSKVLYLGAASGTTVSHVSDIVGPEGVVYAVEFS 107

Query: 64  P----QLVNSS 70
           P    +L+N +
Sbjct: 108 PRVGRELINMA 118


>gnl|CDD|221804 pfam12847, Methyltransf_18, Methyltransferase domain.  Protein in
           this family function as methyltransferases.
          Length = 104

 Score = 35.4 bits (82), Expect = 0.002
 Identities = 21/80 (26%), Positives = 39/80 (48%), Gaps = 9/80 (11%)

Query: 29  ENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLVNSSIQNILHSNARLLTDGHIK 88
             ++VL+IG G+G L   ++ L     +V G++  P+++      +   NA+L     I 
Sbjct: 1   PGARVLDIGCGTGSLAIELARLFP-GARVTGVDLSPEML-----ELARENAKLALGPRIT 54

Query: 89  FVALGMIKRIETVELMMKFD 108
           FV  G     + ++L+  FD
Sbjct: 55  FVQ-GDAP--DALDLLEGFD 71


>gnl|CDD|181382 PRK08317, PRK08317, hypothetical protein; Provisional.
          Length = 241

 Score = 36.5 bits (85), Expect = 0.002
 Identities = 11/61 (18%), Positives = 26/61 (42%), Gaps = 2/61 (3%)

Query: 26 HLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLVNSSIQNI--LHSNARLLT 83
           +    +VL++G G G     ++  +   G+V+GI+    ++  + +    L  N   + 
Sbjct: 16 AVQPGDRVLDVGCGPGNDARELARRVGPEGRVVGIDRSEAMLALAKERAAGLGPNVEFVR 75

Query: 84 D 84
           
Sbjct: 76 G 76


>gnl|CDD|240291 PTZ00146, PTZ00146, fibrillarin; Provisional.
          Length = 293

 Score = 33.9 bits (78), Expect = 0.014
 Identities = 19/58 (32%), Positives = 30/58 (51%), Gaps = 10/58 (17%)

Query: 5   KIGAAI-GGISAILTYLSIIQPHLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIE 61
           K+ AAI GG++ I          +   SKVL +G+ SG   + +S+L+   G V  +E
Sbjct: 116 KLAAAIIGGVANI---------PIKPGSKVLYLGAASGTTVSHVSDLVGPEGVVYAVE 164


>gnl|CDD|234492 TIGR04188, methyltr_grsp, methyltransferase, ATP-grasp peptide
           maturase system.  Members of this protein family are
           predicted SAM-dependent methyltransferases that
           regularly occur in the context of a putative peptide
           modification ATP-grasp enzyme (TIGR04187, related to
           enzymes of microviridin maturation) and a putative
           ribosomal peptide modification target (TIGR04186).
          Length = 363

 Score = 33.5 bits (77), Expect = 0.026
 Identities = 12/36 (33%), Positives = 19/36 (52%), Gaps = 1/36 (2%)

Query: 32  KVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLV 67
           +VLEIG+G+GY   ++   +     V  +E  P L 
Sbjct: 110 RVLEIGTGTGYSAALLCHRLGD-DNVTSVEVDPGLA 144


>gnl|CDD|223574 COG0500, SmtA, SAM-dependent methyltransferases [Secondary
           metabolites biosynthesis, transport, and catabolism /
           General function prediction only].
          Length = 257

 Score = 32.6 bits (71), Expect = 0.042
 Identities = 19/84 (22%), Positives = 33/84 (39%), Gaps = 6/84 (7%)

Query: 25  PHLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLVNSSIQNILHSNARLLTD 84
             L     VL+IG G+G L  +   L      V+G++  P+++      +  + A     
Sbjct: 44  RLLPGGLGVLDIGCGTGRLALLA-RLGGRGAYVVGVDLSPEML-----ALARARAEGAGL 97

Query: 85  GHIKFVALGMIKRIETVELMMKFD 108
           G + FV    +  +   E    FD
Sbjct: 98  GLVDFVVADALGGVLPFEDSASFD 121


>gnl|CDD|233708 TIGR02072, BioC, biotin biosynthesis protein BioC.  This enzyme,
          which is found in biotin biosynthetic gene clusters in
          proteobacteria, firmicutes, green-sulfur bacteria,
          fusobacterium and bacteroides, is believed to carry out
          an enzymatic step prior to the formation of
          pimeloyl-CoA (although attribution of this annotation
          is not traceable). The enzyme appears related to
          methyltransferases by homology [Biosynthesis of
          cofactors, prosthetic groups, and carriers, Biotin].
          Length = 240

 Score = 32.3 bits (74), Expect = 0.046
 Identities = 9/21 (42%), Positives = 14/21 (66%)

Query: 29 ENSKVLEIGSGSGYLTNMISE 49
            + VL+IG G+GYLT  + +
Sbjct: 34 IPASVLDIGCGTGYLTRALLK 54


>gnl|CDD|177970 PLN02336, PLN02336, phosphoethanolamine N-methyltransferase.
          Length = 475

 Score = 32.4 bits (74), Expect = 0.049
 Identities = 22/65 (33%), Positives = 35/65 (53%), Gaps = 8/65 (12%)

Query: 20 LSIIQPHLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLV--NSSIQNILHS 77
          LS++ P+  E   VLE+G+G G  T    EL    GQVI ++ +  ++  N SI N  + 
Sbjct: 30 LSLLPPY--EGKSVLELGAGIGRFT---GELAKKAGQVIALDFIESVIKKNESI-NGHYK 83

Query: 78 NARLL 82
          N + +
Sbjct: 84 NVKFM 88


>gnl|CDD|222171 pfam13489, Methyltransf_23, Methyltransferase domain.  This
          family appears to be a methyltransferase domain.
          Length = 154

 Score = 32.1 bits (73), Expect = 0.049
 Identities = 16/67 (23%), Positives = 28/67 (41%), Gaps = 7/67 (10%)

Query: 25 PHLNENSKVLEIGSGSGYLTNMISELMNSTG-QVIGIEHVPQLVNSSIQNILHSNARLLT 83
          P L    +VL+IG G+G L  ++ E     G  V G++  P      I ++  +    + 
Sbjct: 18 PRLKPGGRVLDIGCGTGILLRLLRER----GFDVTGVD--PSPAAVLIFSLFDAPDPAVL 71

Query: 84 DGHIKFV 90
           G    +
Sbjct: 72 AGKYDLI 78


>gnl|CDD|222236 pfam13578, Methyltransf_24, Methyltransferase domain.  This
          family appears to be a methyltransferase domain.
          Length = 105

 Score = 31.5 bits (72), Expect = 0.053
 Identities = 14/53 (26%), Positives = 24/53 (45%), Gaps = 4/53 (7%)

Query: 34 LEIGSGSGYLTNMISELM--NSTGQVIGIEHVPQLVNSSI--QNILHSNARLL 82
          +EIG  SG  T  +++ +  N  G++  I+  P     +   +  L    RLL
Sbjct: 1  VEIGVYSGVSTLWLAQALKDNGLGRLYSIDPWPGAEAGANLRKAGLADRVRLL 53


>gnl|CDD|188051 TIGR00446, nop2p, NOL1/NOP2/sun family putative RNA methylase.
           [Protein synthesis, tRNA and rRNA base modification].
          Length = 264

 Score = 32.4 bits (74), Expect = 0.055
 Identities = 18/61 (29%), Positives = 30/61 (49%), Gaps = 8/61 (13%)

Query: 21  SIIQP---HLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLVNSSIQNILHS 77
           S+I P      E+ +VL++ +  G  T  IS+LM + G ++  E     ++ S    L S
Sbjct: 60  SMIPPIALEPREDERVLDMAAAPGGKTTQISQLMKNKGCIVANE-----ISKSRTKALIS 114

Query: 78  N 78
           N
Sbjct: 115 N 115


>gnl|CDD|128898 smart00650, rADc, Ribosomal RNA adenine dimethylases. 
          Length = 169

 Score = 31.7 bits (73), Expect = 0.062
 Identities = 15/36 (41%), Positives = 18/36 (50%), Gaps = 3/36 (8%)

Query: 26 HLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIE 61
          +L     VLEIG G G LT    EL+    +V  IE
Sbjct: 10 NLRPGDTVLEIGPGKGALTE---ELLERAKRVTAIE 42


>gnl|CDD|234563 TIGR04364, methyltran_FxLD, methyltransferase, FxLD system.
           Members of this family resemble occur regularly in the
           vicinity of lantibiotic biosynthesis enzymes and their
           probable target, the FxLD family of putative ribosomal
           natural product precursor (TIGR04363). Members resemble
           protein-L-isoaspartate O-methyltransferase (TIGR00080)
           and a predicted methyltranserase, TIGR04188, of another
           putative peptide modification system.
          Length = 394

 Score = 32.0 bits (73), Expect = 0.088
 Identities = 21/61 (34%), Positives = 34/61 (55%), Gaps = 12/61 (19%)

Query: 7   GAAIGGISAILTYLSIIQPHLNENS------KVLEIGSGSGYLTNMISELMNSTGQVIGI 60
           GAA+  +SA       IQ  + E +      +VLEIGSG GY   +++EL+  +G+V  +
Sbjct: 59  GAALSSVSA-----PHIQAMMLEQAGVEPGMRVLEIGSG-GYNAALLAELVGPSGEVTTV 112

Query: 61  E 61
           +
Sbjct: 113 D 113


>gnl|CDD|215900 pfam00398, RrnaAD, Ribosomal RNA adenine dimethylase. 
          Length = 254

 Score = 31.5 bits (72), Expect = 0.099
 Identities = 21/65 (32%), Positives = 30/65 (46%), Gaps = 10/65 (15%)

Query: 26 HLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLVNSSIQNILHSNARLLTDG 85
          +L E+  VLEIG G G LT   +EL     QV+ IE  P+L             +L    
Sbjct: 26 NLQESDTVLEIGPGKGALT---TELAKRAKQVVAIEIDPRLAK-------RLQEKLALHP 75

Query: 86 HIKFV 90
          +++ V
Sbjct: 76 NVEVV 80


>gnl|CDD|233115 TIGR00755, ksgA, dimethyladenosine transferase.  In both E. coli
          and Saccharomyces cerevisiae, this protein is
          responsible for the dimethylation of two adjacent
          adenosine residues in a conserved hairpin of 16S rRNA
          in bacteria, 18S rRNA in eukaryotes. This adjacent
          dimethylation is the only rRNA modification shared by
          bacteria and eukaryotes. A single member of this family
          is present in each of the first 20 completed microbial
          genomes. This protein is essential in yeast, but not in
          E. coli, where its deletion leads to resistance to the
          antibiotic kasugamycin. Alternate name:
          S-adenosylmethionine--6-N',N'-adenosyl (rRNA)
          dimethyltransferase [Protein synthesis, tRNA and rRNA
          base modification].
          Length = 253

 Score = 31.1 bits (71), Expect = 0.15
 Identities = 15/49 (30%), Positives = 21/49 (42%), Gaps = 8/49 (16%)

Query: 13 ISAILTYLSIIQPHLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIE 61
          I  I+   ++      E   VLEIG G G LT  + +      +V  IE
Sbjct: 18 IQKIVEAANV-----LEGDVVLEIGPGLGALTEPLLKR---AKKVTAIE 58


>gnl|CDD|234740 PRK00377, cbiT, cobalt-precorrin-6Y C(15)-methyltransferase;
          Provisional.
          Length = 198

 Score = 30.9 bits (70), Expect = 0.15
 Identities = 14/47 (29%), Positives = 27/47 (57%)

Query: 27 LNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLVNSSIQN 73
          L +   +L+IG G+G +T   S L+  TG+V  ++   + +N + +N
Sbjct: 38 LRKGDMILDIGCGTGSVTVEASLLVGETGKVYAVDKDEKAINLTRRN 84


>gnl|CDD|237852 PRK14896, ksgA, 16S ribosomal RNA methyltransferase KsgA/Dim1
          family protein; Provisional.
          Length = 258

 Score = 31.0 bits (71), Expect = 0.15
 Identities = 14/35 (40%), Positives = 17/35 (48%), Gaps = 3/35 (8%)

Query: 27 LNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIE 61
            +   VLEIG G G LT+   EL     +V  IE
Sbjct: 27 DTDGDPVLEIGPGKGALTD---ELAKRAKKVYAIE 58


>gnl|CDD|237898 PRK15068, PRK15068, tRNA mo(5)U34 methyltransferase; Provisional.
          Length = 322

 Score = 30.6 bits (70), Expect = 0.20
 Identities = 25/96 (26%), Positives = 44/96 (45%), Gaps = 26/96 (27%)

Query: 22  IIQPHLN--ENSKVLEIGSGSGY-LTNMISELMNSTG--QVIGIEHVP-QLVN---SSIQ 72
            + PHL+  +   VL++G G+GY +  M+       G   V+GI+  P QL      +++
Sbjct: 113 RVLPHLSPLKGRTVLDVGCGNGYHMWRMLGA-----GAKLVVGID--PSQLFLCQFEAVR 165

Query: 73  NILHSNARLLTDGHIKFVALGMIKRIETVELMMKFD 108
            +L ++ R         + LG    IE +  +  FD
Sbjct: 166 KLLGNDQR------AHLLPLG----IEQLPALKAFD 191


>gnl|CDD|223109 COG0030, KsgA, Dimethyladenosine transferase (rRNA methylation)
          [Translation, ribosomal structure and biogenesis].
          Length = 259

 Score = 30.3 bits (69), Expect = 0.25
 Identities = 13/36 (36%), Positives = 18/36 (50%), Gaps = 3/36 (8%)

Query: 26 HLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIE 61
          +++    VLEIG G G LT     L+    +V  IE
Sbjct: 27 NISPGDNVLEIGPGLGALT---EPLLERAARVTAIE 59


>gnl|CDD|222295 pfam13659, Methyltransf_26, Methyltransferase domain.  This
          family contains methyltransferase domains.
          Length = 117

 Score = 29.7 bits (67), Expect = 0.27
 Identities = 9/30 (30%), Positives = 16/30 (53%), Gaps = 1/30 (3%)

Query: 32 KVLEIGSGSGYLTNMISELMNSTGQVIGIE 61
          +VL+ G+GSG      +       +V+G+E
Sbjct: 3  RVLDPGAGSGAFLLAAARAGPD-ARVVGVE 31


>gnl|CDD|222313 pfam13679, Methyltransf_32, Methyltransferase domain.  This family
           appears to be a methyltransferase domain.
          Length = 143

 Score = 29.9 bits (68), Expect = 0.29
 Identities = 17/94 (18%), Positives = 35/94 (37%), Gaps = 13/94 (13%)

Query: 19  YLSIIQPHLNENSK------VLEIGSGSGYLTN-MISELMNST--GQVIGIEHVPQLVNS 69
               ++P + E         V++ G+G GYL   +   L       +V+GI+   +LV  
Sbjct: 9   LAEFVEPLIKEVLDKTGCITVVDHGAGKGYLGFILYYLLFKRKYGVRVVGIDARAELVEK 68

Query: 70  SIQNILHSNARLLTDGHIKFVALGMIKRIETVEL 103
           +         +L  +  + F+   +     +  L
Sbjct: 69  A-NAR---AQKLGFNKRMSFLEGTIAGATVSARL 98


>gnl|CDD|222287 pfam13649, Methyltransf_25, Methyltransferase domain.  This
          family appears to be a methyltransferase domain.
          Length = 97

 Score = 28.8 bits (65), Expect = 0.36
 Identities = 7/29 (24%), Positives = 16/29 (55%), Gaps = 2/29 (6%)

Query: 33 VLEIGSGSGYLTNMISELMNSTGQVIGIE 61
          +L++G G+G +   ++    S   V G++
Sbjct: 1  ILDLGCGTGRVLRALARAGPS--SVTGVD 27


>gnl|CDD|217012 pfam02390, Methyltransf_4, Putative methyltransferase.  This is a
          family of putative methyltransferases. The aligned
          region contains the GXGXG S-AdoMet binding site
          suggesting a putative methyltransferase activity.
          Length = 198

 Score = 29.6 bits (67), Expect = 0.39
 Identities = 22/82 (26%), Positives = 31/82 (37%), Gaps = 15/82 (18%)

Query: 28 NENSKVLEIGSGSGYLTNMISELMNSTGQVIGIE-HVP---QLVNSSIQ-------NILH 76
          NE    LEIG G G     +++  N     IGIE  VP   + +   I         IL 
Sbjct: 19 NEQPLFLEIGCGMGDFLVAMAK-KNPDKLFIGIEIRVPGVDKALKKIIALRGLQNLRILC 77

Query: 77 SNARLLTDGHIKFVALGMIKRI 98
           +A  L          G +++I
Sbjct: 78 GDAMKLLP---NLFPDGSLQKI 96


>gnl|CDD|234649 PRK00121, trmB, tRNA (guanine-N(7)-)-methyltransferase; Reviewed.
          Length = 202

 Score = 29.4 bits (67), Expect = 0.48
 Identities = 21/69 (30%), Positives = 29/69 (42%), Gaps = 14/69 (20%)

Query: 28  NENSKVLEIGSGSG-YLTNMISELMNSTGQVIGIE-HVP-------QLVNSSIQN--ILH 76
           N+    LEIG G G +L  M     N     IGIE H P       ++    + N  +L 
Sbjct: 39  NDAPIHLEIGFGKGEFLVEMA--KANPDINFIGIEVHEPGVGKALKKIEEEGLTNLRLLC 96

Query: 77  SNA-RLLTD 84
            +A  +L D
Sbjct: 97  GDAVEVLLD 105


>gnl|CDD|100107 cd02440, AdoMet_MTases, S-adenosylmethionine-dependent
          methyltransferases (SAM or AdoMet-MTase), class I;
          AdoMet-MTases are enzymes that use
          S-adenosyl-L-methionine (SAM or AdoMet) as a substrate
          for methyltransfer, creating the product
          S-adenosyl-L-homocysteine (AdoHcy). There are at least
          five structurally distinct families of AdoMet-MTases,
          class I being the largest and most diverse. Within this
          class enzymes can be classified by different substrate
          specificities (small molecules, lipids, nucleic acids,
          etc.) and different target atoms for methylation
          (nitrogen, oxygen, carbon, sulfur, etc.).
          Length = 107

 Score = 28.5 bits (64), Expect = 0.48
 Identities = 8/30 (26%), Positives = 17/30 (56%), Gaps = 2/30 (6%)

Query: 32 KVLEIGSGSGYLTNMISELMNSTGQVIGIE 61
          +VL++G G+G L   ++       +V G++
Sbjct: 1  RVLDLGCGTGALALALASGPG--ARVTGVD 28


>gnl|CDD|223298 COG0220, COG0220, Predicted S-adenosylmethionine-dependent
           methyltransferase [General function prediction only].
          Length = 227

 Score = 29.2 bits (66), Expect = 0.58
 Identities = 20/68 (29%), Positives = 29/68 (42%), Gaps = 13/68 (19%)

Query: 28  NENSKVLEIGSGSG-YLTNMISELMNSTGQVIGIE-HVP-------QLVNSSIQN--ILH 76
           N    VLEIG G G +L  M     N     +GIE  VP       ++    ++N  +L 
Sbjct: 47  NNAPIVLEIGFGMGEFLVEM--AKKNPEKNFLGIEIRVPGVAKALKKIKELGLKNLRLLC 104

Query: 77  SNARLLTD 84
            +A  + D
Sbjct: 105 GDAVEVLD 112


>gnl|CDD|225136 COG2226, UbiE, Methylase involved in ubiquinone/menaquinone
           biosynthesis [Coenzyme metabolism].
          Length = 238

 Score = 29.2 bits (66), Expect = 0.65
 Identities = 17/70 (24%), Positives = 30/70 (42%), Gaps = 15/70 (21%)

Query: 22  IIQPHLNENSKVLEIGSGSGYLTNMISELMNS--TGQVIGIEHVPQ--------LVNSSI 71
           I    +    KVL++  G+G   +M   L  S  TG+V+G++            L    +
Sbjct: 44  ISLLGIKPGDKVLDVACGTG---DMALLLAKSVGTGEVVGLDISESMLEVAREKLKKKGV 100

Query: 72  QNI--LHSNA 79
           QN+  +  +A
Sbjct: 101 QNVEFVVGDA 110


>gnl|CDD|226472 COG3963, COG3963, Phospholipid N-methyltransferase [Lipid
          metabolism].
          Length = 194

 Score = 28.9 bits (65), Expect = 0.71
 Identities = 20/65 (30%), Positives = 28/65 (43%), Gaps = 7/65 (10%)

Query: 9  AIGGI---SAILT--YLSIIQPHLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHV 63
           +G I   S+IL     S+I P       VLE+G G+G +T  I         +  IE+ 
Sbjct: 25 TVGAILPSSSILARKMASVIDPESGLP--VLELGPGTGVITKAILSRGVRPESLTAIEYS 82

Query: 64 PQLVN 68
          P  V 
Sbjct: 83 PDFVC 87


>gnl|CDD|219760 pfam08242, Methyltransf_12, Methyltransferase domain.  Members of
          this family are SAM dependent methyltransferases.
          Length = 98

 Score = 28.1 bits (63), Expect = 0.74
 Identities = 8/28 (28%), Positives = 14/28 (50%), Gaps = 1/28 (3%)

Query: 34 LEIGSGSGYLTNMISELMNSTGQVIGIE 61
          L+IG G+G L   + E      +  G++
Sbjct: 1  LDIGCGTGTLLRALLEA-LPGLEYTGVD 27


>gnl|CDD|225151 COG2242, CobL, Precorrin-6B methylase 2 [Coenzyme metabolism].
          Length = 187

 Score = 28.7 bits (65), Expect = 0.75
 Identities = 12/47 (25%), Positives = 23/47 (48%), Gaps = 1/47 (2%)

Query: 27 LNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLVNSSIQN 73
               ++ +IG+G+G +T   +    S G+VI IE   + +    +N
Sbjct: 32 PRPGDRLWDIGAGTGSITIEWALAGPS-GRVIAIERDEEALELIERN 77


>gnl|CDD|237872 PRK14968, PRK14968, putative methyltransferase; Provisional.
          Length = 188

 Score = 28.7 bits (65), Expect = 0.84
 Identities = 7/15 (46%), Positives = 12/15 (80%)

Query: 29 ENSKVLEIGSGSGYL 43
          +  +VLE+G+GSG +
Sbjct: 23 KGDRVLEVGTGSGIV 37


>gnl|CDD|233880 TIGR02469, CbiT, precorrin-6Y C5,15-methyltransferase
          (decarboxylating), CbiT subunit.  This model recognizes
          the CbiT methylase which is responsible, in part (along
          with CbiE), for methylating precorrin-6y (or
          cobalt-precorrin-6y) at both the 5 and 15 positions as
          well as the concomitant decarbozylation at C-12. In
          many organisms, this protein is fused to the CbiE
          subunit. The fused protein, when found in organisms
          catalyzing the oxidative version of the cobalamin
          biosynthesis pathway, is called CobL [Biosynthesis of
          cofactors, prosthetic groups, and carriers, Heme,
          porphyrin, and cobalamin].
          Length = 124

 Score = 28.1 bits (63), Expect = 0.87
 Identities = 16/66 (24%), Positives = 35/66 (53%), Gaps = 6/66 (9%)

Query: 20 LSIIQPHLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQ---LVNSSIQNILH 76
          LS ++    +   + +IG+G+G +T   + L+ + G+V  IE  P+   L+  +++    
Sbjct: 12 LSKLRLRPGD--VLWDIGAGTGSVTIEAARLVPN-GRVYAIERNPEALDLIERNLRRFGV 68

Query: 77 SNARLL 82
          SN  ++
Sbjct: 69 SNIVIV 74


>gnl|CDD|241257 cd01223, PH_Vav, Vav pleckstrin homology (PH) domain.  Vav acts as
           a guanosine nucleotide exchange factor (GEF) for Rho/Rac
           proteins. They control processes including T cell
           activation, phagocytosis, and migration of cells. The
           Vav subgroup of Dbl GEFs consists of three family
           members (Vav1, Vav2, and Vav3) in mammals. Vav1 is
           preferentially expressed in the hematopoietic system,
           while Vav2 and Vav3 are described by broader expression
           patterns. Mammalian Vav proteins consist of a calponin
           homology (CH) domain, an acidic region, a catalytic Dbl
           homology (DH) domain, a PH domain, a zinc finger
           cysteine rich domain (C1/CRD), and an SH2 domain,
           flanked by two SH3 domains. In invertebrates such as
           Drosophila and C. elegans, Vav is missing the N-terminal
           SH3 domain. The DH domain is involved in RhoGTPase
           recognition and selectivity and stimulates the
           reorganization of the switch regions for GDP/GTP
           exchange. The PH domain is implicated in directing
           membrane localization, allosteric regulation of guanine
           nucleotide exchange activity, and as a phospholipid-
           dependent regulator of GEF activity. Vavs bind
           RhoGTPases including Rac1, RhoA, RhoG, and Cdc42, while
           other members of the GEF family are specific for a
           single RhoGTPase. This promiscuity is thought to be a
           result of its CRD. PH domains have diverse functions,
           but in general are involved in targeting proteins to the
           appropriate cellular location or in the interaction with
           a binding partner. They share little sequence
           conservation, but all have a common fold, which is
           electrostatically polarized. Less than 10% of PH domains
           bind phosphoinositide phosphates (PIPs) with high
           affinity and specificity. PH domains are distinguished
           from other PIP-binding domains by their specific
           high-affinity binding to PIPs with two vicinal phosphate
           groups: PtdIns(3,4)P2, PtdIns(4,5)P2 or PtdIns(3,4,5)P3
           which results in targeting some PH domain proteins to
           the plasma membrane. A few display strong specificity in
           lipid binding. Any specificity is usually determined by
           loop regions or insertions in the N-terminus of the
           domain, which are not conserved across all PH domains.
           PH domains are found in cellular signaling proteins such
           as serine/threonine kinase, tyrosine kinases, regulators
           of G-proteins, endocytotic GTPases, adaptors, as well as
           cytoskeletal associated molecules and in lipid
           associated enzymes.polarized. PH domains also have
           diverse functions. They are often involved in targeting
           proteins to the plasma membrane, but only a few (less
           than 10%) display strong specificity in binding inositol
           phosphates. Any specificity is usually determined by
           loop regions or insertions in the N-terminus of the
           domain, which are not conserved across all PH domains.
           PH domains are found in cellular signaling proteins such
           as serine/threonine kinases, tyrosine kinases,
           regulators of G-proteins, endocytotic GTPases, adaptors,
           cytoskeletal associated molecules, and in lipid
           associated enzymes.
          Length = 127

 Score = 28.0 bits (63), Expect = 0.91
 Identities = 15/49 (30%), Positives = 21/49 (42%), Gaps = 9/49 (18%)

Query: 65  QLVNSSIQNILHSNARLLTDGHIKFVALGMIKRIETVELMMKFDRYDFL 113
           Q +  ++   L    RL  DG +K         I++ E   K DRY FL
Sbjct: 3   QDLIENLNVSLADYGRLKIDGELK---------IKSHEDQKKKDRYAFL 42


>gnl|CDD|234708 PRK00274, ksgA, 16S ribosomal RNA methyltransferase KsgA/Dim1
          family protein; Reviewed.
          Length = 272

 Score = 28.5 bits (65), Expect = 0.97
 Identities = 12/36 (33%), Positives = 15/36 (41%), Gaps = 3/36 (8%)

Query: 26 HLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIE 61
                 VLEIG G G LT     L+    +V  +E
Sbjct: 39 GPQPGDNVLEIGPGLGALTE---PLLERAAKVTAVE 71


>gnl|CDD|219759 pfam08241, Methyltransf_11, Methyltransferase domain.  Members of
          this family are SAM dependent methyltransferases.
          Length = 92

 Score = 27.2 bits (61), Expect = 1.2
 Identities = 7/28 (25%), Positives = 15/28 (53%), Gaps = 2/28 (7%)

Query: 34 LEIGSGSGYLTNMISELMNSTGQVIGIE 61
          L++G G+G L   ++       +V G++
Sbjct: 1  LDVGCGTGLLAEALARRGG--ARVTGVD 26


>gnl|CDD|225100 COG2189, COG2189, Adenine specific DNA methylase Mod [DNA
           replication, recombination, and repair].
          Length = 590

 Score = 28.6 bits (64), Expect = 1.3
 Identities = 18/91 (19%), Positives = 42/91 (46%), Gaps = 11/91 (12%)

Query: 28  NENSKVLEIGSGSGYLTNMISELMNSTGQVIGIE---HVPQLVNSSIQNIL----HSNAR 80
           N+N  VL+  +GSG  T +  +L     + IGIE   ++  ++   ++ ++       ++
Sbjct: 452 NKNDLVLDFFAGSGTTTAVAHKL---NRKYIGIEQMEYIENIILERLKKVIDGEQGGISK 508

Query: 81  LLT-DGHIKFVALGMIKRIETVELMMKFDRY 110
                G + FV   +++  E +  ++   +Y
Sbjct: 509 DFDWKGGLSFVYAELLEENEELIQLILLSKY 539


>gnl|CDD|177877 PLN02233, PLN02233, ubiquinone biosynthesis methyltransferase.
          Length = 261

 Score = 28.3 bits (63), Expect = 1.3
 Identities = 15/48 (31%), Positives = 29/48 (60%), Gaps = 1/48 (2%)

Query: 33  VLEIGSGSGYLTNMISELMNSTGQVIGIEH-VPQLVNSSIQNILHSNA 79
           VL++  GSG L  ++SE + S G+V+G++    QL  ++ +  L + +
Sbjct: 77  VLDLCCGSGDLAFLLSEKVGSDGKVMGLDFSSEQLAVAASRQELKAKS 124


>gnl|CDD|219269 pfam07021, MetW, Methionine biosynthesis protein MetW.  This
          family consists of several bacterial and one archaeal
          methionine biosynthesis MetW proteins. Biosynthesis of
          methionine from homoserine in Pseudomonas putida takes
          place in three steps. The first step is the acylation
          of homoserine to yield an acyl-L-homoserine. This
          reaction is catalyzed by the products of the metXW
          genes and is equivalent to the first step in
          enterobacteria, gram-positive bacteria and fungi,
          except that in these microorganisms the reaction is
          catalyzed by a single polypeptide (the product of the
          metA gene in Escherichia coli and the met5 gene product
          in Neurospora crassa). In Pseudomonas putida, as in
          gram-positive bacteria and certain fungi, the second
          and third steps are a direct sulfhydrylation that
          converts the O-acyl-L-homoserine into homocysteine and
          further methylation to yield methionine. The latter
          reaction can be mediated by either of the two
          methionine synthetases present in the cells.
          Length = 193

 Score = 27.8 bits (62), Expect = 1.5
 Identities = 15/57 (26%), Positives = 23/57 (40%), Gaps = 2/57 (3%)

Query: 20 LSIIQPHLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLVNSSIQNILH 76
             I   +   S+VL++G G G L  ++ E     G+  GIE     V   +   L 
Sbjct: 4  FRYIARLIPPGSRVLDLGCGDGSLLYLLQEEKQVDGR--GIELDAAGVAECVAKGLS 58


>gnl|CDD|226607 COG4122, COG4122, Predicted O-methyltransferase [General function
          prediction only].
          Length = 219

 Score = 27.7 bits (62), Expect = 2.0
 Identities = 9/31 (29%), Positives = 17/31 (54%)

Query: 32 KVLEIGSGSGYLTNMISELMNSTGQVIGIEH 62
          ++LEIG+  GY    ++  +   G++  IE 
Sbjct: 62 RILEIGTAIGYSALWMALALPDDGRLTTIER 92


>gnl|CDD|240367 PTZ00338, PTZ00338, dimethyladenosine transferase-like protein;
          Provisional.
          Length = 294

 Score = 27.7 bits (62), Expect = 2.2
 Identities = 16/36 (44%), Positives = 22/36 (61%), Gaps = 3/36 (8%)

Query: 33 VLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLVN 68
          VLEIG G+G LT    +L+    +VI IE  P++V 
Sbjct: 40 VLEIGPGTGNLT---EKLLQLAKKVIAIEIDPRMVA 72


>gnl|CDD|225173 COG2264, PrmA, Ribosomal protein L11 methylase [Translation,
           ribosomal structure and biogenesis].
          Length = 300

 Score = 27.6 bits (62), Expect = 2.3
 Identities = 17/60 (28%), Positives = 31/60 (51%), Gaps = 4/60 (6%)

Query: 18  TYLSI--IQPHLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLVNSSIQNIL 75
           T L +  ++  L +   VL++G GSG L   I+       +V+G++  PQ V ++ +N  
Sbjct: 149 TSLCLEALEKLLKKGKTVLDVGCGSGIL--AIAAAKLGAKKVVGVDIDPQAVEAARENAR 206


>gnl|CDD|237007 PRK11873, arsM, arsenite S-adenosylmethyltransferase; Reviewed.
          Length = 272

 Score = 27.2 bits (61), Expect = 2.5
 Identities = 13/41 (31%), Positives = 25/41 (60%)

Query: 27  LNENSKVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLV 67
           L     VL++GSG G+   + +  +  TG+VIG++  P+++
Sbjct: 75  LKPGETVLDLGSGGGFDCFLAARRVGPTGKVIGVDMTPEML 115


>gnl|CDD|144813 pfam01358, PARP_regulatory, Poly A polymerase regulatory subunit.
          
          Length = 295

 Score = 26.7 bits (59), Expect = 4.6
 Identities = 12/44 (27%), Positives = 21/44 (47%), Gaps = 2/44 (4%)

Query: 17 LTYLSIIQPH--LNENSKVLEIGSGSGYLTNMISELMNSTGQVI 58
          L +L+ +  H  L+    V+ IGS  GY    + +  +  G +I
Sbjct: 44 LFFLNKLIRHNILDGIITVVYIGSAPGYHIRYLRDHFSDLGVII 87


>gnl|CDD|149213 pfam08003, Methyltransf_9, Protein of unknown function (DUF1698).
           This family contains many hypothetical proteins. It also
           includes two putative methyltransferase proteins.
          Length = 315

 Score = 26.6 bits (59), Expect = 5.1
 Identities = 13/40 (32%), Positives = 23/40 (57%), Gaps = 6/40 (15%)

Query: 25  PHLN--ENSKVLEIGSGSGY-LTNMISELMNSTGQVIGIE 61
           PHL+  +   +L++G G+GY +  M+ E       V+GI+
Sbjct: 109 PHLSPLKGRTILDVGCGNGYHMWRMLGE---GAALVVGID 145


>gnl|CDD|234637 PRK00107, gidB, 16S rRNA methyltransferase GidB; Reviewed.
          Length = 187

 Score = 25.9 bits (58), Expect = 7.8
 Identities = 9/19 (47%), Positives = 14/19 (73%)

Query: 23 IQPHLNENSKVLEIGSGSG 41
          I P+L    +VL++GSG+G
Sbjct: 39 IAPYLPGGERVLDVGSGAG 57


>gnl|CDD|161703 TIGR00091, TIGR00091, tRNA (guanine-N(7)-)-methyltransferase.
          This predicted S-adenosylmethionine-dependent
          methyltransferase is found in a single copy in most
          Bacteria. It is also found, with a short amino-terminal
          extension in eukaryotes. Its function is unknown. In E.
          coli, this protein flanks the DNA repair protein MutY,
          also called micA [Protein synthesis, tRNA and rRNA base
          modification].
          Length = 194

 Score = 25.8 bits (57), Expect = 8.2
 Identities = 21/82 (25%), Positives = 35/82 (42%), Gaps = 16/82 (19%)

Query: 28 NENSKVLEIGSGSG-YLTNMISELMNSTGQVIGIEHVPQLVNSSIQ----------NILH 76
          N+    LEIG G G +L +M  +  N     +GIE    +V ++            ++L 
Sbjct: 15 NKAPLHLEIGCGKGRFLIDMAKQ--NPDKNFLGIEIHTPIVLAANNKANKLGLKNLHVLC 72

Query: 77 SNARLLTDGHIKFVALGMIKRI 98
           +A  L D   KF   G + ++
Sbjct: 73 GDANELLD---KFFPDGSLSKV 91


>gnl|CDD|168661 PRK06753, PRK06753, hypothetical protein; Provisional.
          Length = 373

 Score = 25.8 bits (57), Expect = 8.6
 Identities = 28/91 (30%), Positives = 41/91 (45%), Gaps = 5/91 (5%)

Query: 1  MNLVKIGAAIGGISAILTYLSI-IQPHLNE-NSKVLEIGSGSGYLTNMISELMNSTGQVI 58
          M +  IGA IGG++A         +  + E N  V E+G+G G   N+I +L N      
Sbjct: 1  MKIAIIGAGIGGLTAAALLQEQGHEVKVFEKNESVKEVGAGIGIGDNVIKKLGNH-DLAK 59

Query: 59 GIEHVPQLVNSSIQNILHSNARLLTDGHIKF 89
          GI++  Q    S  N+L     LL    +K 
Sbjct: 60 GIKNAGQ--ILSTMNLLDDKGTLLNKVKLKS 88


>gnl|CDD|110227 pfam01209, Ubie_methyltran, ubiE/COQ5 methyltransferase family. 
          Length = 233

 Score = 25.5 bits (56), Expect = 9.4
 Identities = 10/33 (30%), Positives = 21/33 (63%)

Query: 29 ENSKVLEIGSGSGYLTNMISELMNSTGQVIGIE 61
            +K L++  G+G  T  +S+   S+G+V+G++
Sbjct: 47 RGNKFLDVAGGTGDWTFGLSDSAGSSGKVVGLD 79


>gnl|CDD|225139 COG2230, Cfa, Cyclopropane fatty acid synthase and related
           methyltransferases [Cell envelope biogenesis, outer
           membrane].
          Length = 283

 Score = 25.7 bits (57), Expect = 9.9
 Identities = 11/35 (31%), Positives = 16/35 (45%), Gaps = 2/35 (5%)

Query: 26  HLNENSKVLEIGSGSGYLTNMISELMNSTGQVIGI 60
            L     +L+IG G G L    +E    T  V+G+
Sbjct: 69  GLKPGMTLLDIGCGWGGLAIYAAEEYGVT--VVGV 101


>gnl|CDD|218990 pfam06325, PrmA, Ribosomal protein L11 methyltransferase (PrmA).
           This family consists of several Ribosomal protein L11
           methyltransferase (EC:2.1.1.-) sequences.
          Length = 294

 Score = 25.7 bits (57), Expect = 10.0
 Identities = 16/52 (30%), Positives = 27/52 (51%), Gaps = 4/52 (7%)

Query: 32  KVLEIGSGSGYLTNMISELMNSTGQVIGIEHVPQLVNSSIQNIL--HSNARL 81
            VL++G GSG L   I+ L     +V+G++  P  V ++ +N       A+L
Sbjct: 163 TVLDVGCGSGILA--IAALKLGAKKVVGVDIDPVAVRAAKENAELNGVEAQL 212


  Database: CDD.v3.10
    Posted date:  Mar 20, 2013  7:55 AM
  Number of letters in database: 10,937,602
  Number of sequences in database:  44,354
  
Lambda     K      H
   0.320    0.136    0.391 

Gapped
Lambda     K      H
   0.267   0.0829    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Sequences: 44354
Number of Hits to DB: 7,030,891
Number of extensions: 620047
Number of successful extensions: 689
Number of sequences better than 10.0: 1
Number of HSP's gapped: 687
Number of HSP's successfully gapped: 82
Length of query: 139
Length of database: 10,937,602
Length adjustment: 87
Effective length of query: 52
Effective length of database: 7,078,804
Effective search space: 368097808
Effective search space used: 368097808
Neighboring words threshold: 11
Window for multiple hits: 40
X1: 16 ( 7.4 bits)
X2: 38 (14.6 bits)
X3: 64 (24.7 bits)
S1: 41 (21.8 bits)
S2: 54 (24.4 bits)