Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,235,399 G→A L159L (CTC→CTT rlmG ← 23S rRNA m(2)G1835 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,235,3990GA100.0% 182.3 / NA 47L159L (CTC→CTTrlmG23S rRNA m(2)G1835 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/27);  total (20/27)

AGACATTCGCATGGTTGTGGATAGTCCAGTCAGTACCTTCCAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCGAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTGCGCAGCGCACGCAGTTGCTGTTCCAGCAATGCCAGTGTTTTCGGCACTTTGAT  >  NC_000913/3235254‑3235542
                                                                                                                                                 |                                                                                                                                               
agaCATTCGCATGGTTGTGGATAGTCCAGGCAGTACCTTCCAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGCTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCAc                                                                                                                                            >  2:542285/1‑151 (MQ=255)
    aTTCGCATGGTTGTGGATAGTCCAGTCAGTACCTTCCAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACtttt                                                                                                                                        <  1:92458/151‑1 (MQ=255)
      tCGCATGGTTGTGGATAGTCCAGTCAGTACCTTCCAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACttttt                                                                                                                                       >  1:605744/1‑150 (MQ=255)
       cGCATGGTTGTGGATAGTCCAGTCAGTACCTTCCAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCTCTTTTTc                                                                                                                                      <  2:358779/150‑1 (MQ=255)
        gCATGGTTGTGGATAGTCCAGTCAGTACCTTCCAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGa                                                                                                                                    <  2:529079/151‑1 (MQ=255)
               tgtgGATAGTCCAGTCAGTACCTTCCAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGtt                                                                                                                              <  1:311385/150‑1 (MQ=255)
               tgtgGATAGTCCAGTCAGTACCTTCCAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTc                                                                                                                             >  1:306501/1‑151 (MQ=255)
                  ggATAGTCCAGTCAGTACCTTCCAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCa                                                                                                                           <  2:112209/150‑1 (MQ=255)
                   gATAGTCCAGTCAGTACCTTCCAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGc                                                                                                                         >  1:227810/1‑151 (MQ=255)
                    aTAGTCCAGTCAGTACCTTCCAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCg                                                                                                                        >  2:430328/1‑151 (MQ=255)
                       gTCCAGTCAGTACCTTCCAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGGGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTgg                                                                                                                     >  2:167712/1‑151 (MQ=255)
                          cagtcagtACCTTCCAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAg                                                                                                                  <  1:108260/151‑1 (MQ=255)
                                  aCCTTCCAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAAt                                                                                                           <  1:228477/150‑1 (MQ=255)
                                     ttccagcttccagcTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTc                                                                                                        >  2:110257/1‑150 (MQ=255)
                                       ccagcttccagctAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACg                                                                                                     <  2:299230/151‑1 (MQ=255)
                                             tgcagctGACGGTCTGCGGGGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCATGGTGGTGGTTGGGCCAAGCACTGTTTCGCGCAGCTCGAGCGTGGAAGTGTGAATGTCACggggc                                                                                                 <  2:60086/147‑3 (MQ=255)
                                               cagctAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTgg                                                                                             >  1:210302/1‑151 (MQ=255)
                                                    aaCGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCg                                                                                        <  2:429399/151‑1 (MQ=255)
                                                      cGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGc                                                                                      <  1:655971/151‑1 (MQ=255)
                                                       ggTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCa                                                                                     <  1:633585/151‑1 (MQ=255)
                                                                 gcATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACgtg                                                                           >  2:71342/1‑151 (MQ=255)
                                                                 gcATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACgtg                                                                           <  1:333804/151‑1 (MQ=255)
                                                                   aTCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTa                                                                         <  1:504525/151‑1 (MQ=255)
                                                                     cGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATc                                                                       >  1:172203/1‑151 (MQ=255)
                                                                             tGCGGTTGATTGAAAGTGGAATTAATCAGGCGAGCTTATTTCCATGGCAGCGTGGTGGTGGTGGGGCCAAGCACTTTTTCGAACAGTTCCGGCGTGGAAGTGTAAATGTCACGGGCCTTGGCACCGGCAATAACACGTGTAACCGACGTg                                                                <  2:564501/150‑1 (MQ=255)
                                                                                     aTTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTg                                                       >  2:217639/1‑151 (MQ=255)
                                                                                       tGAAAGAGCAATTAATTAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAGGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTgcg                                                     >  2:679510/1‑151 (MQ=255)
                                                                                        gAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTgcgc                                                    >  2:606256/1‑151 (MQ=255)
                                                                                              gCAATTAATCAGGCGCGCCTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGTGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCCCTTTgcgcagcgc                                               <  2:229369/150‑1 (MQ=255)
                                                                                                aaTTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTGCGCAGcgcacg                                            >  1:337726/1‑151 (MQ=255)
                                                                                                 aTTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTGCGCAGcgcacgc                                           >  2:351167/1‑151 (MQ=255)
                                                                                                     aTCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTGCGCAGCGCACGCAGtt                                       <  2:233071/151‑1 (MQ=255)
                                                                                                     aTCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTGCGCAGCGCACGCAGtt                                       >  1:131472/1‑151 (MQ=255)
                                                                                                     aTCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTGCGCAGCGCACGCAGtg                                       <  2:76237/151‑2 (MQ=255)
                                                                                                      tCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTGCGCAGCGCACGCAGTTg                                      >  1:623766/1‑151 (MQ=255)
                                                                                                         ggCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTGCGCAGCGCACGCAGTTGCTg                                   <  2:340350/151‑1 (MQ=255)
                                                                                                          gcgcgcTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTGCGCAGCGCACGCAGTTGCTGt                                  <  2:369689/151‑1 (MQ=255)
                                                                                                          gcgcgcTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTGCGCAGCGCACGCAGTTGCTGt                                  >  2:151986/1‑151 (MQ=255)
                                                                                                           cgcgcTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTGCGCAGCGCACGCAGTTGCTGtt                                 <  2:249127/151‑1 (MQ=255)
                                                                                                              gcTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTGCGCAGCGCACGCAGTTGCTGTTCCa                              <  2:556549/151‑1 (MQ=255)
                                                                                                                ttttttCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTGCGCAGCGCACGCAGTTGCTGTTCCAGc                            >  1:163578/1‑151 (MQ=255)
                                                                                                                    ttCCATGCCAGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTGCGCAGCGCACGCAGTTGCTGTTCCAGCAAt                         <  2:390053/150‑1 (MQ=255)
                                                                                                                             aGCGTGGTGGTGGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTGCGCAGCGCACGCAGTTGCTGTTCCAGCAATGCCAGTGttt               >  2:397820/1‑151 (MQ=255)
                                                                                                                                gtggtggtgGTTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTGCGCAGCGCACGCAGTTGCTGTTCCAGCAATGCCAGTGTTTTCg            <  1:563272/151‑1 (MQ=255)
                                                                                                                                      gtggtTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTGCGCAGCGCACGCAGTTGCTGTTCCAGCAATGCCAGTGTTTTCGGCACtt      <  1:414314/151‑1 (MQ=255)
                                                                                                                                        ggtTGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTGCGCAGCGCACGCAGTTGCTGTTCCAGCAATGCCAGTGTTTTCGGCACTTTg    <  2:216779/151‑1 (MQ=255)
                                                                                                                                           tGGGCCAAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTGCGCAGCGCACGCAGTTGCTGTTCCAGCAATGCCAGTGTTTTCGGCACTTTGAt  <  2:672533/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
AGACATTCGCATGGTTGTGGATAGTCCAGTCAGTACCTTCCAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCGAGCACTTTTTCGAACAGTTCCAGCGTGGAAGTGTGAATGTCACGGGCCTTGGCACCGGCAATAATACGTGTATCCGACGTGACCACTTTGCGCAGCGCACGCAGTTGCTGTTCCAGCAATGCCAGTGTTTTCGGCACTTTGAT  >  NC_000913/3235254‑3235542

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