Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,853,037 A→G 9.2% S688S (AGT→AGC 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,0370AG 9.2% 142.3 / 5.5 54S688S (AGT→AGCdnaXDNA polymerase III subunit gamma
Reads supporting (aligned to +/- strand):  ref base A (28/21);  new base G (3/2);  total (31/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 = 2.21e-01

TGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCGCTCGGCTCGCTTGCGTGCCGCCGCCT  >  NC_002947/4852891‑4853179
                                                                                                                                                  |                                                                                                                                              
tGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGAc                                                                                                                                               >  2:29613/1‑148 (MQ=255)
       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                                                                                   >  5:14816/1‑149 (MQ=255)
                                                           aCCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGa                                                                                   >  3:194660/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                                                                                    >  1:224539/1‑108 (MQ=255)
                                                                                                   gCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTg                                                                                    <  2:224539/108‑1 (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)
                                                                                                                                                  |                                                                                                                                              
TGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCGCTCGGCTCGCTTGCGTGCCGCCGCCT  >  NC_002947/4852891‑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: