Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 90,450 G→T M119I (ATG→ATT rsmH → 16S rRNA m(4)C1402 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091390,4500GT98.5% 242.2 / ‑2.7 68M119I (ATG→ATTrsmH16S rRNA m(4)C1402 methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (34/33);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

GTCTCTTCACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGT  >  NC_000913/90403‑90496
                                               |                                              
gTCTCTTCACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATtcg                                              >  1:4048871/1‑50 (MQ=255)
gTCTCTTCACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATtcg                                              >  1:6440734/1‑50 (MQ=255)
gTCTCTTCACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATtcg                                              >  1:2213391/1‑50 (MQ=255)
 tctcTTCACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATtcgc                                             >  1:5261993/1‑50 (MQ=255)
 tctcTTCACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATtcgc                                             >  1:5035534/1‑50 (MQ=255)
  ctctTCACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATtcgc                                             >  1:2010405/1‑49 (MQ=255)
   tctTCACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATTCGCGa                                           >  1:2744463/1‑50 (MQ=255)
   tctTCACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATTCGCGa                                           >  1:2695807/1‑50 (MQ=255)
     ttCACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATTcgcg                                            <  1:1346195/47‑1 (MQ=255)
       cACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATTCGCGATg                                         <  1:4993504/48‑1 (MQ=255)
       cACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATTCGCGATGGt                                       >  1:3039100/1‑50 (MQ=255)
       cACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATTCGCGATGGt                                       >  1:2187848/1‑50 (MQ=255)
        aCCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATTCGCGATGGt                                       <  1:1509619/49‑1 (MQ=255)
        aCCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATTCGCGATGGt                                       <  1:2310967/49‑1 (MQ=255)
        aCCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATTCGCGATGGTc                                      <  1:2281862/50‑1 (MQ=255)
        aCCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATTCGCGATGGTc                                      >  1:3594723/1‑50 (MQ=255)
         ccGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATTCGCGATGGTcc                                     >  1:2012391/1‑50 (MQ=255)
         ccGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATTCGCGATGGTcc                                     <  1:6671932/50‑1 (MQ=255)
           gCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGc                                   >  1:6496012/1‑50 (MQ=255)
           gCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGc                                   <  1:1354849/50‑1 (MQ=255)
           gCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGc                                   >  1:1827354/1‑50 (MQ=255)
            cAACTTGATGATGCTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCt                                  <  1:737189/50‑1 (MQ=255)
            cAACTTGATGATGCTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCt                                  >  1:4780175/1‑50 (MQ=255)
            cAACTTGATGATGCTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCt                                  <  1:3111280/50‑1 (MQ=255)
             aaCTTGATGATGCTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCTg                                 >  1:6017336/1‑50 (MQ=255)
             aaCTTGATGAGGCTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCTg                                 >  1:5433818/1‑50 (MQ=255)
              aCTTGATGATGCTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCTgg                                <  1:6228626/50‑1 (MQ=255)
              aCTTGATGATGCTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCTgg                                <  1:5193788/50‑1 (MQ=255)
                  gatgatGCTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACa                             <  1:1351689/49‑1 (MQ=255)
                  gatgatGCTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACAt                            >  1:3690574/1‑50 (MQ=255)
                   atgatgCTGAACGTGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATg                           <  1:1774258/50‑1 (MQ=255)
                   atgatgCTAAACGTGTCTTTTCCTTTATTCGCGATGGTCCGCTGGACATg                           <  1:5941928/50‑1 (MQ=255)
                    tgatgCTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGc                          >  1:1600512/1‑50 (MQ=255)
                      atgCTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCg                         <  1:3048631/49‑1 (MQ=255)
                      atgCTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCg                         <  1:2844276/49‑1 (MQ=255)
                      atgCTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCg                         >  1:2640680/1‑49 (MQ=255)
                      atgCTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCg                         <  1:2432975/49‑1 (MQ=255)
                       tgCTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTa                       >  1:6507062/1‑50 (MQ=255)
                        gCTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTAt                      >  1:1391920/1‑50 (MQ=255)
                        gCTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTAt                      >  1:5557497/1‑50 (MQ=255)
                         cTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATg                     >  1:2595658/1‑50 (MQ=255)
                         cTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATg                     >  1:2546857/1‑50 (MQ=255)
                         cTGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATg                     <  1:4732722/50‑1 (MQ=255)
                          tGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATg                     >  1:4306483/1‑49 (MQ=255)
                          tGAACGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATg                     >  1:4651136/1‑49 (MQ=255)
                            aaCGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGa                   <  1:6475321/49‑1 (MQ=255)
                            aaCGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGa                   <  1:4038535/49‑1 (MQ=255)
                            aaCGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGa                   <  1:2510490/49‑1 (MQ=255)
                            aaCGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGa                   <  1:2362049/49‑1 (MQ=255)
                            aaCGTGGCTTTTCCTTTATTCGCGACGGTCCGCTGGACATGCGTATGGAc                  >  1:1052926/1‑50 (MQ=255)
                             aCGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGAcc                 <  1:5108836/50‑1 (MQ=255)
                             aCGTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGAcc                 <  1:1684646/50‑1 (MQ=255)
                               gTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGAccc                >  1:2574613/1‑49 (MQ=255)
                               gTGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGACCCa               <  1:853142/50‑1 (MQ=255)
                                tGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGACCCaa              <  1:3783671/50‑1 (MQ=255)
                                tGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGACCCaa              >  1:2553690/1‑50 (MQ=255)
                                 ggCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGACCCAAc             >  1:5393979/1‑50 (MQ=255)
                                 ggCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGACCCAAc             >  1:2373786/1‑50 (MQ=255)
                                 ggCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGACCCAAc             >  1:5790227/1‑50 (MQ=255)
                                  gCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGACCCAAcc            <  1:788735/50‑1 (MQ=255)
                                  gCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGACCCAAcc            >  1:3588261/1‑50 (MQ=255)
                                  gCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGACCCAAcc            <  1:2971905/50‑1 (MQ=255)
                                  gCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGACCCAAc             <  1:361630/49‑1 (MQ=255)
                                   cTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGACCCAAccc           <  1:2909960/50‑1 (MQ=255)
                                    ttttCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCg          >  1:4254028/1‑50 (MQ=255)
                                     tttCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGt         >  1:567971/1‑50 (MQ=255)
                                       tCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTgg       >  1:5340621/1‑50 (MQ=255)
                                       tCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTgg       >  1:2412760/1‑50 (MQ=255)
                                       tCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTg        <  1:1958217/49‑1 (MQ=255)
                                         cTTTATTCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGc     <  1:6174836/50‑1 (MQ=255)
                                         cTTTATTCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGc     <  1:5764001/50‑1 (MQ=255)
                                         cTTTATTCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGc     <  1:4812202/50‑1 (MQ=255)
                                          tttATTCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTggg      >  1:6494589/1‑48 (MQ=255)
                                            tATTCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGt  >  1:5687010/1‑50 (MQ=255)
                                               |                                              
GTCTCTTCACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGT  >  NC_000913/90403‑90496

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