Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 208,575 C→T 33.8% I195I (ATC→ATT 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,5750CT33.8% 101.2 / 67.6 124I195I (ATC→ATTdkgB2,5‑diketo‑D‑gluconate reductase B
Reads supporting (aligned to +/- strand):  ref base C (46/36);  new base T (12/30);  total (58/66)
Fisher's exact test for biased strand distribution p-value = 4.38e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-01

Only 100 of 124 total aligned reads displayed.

CATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATGGGGGA  >  minE/208510‑208631
                                                                 |                                                        
cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCa                                                       >  1:1003909/1‑69 (MQ=255)
cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCa                                                       >  1:12262/1‑69 (MQ=255)
    ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCt                                                    >  1:1326627/1‑68 (MQ=255)
    ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                                   >  1:1389836/1‑69 (MQ=255)
    ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                                   >  1:2268360/1‑69 (MQ=255)
    ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                                   >  1:1498611/1‑69 (MQ=255)
    ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                                   >  1:705270/1‑69 (MQ=255)
    ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                                   >  1:2101360/1‑69 (MQ=255)
    ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                                   >  1:3216787/1‑69 (MQ=255)
    ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                                   >  1:2060617/1‑69 (MQ=255)
    ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                                   >  1:2676355/1‑69 (MQ=255)
    ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                                   >  1:1948819/1‑69 (MQ=255)
    ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                                   >  1:2611728/1‑69 (MQ=255)
    ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                                   >  1:1968145/1‑69 (MQ=255)
    ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                                   >  1:250276/1‑69 (MQ=255)
      aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                                                 <  1:3131083/69‑1 (MQ=255)
      aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                                                 <  1:2150722/69‑1 (MQ=255)
      aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                                                 <  1:3117586/69‑1 (MQ=255)
      aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                                                 <  1:1032590/69‑1 (MQ=255)
      aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                                                 <  1:2822245/69‑1 (MQ=255)
      aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                                                 <  1:681717/69‑1 (MQ=255)
      aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                                                 <  1:2444844/69‑1 (MQ=255)
      aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                                                 <  1:1088427/69‑1 (MQ=255)
       cTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                                                 <  1:627091/68‑1 (MQ=255)
       cTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                                                 <  1:1488852/68‑1 (MQ=255)
       cTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                                                 <  1:2968192/68‑1 (MQ=255)
         tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                                               <  1:2940235/68‑1 (MQ=255)
         tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                                               <  1:2023031/68‑1 (MQ=255)
         tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                                               <  1:1814549/68‑1 (MQ=255)
         tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                                               <  1:1787404/68‑1 (MQ=255)
         tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                                               <  1:1737558/68‑1 (MQ=255)
         tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                                               <  1:1088885/68‑1 (MQ=255)
         tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                                               <  1:2362900/68‑1 (MQ=255)
         tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                                               <  1:3055588/68‑1 (MQ=255)
         tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                                               <  1:3368111/68‑1 (MQ=255)
         tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                                               <  1:1165364/68‑1 (MQ=255)
         tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                                               <  1:1402167/68‑1 (MQ=255)
              tatGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTACACACAATGc                                         >  1:17073/1‑69 (MQ=255)
              tatGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAacac                                              <  1:2597360/64‑1 (MQ=255)
              tatGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATg                                          >  1:2415118/1‑68 (MQ=255)
              tatGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGc                                         >  1:1663353/1‑69 (MQ=255)
              tatGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGc                                         >  1:2606984/1‑69 (MQ=255)
              tatGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGc                                         >  1:552153/1‑69 (MQ=255)
              tatGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGc                                         >  1:1375014/1‑69 (MQ=255)
                tGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                                                 <  1:630564/59‑1 (MQ=255)
                          gTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGcaca                             <  1:1441760/69‑1 (MQ=255)
                          gTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGcaca                             <  1:3212439/69‑1 (MQ=255)
                          gTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGcaca                             <  1:1082987/69‑1 (MQ=255)
                                  aGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTc                     >  1:1674433/1‑69 (MQ=255)
                                  aGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTc                     >  1:2383945/1‑69 (MQ=255)
                                       cTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCg                >  1:2706009/1‑69 (MQ=255)
                                       cTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCg                <  1:2300912/69‑1 (MQ=255)
                                       cTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCg                <  1:1027493/69‑1 (MQ=255)
                                       cTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCg                <  1:91375/69‑1 (MQ=255)
                                       cTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCg                <  1:1080711/69‑1 (MQ=255)
                                       cTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCg                <  1:876586/69‑1 (MQ=255)
                                       cTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCg                <  1:34977/69‑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:2223739/67‑1 (MQ=255)
                                               ttaGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTAt        <  1:2482348/67‑1 (MQ=255)
                                               tGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTAt        <  1:916283/69‑1 (MQ=255)
                                               tGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTAt        <  1:167051/69‑1 (MQ=255)
                                               tGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTAt        <  1:2538913/69‑1 (MQ=255)
                                               tGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTAt        <  1:1987090/69‑1 (MQ=255)
                                               tGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTAt        <  1:2742195/69‑1 (MQ=255)
                                                 aGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATgg      >  1:223263/1‑69 (MQ=255)
                                                 aGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATgg      >  1:430610/1‑69 (MQ=255)
                                                 aGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATgg      >  1:1359603/1‑69 (MQ=255)
                                                 aGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATgg      >  1:640211/1‑69 (MQ=255)
                                                 aGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATgg      >  1:2820492/1‑69 (MQ=255)
                                                 aGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATgg      >  1:739788/1‑69 (MQ=255)
                                                 aGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATgg      >  1:781231/1‑69 (MQ=255)
                                                 aGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATgg      >  1:785242/1‑69 (MQ=255)
                                                 aGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATgg      >  1:2624774/1‑69 (MQ=255)
                                                 aGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATgg      >  1:1005314/1‑69 (MQ=255)
                                                     tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa  >  1:1701040/1‑69 (MQ=255)
                                                     tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa  >  1:1185317/1‑69 (MQ=255)
                                                     tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa  >  1:1449512/1‑69 (MQ=255)
                                                     tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa  >  1:568438/1‑69 (MQ=255)
                                                     tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa  >  1:479746/1‑69 (MQ=255)
                                                     tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa  >  1:1129860/1‑69 (MQ=255)
                                                     tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa  >  1:1815465/1‑69 (MQ=255)
                                                     tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa  >  1:2090795/1‑69 (MQ=255)
                                                     tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa  >  1:3023188/1‑69 (MQ=255)
                                                     tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCTCAAGTGATTCTGGCGTGGGCTATGGGGGa  >  1:2642548/1‑69 (MQ=255)
                                                     tATTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATGGGGGa  >  1:2394173/1‑69 (MQ=255)
                                                      aTTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc           <  1:3371596/59‑1 (MQ=255)
                                                      aTTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc           <  1:1946676/59‑1 (MQ=255)
                                                      aTTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc           <  1:2208407/59‑1 (MQ=255)
                                                      aTTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc           <  1:2288097/59‑1 (MQ=255)
                                                      aTTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc           <  1:2756094/59‑1 (MQ=255)
                                                      aTTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc           <  1:2436012/59‑1 (MQ=255)
                                                      aTTGCTCGTATTGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc           <  1:2667459/59‑1 (MQ=255)
                                                      aTTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc           <  1:25109/59‑1 (MQ=255)
                                                      aTTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc           <  1:357582/59‑1 (MQ=255)
                                                      aTTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc           >  1:1036083/1‑59 (MQ=255)
                                                      aTTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc           >  1:1348831/1‑59 (MQ=255)
                                                      aTTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGc           >  1:1537295/1‑59 (MQ=255)
                                                                 |                                                        
CATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTCTGGCGTGGGCTATGGGGGA  >  minE/208510‑208631

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: