Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,887,033 G→A G69D (GGC→GAC)  mnmE → tRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,887,0330GA100.0% 157.0 / NA 43G69D (GGC→GAC) mnmEtRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/18);  total (25/18)

GGCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGGCCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGTTTTCCGAACGCGCGTTTCTT  >  NC_000913/3886885‑3887178
                                                                                                                                                    |                                                                                                                                                 
ggCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGAcc                                                                                                                                                 >  1:345927/1‑151 (MQ=255)
  cATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCg                                                                                                                                               <  2:55111/151‑1 (MQ=255)
         cgcATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGt                                                                                                                                        >  2:239514/1‑151 (MQ=255)
                 cGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGc                                                                                                                                <  1:145458/151‑1 (MQ=255)
                      tCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAg                                                                                                                            <  1:556141/150‑1 (MQ=255)
                       cAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGa                                                                                                                              >  1:850799/1‑147 (MQ=255)
                          aGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGAtgtg                                                                                                                       >  1:599433/1‑151 (MQ=255)
                            cccGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCt                                                                                                                     >  2:11167/1‑151 (MQ=255)
                               gTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTg                                                                                                                    <  2:836529/149‑1 (MQ=255)
                               gTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTg                                                                                                                    >  1:836471/1‑149 (MQ=255)
                                  aaGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACt                                                                                                               <  2:548797/151‑1 (MQ=255)
                                   aGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTg                                                                                                              <  2:427986/151‑1 (MQ=255)
                                     ttGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGc                                                                                                             >  1:759750/1‑150 (MQ=255)
                                       gCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAg                                                                                                          <  2:405026/151‑1 (MQ=255)
                                         cGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGt                                                                                                        <  1:47824/151‑1 (MQ=255)
                                          gAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTc                                                                                                       >  1:283563/1‑151 (MQ=255)
                                            aaCCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCAt                                                                                                     >  1:152982/1‑151 (MQ=255)
                                                gTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCg                                                                                                 <  1:829426/151‑1 (MQ=255)
                                                             tGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATcc                                                                                     >  1:586318/1‑150 (MQ=255)
                                                                  aaGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGAcc                                                                               >  2:273891/1‑151 (MQ=255)
                                                                    gCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACctg                                                                             <  1:443888/151‑1 (MQ=255)
                                                                        cgcTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGGGATCCTAGACCTGCTGt                                                                         >  2:626635/1‑151 (MQ=255)
                                                                            aCGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGCCCTGCTGTTa                                                                       <  1:225113/149‑1 (MQ=255)
                                                                              gCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTaaaa                                                                    >  2:607115/1‑150 (MQ=255)
                                                                                cGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGAGATCCTCGACCTGCTGTTAAAACGc                                                                 >  1:332226/1‑151 (MQ=255)
                                                                                         tCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCt                                                            <  1:307325/147‑1 (MQ=255)
                                                                                         tCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGAcc                                                        >  1:14567/1‑151 (MQ=255)
                                                                                             tttAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAAGTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGACCGGGGAACCTCCAACTGCTGTTTAAACGAACTATGACCAtt                                                     >  1:188120/1‑150 (MQ=255)
                                                                                              ttAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTcc                                                   >  2:120422/1‑151 (MQ=255)
                                                                                                aaaGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTccc                                                  <  1:484573/150‑1 (MQ=255)
                                                                                                 aaGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGGGATTCTCGACCTGCTGTTAAAACGCATTCTGACCATTTCCgg                                                >  1:129411/1‑151 (MQ=255)
                                                                                                    aCGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCt                                             <  2:383535/151‑1 (MQ=255)
                                                                                                         gACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCgg                                         >  2:129123/1‑150 (MQ=255)
                                                                                                                    gCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGGGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGccc                              >  1:773949/1‑150 (MQ=255)
                                                                                                                       cGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTgg                           >  2:421918/1‑150 (MQ=255)
                                                                                                                          tCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTg                         >  2:486683/1‑149 (MQ=255)
                                                                                                                            aGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGtt                     <  1:392778/151‑1 (MQ=255)
                                                                                                                                gATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGGCCGGTGCTCCTCGACCTGCTGTTAAAACGCATTCTGACCCTTCCCGGCCTGCGGATTGCTCGCCCCGGTTGGGTTTcc                 >  1:827306/1‑151 (MQ=255)
                                                                                                                                    gcgcTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGTTTTCCGAAc             <  2:606039/151‑1 (MQ=255)
                                                                                                                                      gcTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGTTTTCCGAAcgc           <  2:425535/151‑1 (MQ=255)
                                                                                                                                         aTGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGTTTTCCGAAcgcgcg        <  1:667307/151‑1 (MQ=255)
                                                                                                                                         aTGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCCGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGGGAGTTTTCCGAAcgcgcg        >  1:752963/1‑151 (MQ=255)
                                                                                                                                                 cTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGGGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGTTTTCCGAACGCGCGTTTCtt  >  1:413452/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGGCCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGTTTTCCGAACGCGCGTTTCTT  >  NC_000913/3886885‑3887178

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑