Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 208,602 A→T 27.1% A204A (GCA→GCT dkgB → 2,5‑diketo‑D‑gluconate reductase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE208,6020AT27.1% 95.9 / 37.8 70A204A (GCA→GCTdkgB2,5‑diketo‑D‑gluconate reductase B
Reads supporting (aligned to +/- strand):  ref base A (20/31);  new base T (11/8);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

GTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATGGGGGAAGGTTACTCAGTAATTCCTTCTTCTACTAAACGTAAAAA  >  minE/208536‑208670
                                                                  |                                                                    
gTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGcaca                                                                    <  1:3212439/69‑1 (MQ=255)
gTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGcaca                                                                    <  1:1082987/69‑1 (MQ=255)
gTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGcaca                                                                    <  1:1441760/69‑1 (MQ=255)
        aGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTc                                                            >  1:2383945/1‑69 (MQ=255)
        aGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTc                                                            >  1:1674433/1‑69 (MQ=255)
             cTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCg                                                       >  1:2706009/1‑69 (MQ=255)
             cTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCg                                                       >  1:1261858/1‑69 (MQ=255)
             cTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCg                                                       <  1:1027493/69‑1 (MQ=255)
             cTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCg                                                       <  1:34977/69‑1 (MQ=255)
             cTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCg                                                       <  1:2300912/69‑1 (MQ=255)
             cTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCg                                                       <  1:876586/69‑1 (MQ=255)
             cTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCg                                                       <  1:1080711/69‑1 (MQ=255)
             cTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCg                                                       <  1:91375/69‑1 (MQ=255)
              tGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCg                                                       <  1:864947/68‑1 (MQ=255)
                  aGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTgg                                                    <  1:297059/67‑1 (MQ=255)
                  aGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTgg                                                    <  1:962383/67‑1 (MQ=255)
                  aGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTgg                                                    <  1:1570975/67‑1 (MQ=255)
                  aGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTgg                                                    <  1:2378600/67‑1 (MQ=255)
                  aGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTgg                                                    <  1:2223739/67‑1 (MQ=255)
                  aGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTgg                                                    <  1:1850376/67‑1 (MQ=255)
                  aGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTgg                                                    <  1:2015316/67‑1 (MQ=255)
                     ttaGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTAt                                               <  1:2482348/67‑1 (MQ=255)
                     tGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTAt                                               <  1:1987090/69‑1 (MQ=255)
                     tGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTAt                                               <  1:2538913/69‑1 (MQ=255)
                     tGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTAt                                               <  1:916283/69‑1 (MQ=255)
                     tGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTAt                                               <  1:271161/69‑1 (MQ=255)
                     tGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTAt                                               <  1:2742195/69‑1 (MQ=255)
                     tGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTAt                                               <  1:171100/69‑1 (MQ=255)
                     tGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTAt                                               <  1:3299781/69‑1 (MQ=255)
                     tGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTAt                                               <  1:167051/69‑1 (MQ=255)
                       aGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATgg                                             >  1:2624774/1‑69 (MQ=255)
                       aGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATgg                                             >  1:430610/1‑69 (MQ=255)
                       aGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATgg                                             >  1:1359603/1‑69 (MQ=255)
                       aGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATgg                                             >  1:1005314/1‑69 (MQ=255)
                       aGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATgg                                             >  1:785242/1‑69 (MQ=255)
                       aGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATgg                                             >  1:2820492/1‑69 (MQ=255)
                       aGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATgg                                             >  1:640211/1‑69 (MQ=255)
                       aGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATgg                                             >  1:223263/1‑69 (MQ=255)
                       aGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATgg                                             >  1:739788/1‑69 (MQ=255)
                       aGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATgg                                             >  1:781231/1‑69 (MQ=255)
                           tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa                                         >  1:3023188/1‑69 (MQ=255)
                           tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa                                         >  1:2642548/1‑69 (MQ=255)
                           tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa                                         >  1:1185317/1‑69 (MQ=255)
                           tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa                                         >  1:479746/1‑69 (MQ=255)
                           tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa                                         >  1:568438/1‑69 (MQ=255)
                           tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa                                         >  1:2090795/1‑69 (MQ=255)
                           tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa                                         >  1:1815465/1‑69 (MQ=255)
                           tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa                                         >  1:1701040/1‑69 (MQ=255)
                           tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa                                         >  1:1449512/1‑69 (MQ=255)
                           tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa                                         >  1:1129860/1‑69 (MQ=255)
                           tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa                                         >  1:1376429/1‑69 (MQ=255)
                           tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATGGGGGa                                         >  1:2394173/1‑69 (MQ=255)
                            aTTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc                                                  <  1:3371596/59‑1 (MQ=255)
                            aTTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc                                                  <  1:967182/59‑1 (MQ=255)
                            aTTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc                                                  <  1:2901320/59‑1 (MQ=255)
                            aTTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc                                                  <  1:2756094/59‑1 (MQ=255)
                            aTTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc                                                  <  1:2667459/59‑1 (MQ=255)
                            aTTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc                                                  <  1:2436012/59‑1 (MQ=255)
                            aTTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc                                                  <  1:2288097/59‑1 (MQ=255)
                            aTTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc                                                  <  1:2208407/59‑1 (MQ=255)
                            aTTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc                                                  <  1:1946676/59‑1 (MQ=255)
                            aTTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc                                                  >  1:1036083/1‑59 (MQ=255)
                            aTTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc                                                  >  1:1054642/1‑59 (MQ=255)
                            aTTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc                                                  >  1:1099813/1‑59 (MQ=255)
                            aTTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc                                                  >  1:1348831/1‑59 (MQ=255)
                            aTTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc                                                  >  1:1537295/1‑59 (MQ=255)
                            aTTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc                                                  <  1:25109/59‑1 (MQ=255)
                            aTTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc                                                  <  1:357582/59‑1 (MQ=255)
                                      ttgCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc                                                  <  1:3206116/47‑1 (MQ=255)
                                        gCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATGGGGGAAGGTTACTCAGt                             <  1:2617151/68‑1 (MQ=255)
                                              aaaCACAATGCGACTCCGGCTCAAGTGATTCTGGCg                                                       <  1:1207830/36‑1 (MQ=255)
                                                      tGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATGGGGGAAGGTTACTCAGt                             <  1:3055347/54‑1 (MQ=255)
                                                         gACTCCGGCACAAGTGATTCTGGCGTGGGCTATGGGGGAAGGTTACTCAGt                             <  1:781029/51‑1 (MQ=255)
                                                                  acaAGTGATTCTGGCGTGGGCTATGGGGGAAGGTTACTCAGTAATTCCTTCTTCTACTAAACGTaaaaa  <  1:824086/69‑1 (MQ=255)
                                                                  acaAGTGATTCTGGCGTGGGCTATGGGGGAAGGTTACTCAGTAATTCCTTCTTCTACTAAACGTaaaaa  <  1:1127238/69‑1 (MQ=255)
                                                                  |                                                                    
GTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATGGGGGAAGGTTACTCAGTAATTCCTTCTTCTACTAAACGTAAAAA  >  minE/208536‑208670

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