Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 208,542 T→C 28.5% G184G (GGT→GGC 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,5420TC28.5% 105.7 / 46.6 105G184G (GGT→GGCdkgB2,5‑diketo‑D‑gluconate reductase B
Reads supporting (aligned to +/- strand):  ref base T (49/26);  new base C (16/14);  total (65/40)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 105 total aligned reads displayed.

GTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACA  >  minE/208483‑208604
                                                           |                                                              
gTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                                        >  1:3033103/1‑68 (MQ=255)
gTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                                        >  1:1865849/1‑68 (MQ=255)
gTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                                        >  1:199192/1‑68 (MQ=255)
gTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                                        >  1:3145774/1‑68 (MQ=255)
gTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                                        >  1:2607407/1‑68 (MQ=255)
gTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                                        >  1:2456026/1‑68 (MQ=255)
gTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                                        >  1:2560279/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:2023740/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                                         >  1:3419316/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                                         >  1:2633057/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                                         >  1:1800249/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                                         >  1:435871/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                                         >  1:60287/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                                         >  1:2795393/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                                         >  1:2916100/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                                         >  1:1398986/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:1389836/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                        >  1:2268360/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                        >  1:25895/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                        >  1:250276/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                        >  1:2611728/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                        >  1:705270/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                        >  1:3216787/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                        >  1:2101360/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:1498611/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                        >  1:397562/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa                        >  1:2676355/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: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:2444844/69‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                      <  1:3131083/69‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                      <  1:2150722/69‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                      <  1:681717/69‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                      <  1:1088427/69‑1 (MQ=255)
                                  cTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                      <  1:2476926/68‑1 (MQ=255)
                                  cTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                      <  1:627091/68‑1 (MQ=255)
                                  cTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                      <  1:1488852/68‑1 (MQ=255)
                                  cTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa                      <  1:2968192/68‑1 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTc                                    >  1:1893364/1‑53 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                   >  1:3410963/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                   >  1:432478/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                   >  1:648755/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                   >  1:838319/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                   >  1:2584853/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                   >  1:183921/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                   >  1:2061705/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                   >  1:1683750/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                   >  1:1625706/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                   >  1:1493334/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                   >  1:2898730/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                   >  1:2927005/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                                   >  1:2962501/1‑54 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                    <  1:1402167/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                    <  1:1165364/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                    <  1:2940235/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                    <  1:161988/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:2023031/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                    <  1:2109697/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                    <  1:2627024/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                    <  1:3055588/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                    <  1:2362900/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca                    <  1:3368111/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:1375014/1‑69 (MQ=255)
                                         tatGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGc              >  1:552153/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)
                                                           |                                                              
GTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGCGACTCCGGCACA  >  minE/208483‑208604

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: