Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,853,042 C→T 9.4% A687T (GCT→ACT)  dnaX ← DNA polymerase III subunit gamma

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,853,0420CT 9.4% 146.9 / 5.3 53A687T (GCT→ACT) dnaXDNA polymerase III subunit gamma
Reads supporting (aligned to +/- strand):  ref base C (27/21);  new base T (3/2);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01

TCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCGCTCGGCTCGCTTGCGTGCCGCCGCCT  >  NC_002947/4852898‑4853179
                                                                                                                                                |                                                                                                                                         
tctGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAgg                                                                                                                                       >  3:184315/1‑149 (MQ=255)
     aTTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGTGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGGCTAGTCAGGGccag                                                                                                                                  <  5:147639/149‑4 (MQ=255)
     aTTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTAACCATTACTGGCCCTGACTAGCCAGGGCTTc                                                                                                                                  <  1:29613/149‑1 (MQ=255)
      ttttttCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGGCTAGTCAGGGccagt                                                                                                                                 <  8:81707/149‑5 (MQ=255)
               gCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTc                                                                                                                         >  4:196486/1‑148 (MQ=255)
                          gCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGGCTAGTCAGGGc                                                                                                                                     >  5:22567/1‑125 (MQ=255)
                          gCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTg                                                                                                             <  3:247195/149‑1 (MQ=255)
                            cTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGcc                                                                                                           >  2:25549/1‑149 (MQ=255)
                                ttCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCAc                                                                                                       >  4:147246/1‑149 (MQ=255)
                                        gccggccATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGgcgc                                                                                               >  4:251857/1‑149 (MQ=255)
                                           ggccATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCg                                                                                             >  2:130643/1‑148 (MQ=255)
                                                 gCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACt                                                                                      >  3:30383/1‑149 (MQ=255)
                                                   cACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGaacaac                                                                                       >  5:5947/1‑146 (MQ=255)
                                                   cACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGaacaac                                                                                       <  6:5947/146‑1 (MQ=255)
                                                    aCCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGa                                                                                   >  3:194660/1‑149 (MQ=255)
                                                    aCCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGa                                                                                   >  5:14816/1‑149 (MQ=255)
                                                     ccTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGat                                                                                  <  1:25549/149‑1 (MQ=255)
                                                        ttCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGACGGCCGCATGGTTATCCATTACTGGCCATGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACATTGGCGCCGAACAACTTGatcat                                                                               >  6:27422/1‑149 (MQ=255)
                                                           atcatGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGatc                                                                                 <  8:247035/144‑1 (MQ=255)
                                                           atcatGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGatc                                                                                 >  7:247035/1‑144 (MQ=255)
                                                              atGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGatc                                                                                 >  3:32728/1‑141 (MQ=255)
                                                              atGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGatc                                                                                 <  4:32728/141‑1 (MQ=255)
                                                                      taCCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAAc                                                                  <  3:147246/148‑1 (MQ=255)
                                                                                       gtgATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATGATCTGCTGGATCAACGGCTCCTGCTCGATcgac                                                >  8:210365/1‑149 (MQ=255)
                                                                                            gCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTg                                                                                    <  2:224539/108‑1 (MQ=255)
                                                                                            gCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTg                                                                                    >  1:224539/1‑108 (MQ=255)
                                                                                             cGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATcgacgcgacg                                          <  1:52571/149‑1 (MQ=255)
                                                                                               gCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGc                                        >  7:297994/1‑149 (MQ=255)
                                                                                               gCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGc                                        >  1:152008/1‑149 (MQ=255)
                                                                                                    gCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCct                                   >  4:195718/1‑149 (MQ=255)
                                                                                                      aCCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCctgc                                 <  1:130643/149‑1 (MQ=255)
                                                                                                      aCCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCctgc                                 <  3:196486/149‑1 (MQ=255)
                                                                                                        cGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGAt                                                              <  5:11532/118‑1 (MQ=255)
                                                                                                        cGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGAt                                                              >  6:11532/1‑118 (MQ=255)
                                                                                                              ggccgCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGatcat                                                                               >  7:49729/1‑95 (MQ=255)
                                                                                                              ggccgCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGatcat                                                                               <  8:49729/95‑1 (MQ=255)
                                                                                                              ggccgCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCgctc                         >  6:208423/1‑149 (MQ=255)
                                                                                                               gccgCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCgctcg                        >  8:110925/1‑149 (MQ=255)
                                                                                                                 cgCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCgctcggc                      <  5:208423/149‑1 (MQ=255)
                                                                                                                   cATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCgctcggctc                    >  6:223680/1‑149 (MQ=255)
                                                                                                                    aTGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCgctcggctcg                   >  3:295136/1‑149 (MQ=255)
                                                                                                                        ttATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCGCTCGGCTCGCtta               <  4:30383/149‑2 (MQ=255)
                                                                                                                        ttATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCGCTCGGCTCGCTTg               <  3:195718/149‑1 (MQ=255)
                                                                                                                        ttATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCGCTCGGCTCGCTTg               >  1:216705/1‑149 (MQ=255)
                                                                                                                          aTCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCGCTCGGCTCGCTTGc              <  3:251857/148‑1 (MQ=255)
                                                                                                                            ccATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCtcgatc                                                   >  1:216638/1‑109 (MQ=255)
                                                                                                                            ccATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCtcgatc                                                   <  2:216638/109‑1 (MQ=255)
                                                                                                                               ttACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCCCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCGCTCGGCTCGCTTGCGTgccg        >  8:182316/1‑149 (MQ=255)
                                                                                                                                  cTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCGCTCGGCTCGCTTGCGTgccgcc      <  4:194660/148‑1 (MQ=255)
                                                                                                                                  cTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCGCGCGGCACGCAAGCGAgccgagc     <  1:90850/149‑4 (MQ=255)
                                                                                                                                     gCCCTGGCTAGTCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCGCTCGGCTCGCTTGCGTGCCGCCGCCt  >  4:269287/1‑149 (MQ=255)
                                                                                                                                     gCCCTGGCTAGTCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACCCGACGGCGTCCTGCTGGCGCTCGGCTCGCTTGAGTGCCGCCGCCt  >  7:114323/1‑149 (MQ=255)
                                                                                                                                     gCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCGCTCGGCTCGCTTGCGTgccgccgc    <  8:297994/147‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                         
TCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCGCTCGGCTCGCTTGCGTGCCGCCGCCT  >  NC_002947/4852898‑4853179

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑