Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,493,697 T→C 19.4% intergenic (+79/‑296) gyrB → / → ibpA DNA gyrase, subunit B/heat shock chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,493,6970TC19.4% 152.9 / 26.4 67intergenic (+79/‑296)gyrB/ibpADNA gyrase, subunit B/heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base T (32/22);  new base C (6/7);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 5.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

ATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCTATAAGAAAAGGGCGGAGATCA  >  minE/2493630‑2493766
                                                                   |                                                                     
aTGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGAt                                                                    <  1:2327127/71‑1 (MQ=255)
    gAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTg                                                                <  1:4009331/71‑1 (MQ=255)
    gAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTg                                                                <  1:4002687/71‑1 (MQ=255)
           ccTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTg                                                                <  1:1247188/64‑1 (MQ=255)
               aTGCGCTCCGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGAt                                                      <  1:254560/70‑1 (MQ=14)
               aTGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGAt                                                      <  1:4306832/70‑1 (MQ=255)
                tGCGCTCCGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGAtt                                                                   >  1:1173598/1‑56 (MQ=14)
                tGCGCTCCGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGAtt                                                                   >  1:4838595/1‑56 (MQ=14)
                tGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGAtttt                                                                 >  1:3730010/1‑58 (MQ=255)
                tGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGAtttt                                                                 >  1:3610521/1‑58 (MQ=255)
                 gcgcTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATaa                                                    <  1:3365579/70‑1 (MQ=255)
                    cTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCg                                                >  1:3665076/1‑71 (MQ=255)
                    cTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCg                                                <  1:2182818/71‑1 (MQ=255)
                    cTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCg                                                <  1:2305677/71‑1 (MQ=255)
                     tACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGt                                               >  1:3218808/1‑71 (MQ=255)
                      ccGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGc                                                           >  1:2676404/2‑58 (MQ=14)
                       cGCTTATCAGGCCTACGAAAATTCTGCAATTTATTGAATTTGCACGATTTTGTa                                                              <  1:983921/54‑1 (MQ=11)
                       cGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTa                                                              <  1:1981738/54‑1 (MQ=14)
                            aTCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGATAAGGCGTTAAcg                                          <  1:3397947/69‑1 (MQ=14)
                              cAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTg                                                                >  1:3333938/1‑45 (MQ=14)
                              cAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTg                                                                >  1:5193231/1‑45 (MQ=14)
                               aGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGATAAGGCGTTAACGCCGCa                                     <  1:4731642/71‑1 (MQ=14)
                               aGGCCTACGAAAATCTTGCAATGTATTGAATTTTTATGATTTTGTAGGc                                                           >  1:4527492/1‑49 (MQ=255)
                                   cTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGATAAGGCGtt                                              <  1:3876547/58‑1 (MQ=11)
                                   cTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGATAAGGCGTTAACGCCGCATcc                                  <  1:1039345/70‑1 (MQ=11)
                                   cTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCCg                                 >  1:2040621/1‑71 (MQ=11)
                                   cTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTa                                                              >  1:3301491/1‑42 (MQ=255)
                                   cTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                  <  1:5466208/70‑1 (MQ=255)
                                    tACGAAAATCTTGCAATGTATTGAATTTTCATGAttt                                                                  <  1:4736070/37‑1 (MQ=255)
                                    tACGAAAATCTTGCAATGTATTGAATTTTCATGAttt                                                                  <  1:679215/37‑1 (MQ=255)
                                    tACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCAcgccgc                                      <  1:4336483/65‑1 (MQ=255)
                                    tACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCAc                                           <  1:3905775/60‑1 (MQ=255)
                                      cGAAAATCTTGCAATGTATTGAATTTTCATGAtttt                                                                 >  1:1183731/1‑36 (MQ=255)
                                      cGAAAATCTTGCAATGTATTGAATTTTCATGAtttt                                                                 >  1:87395/1‑36 (MQ=255)
                                      cGAAAATCTTGCAATGTATTGAATTTTCATGAtttt                                                                 >  1:2569587/1‑36 (MQ=255)
                                      cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCAcgc                                         >  1:4238418/1‑60 (MQ=255)
                                      cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCAcgc                                         >  1:4880723/1‑60 (MQ=255)
                                      cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCa                                     <  1:2333479/64‑1 (MQ=255)
                                      cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                               >  1:4910070/1‑70 (MQ=255)
                                      cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                               >  1:2170843/1‑70 (MQ=255)
                                      cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                               >  1:3010509/1‑70 (MQ=255)
                                      cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                               >  1:3104643/1‑70 (MQ=255)
                                        aaaaTCTTGCAATGTATTGAATTTTCATGATTTTGTa                                                              >  1:2625088/1‑37 (MQ=255)
                                        aaaaTCTTGCAATGTATTGAATTTTCATGATTTTGTa                                                              >  1:2374485/1‑37 (MQ=255)
                                        aaaaTCTTGCAATGTATTGAATTTTCATGATTTTGTa                                                              >  1:1563655/1‑37 (MQ=255)
                                         aaaTCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaaga                           <  1:1621509/71‑1 (MQ=255)
                                           aTCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCg                                                <  1:3741076/48‑1 (MQ=255)
                                           aTCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCAcgcc                                        >  1:4676899/1‑56 (MQ=255)
                                             cTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGtt                                              >  1:2264722/1‑48 (MQ=255)
                                              ttGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaagaag                         >  1:2946024/1‑68 (MQ=255)
                                              ttGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaagaag                         >  1:4330642/1‑68 (MQ=255)
                                              ttGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaagaa                          >  1:462593/1‑67 (MQ=255)
                                                 cAATGTATTGAATTTGCATGATTTTGTAGGCCGGATAAGGCGTTCAc                                           <  1:3526024/47‑1 (MQ=255)
                                                    tGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCTATAAGaa                <  1:3060236/71‑1 (MQ=255)
                                                      tATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCTATAAGaa                >  1:3563597/1‑69 (MQ=255)
                                                        ttGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCAc                                           <  1:3786309/40‑1 (MQ=40)
                                                        ttGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCTATAAGaaa               <  1:244967/68‑1 (MQ=255)
                                                        ttGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCTATAAGAAAAgg            >  1:2694043/1‑71 (MQ=255)
                                                          gAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCTATAAGaaa               >  1:426694/1‑66 (MQ=255)
                                                          gAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCTATAAGAAAAgg            >  1:3792377/1‑69 (MQ=255)
                                                          gAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCTATAAGAAAAgg            >  1:3051503/1‑69 (MQ=255)
                                                           aaTTTTCATGATTTTGTAGGCCGGATAAGGCGTTCAc                                           <  1:2969355/37‑1 (MQ=40)
                                                               ttCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCTATAAGaa                >  1:331819/1‑60 (MQ=255)
                                                               ttCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCTATAAGAAAAGGGCg         >  1:4111671/1‑67 (MQ=255)
                                                               ttCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCTATAAGAAAAGGGCg         >  1:513090/1‑67 (MQ=255)
                                                                 cATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCTATAAGAAAAGGGCGGAGat    <  1:3827537/70‑1 (MQ=255)
                                                                  aTGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCTATAAGAAAAGGGCGGAGatca  >  1:2543330/1‑71 (MQ=255)
                                                                   tGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCTATAAGaaa               >  1:1319425/1‑57 (MQ=255)
                                                                   |                                                                     
ATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCTATAAGAAAAGGGCGGAGATCA  >  minE/2493630‑2493766

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