Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,493,694 T→G 22.4% intergenic (+76/‑299) 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,6940TG22.4% 127.7 / 38.0 67intergenic (+76/‑299)gyrB/ibpADNA gyrase, subunit B/heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base T (28/24);  new base G (7/8);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 7.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCTATAAGAAAAGGGCG  >  minE/2493630‑2493759
                                                                |                                                                 
aTGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGAt                                                             <  1:2327127/71‑1 (MQ=255)
 tGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTc                                                                  <  1:3570730/65‑1 (MQ=255)
 tGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTc                                                                  <  1:633648/65‑1 (MQ=255)
    gAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTg                                                         <  1:4009331/71‑1 (MQ=255)
    gAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTg                                                         <  1:4002687/71‑1 (MQ=255)
         tGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTc                                                                  <  1:5127138/57‑1 (MQ=255)
         tGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTc                                                                  <  1:3892662/57‑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:2305677/71‑1 (MQ=255)
                    cTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCg                                         <  1:2182818/71‑1 (MQ=255)
                    cTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCg                                         >  1:3665076/1‑71 (MQ=255)
                     tACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGt                                        >  1:3218808/1‑71 (MQ=255)
                      ccGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGc                                                    >  1:2676404/2‑58 (MQ=14)
                       cGCTTATCAGGCCTACGAAAATTCTGCAATTTATTGAATTTGCACGATTTTGTa                                                       <  1:983921/54‑1 (MQ=11)
                       cGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCa                                                                 >  1:3049601/1‑44 (MQ=11)
                       cGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTa                                                       <  1:1981738/54‑1 (MQ=14)
                            aTCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGATAAGGCGTTAAcg                                   <  1:3397947/69‑1 (MQ=14)
                              cAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTg                                                         >  1:5193231/1‑45 (MQ=14)
                              cAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTg                                                         >  1:3333938/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:679215/37‑1 (MQ=255)
                                    tACGAAAATCTTGCAATGTATTGAATTTTCATGAttt                                                           <  1:4736070/37‑1 (MQ=255)
                                    tACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCAcgccgc                               <  1:4336483/65‑1 (MQ=255)
                                    tACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCAc                                    <  1:3905775/60‑1 (MQ=255)
                                      cGAAAATCTTGCAATGTATTGAATTTTCATGAtttt                                                          >  1:87395/1‑36 (MQ=255)
                                      cGAAAATCTTGCAATGTATTGAATTTTCATGAtttt                                                          >  1:1183731/1‑36 (MQ=255)
                                      cGAAAATCTTGCAATGTATTGAATTTTCATGAtttt                                                          >  1:2569587/1‑36 (MQ=255)
                                      cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCAcgc                                  >  1:4880723/1‑60 (MQ=255)
                                      cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCAcgc                                  >  1:4238418/1‑60 (MQ=255)
                                      cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCa                              <  1:2333479/64‑1 (MQ=255)
                                      cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                        >  1:3010509/1‑70 (MQ=255)
                                      cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                        >  1:4910070/1‑70 (MQ=255)
                                      cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                        >  1:2170843/1‑70 (MQ=255)
                                      cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                        >  1:3104643/1‑70 (MQ=255)
                                        aaaaTCTTGCAATGTATTGAATTTTCATGATTTTGTa                                                       >  1:2374485/1‑37 (MQ=255)
                                        aaaaTCTTGCAATGTATTGAATTTTCATGATTTTGTa                                                       >  1:1563655/1‑37 (MQ=255)
                                        aaaaTCTTGCAATGTATTGAATTTTCATGATTTTGTa                                                       >  1:2625088/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:4330642/1‑68 (MQ=255)
                                              ttGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaagaag                  >  1:2946024/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)
                                                                |                                                                 
ATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAGCTATAAGAAAAGGGCG  >  minE/2493630‑2493759

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: