Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 208,545 G→T 27.9% K185N (AAG→AAT 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,5450GT27.9% 158.8 / 56.3 107K185N (AAG→AATdkgB2,5‑diketo‑D‑gluconate reductase B
Reads supporting (aligned to +/- strand):  ref base G (51/26);  new base T (16/14);  total (67/40)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 107 total aligned reads displayed.

GTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTC  >  minE/208483‑208612
                                                              |                                                                   
gTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                                                >  1:2607407/1‑68 (MQ=255)
gTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                                                >  1:2560279/1‑68 (MQ=255)
gTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                                                >  1:2456026/1‑68 (MQ=255)
gTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                                                >  1:3033103/1‑68 (MQ=255)
gTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                                                >  1:199192/1‑68 (MQ=255)
gTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                                                >  1:3145774/1‑68 (MQ=255)
gTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                                                >  1:3330064/1‑68 (MQ=255)
gTTGCCTGGGCTAAACAGCACGGCATCAATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                                                >  1:2432364/1‑68 (MQ=255)
   gCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGaaa                                                            <  1:2543749/69‑1 (MQ=255)
   gCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGaaa                                                            <  1:698639/69‑1 (MQ=255)
   gCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGaaa                                                            <  1:2925393/69‑1 (MQ=255)
   gCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGaaa                                                            <  1:3363367/69‑1 (MQ=255)
           tAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGt                                                    <  1:3038428/69‑1 (MQ=255)
           tAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGt                                                    <  1:1732815/69‑1 (MQ=255)
                    cGGCTTCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGaaa                                                            <  1:3122265/52‑1 (MQ=255)
                       cATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTAt                                        >  1:783804/1‑69 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                                                 >  1:435871/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                                                 >  1:2795393/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                                                 >  1:2916100/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                                                 >  1:2023740/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                                                 >  1:1398986/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                                                 >  1:1800249/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                                                 >  1:60287/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                                                 >  1:3419316/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCa                                    >  1:1003909/1‑69 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCa                                    >  1:12262/1‑69 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCa                                    >  1:1484985/1‑69 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCATGAAAGATGAGGTTAt                                                 >  1:567644/1‑56 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCg                                           >  1:3277643/1‑58 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCt                                 >  1:1326627/1‑68 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                >  1:250276/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                >  1:2676355/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                >  1:705270/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                >  1:397562/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                >  1:1498611/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                >  1:2611728/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                >  1:25895/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                >  1:2268360/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                >  1:2255164/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                >  1:3216787/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                >  1:2101360/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                >  1:1389836/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                >  1:2060617/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                >  1:1968145/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                >  1:1948819/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                >  1:58341/1‑69 (MQ=255)
                               ttACTTCCTAAATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                                >  1:691795/1‑69 (MQ=255)
                                tACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGaaa                                                            <  1:3099813/40‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                              <  1:681717/69‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                              <  1:2822245/69‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                              <  1:3044561/69‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                              <  1:3117586/69‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                              <  1:1032590/69‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                              <  1:3131083/69‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                              <  1:2444844/69‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                              <  1:1088427/69‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                              <  1:2150722/69‑1 (MQ=255)
                                  cTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                              <  1:1488852/68‑1 (MQ=255)
                                  cTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                              <  1:627091/68‑1 (MQ=255)
                                  cTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                              <  1:2476926/68‑1 (MQ=255)
                                  cTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                              <  1:2968192/68‑1 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTc                                            >  1:1893364/1‑53 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                           >  1:838319/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                           >  1:1625706/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                           >  1:648755/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                           >  1:3410963/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                           >  1:1679642/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                           >  1:1683750/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                           >  1:183921/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                           >  1:2061705/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                           >  1:2962501/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                           >  1:2584853/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                           >  1:281381/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                           >  1:2898730/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                           >  1:2927005/1‑54 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                            <  1:2362900/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                            <  1:1402167/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                            <  1:161988/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                            <  1:2627024/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                            <  1:3368111/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                            <  1:1165364/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                            <  1:2109697/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                            <  1:1088885/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                            <  1:1737558/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                            <  1:1787404/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                            <  1:1814549/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                            <  1:2940235/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                            <  1:2023031/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:552153/1‑69 (MQ=255)
                                         tatGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGc                      >  1:1663353/1‑69 (MQ=255)
                                         tatGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGc                      >  1:1375014/1‑69 (MQ=255)
                                         tatGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGc                      >  1:2606984/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)
                                                              |                                                                   
GTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACAAGTGATTC  >  minE/208483‑208612

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: