Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 208,530 G→T 29.3% T180T (ACG→ACT 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,5300GT29.3% 143.7 / 56.8 102T180T (ACG→ACTdkgB2,5‑diketo‑D‑gluconate reductase B
Reads supporting (aligned to +/- strand):  ref base G (49/23);  new base T (16/14);  total (65/37)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 102 total aligned reads displayed.

GTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGC  >  minE/208483‑208592
                                               |                                                              
gTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                            >  1:2607407/1‑68 (MQ=255)
gTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                            >  1:3330064/1‑68 (MQ=255)
gTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                            >  1:3145774/1‑68 (MQ=255)
gTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                            >  1:1865849/1‑68 (MQ=255)
gTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                            >  1:199192/1‑68 (MQ=255)
gTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCt                                            >  1:3033103/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:2925393/69‑1 (MQ=255)
   gCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGaaa                                        <  1:698639/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:2916100/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                             >  1:2023740/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                             >  1:2795393/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                             >  1:1800249/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                             >  1:2633057/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                             >  1:3419316/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                             >  1:435871/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                             >  1:60287/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTAt                             >  1:1398986/1‑56 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCa                >  1:1003909/1‑69 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCa                >  1:1484985/1‑69 (MQ=255)
                           cATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCa                >  1:12262/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:2268360/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa            >  1:3216787/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa            >  1:1389836/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa            >  1:250276/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa            >  1:25895/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa            >  1:705270/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa            >  1:2676355/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa            >  1:2611728/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa            >  1:2255164/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa            >  1:2101360/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa            >  1:2060617/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa            >  1:58341/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa            >  1:1498611/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa            >  1:1968145/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa            >  1:1948819/1‑69 (MQ=255)
                               ttACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa            >  1:397562/1‑69 (MQ=255)
                               ttACTTCCTAAATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTa            >  1:691795/1‑69 (MQ=255)
                                tACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGaaa                                        <  1:3099813/40‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa          <  1:3044561/69‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa          <  1:681717/69‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa          <  1:2822245/69‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa          <  1:3117586/69‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa          <  1:3131083/69‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa          <  1:1032590/69‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa          <  1:1088427/69‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa          <  1:2444844/69‑1 (MQ=255)
                                 aCTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa          <  1:2150722/69‑1 (MQ=255)
                                  cTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa          <  1:2476926/68‑1 (MQ=255)
                                  cTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa          <  1:2968192/68‑1 (MQ=255)
                                  cTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa          <  1:1488852/68‑1 (MQ=255)
                                  cTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa          <  1:627091/68‑1 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTc                        >  1:1893364/1‑53 (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:432478/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                       >  1:1625706/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                       >  1:1493334/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                       >  1:183921/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                       >  1:2061705/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                       >  1:648755/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                       >  1:2962501/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                       >  1:2927005/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                       >  1:2898730/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                       >  1:281381/1‑54 (MQ=255)
                                   ttCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCg                       >  1:838319/1‑54 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca        <  1:1402167/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca        <  1:1165364/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:3368111/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca        <  1:2627024/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:3055588/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca        <  1:2109697/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca        <  1:2362900/68‑1 (MQ=255)
                                    tCCTATATGACTCTGGCGTATGGCAATGCACTGAAAGATAAGGTTATTGCTCGTATTGCAGCTAAaca        <  1:2940235/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:2606984/1‑69 (MQ=255)
                                         tatGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGc  >  1:1375014/1‑69 (MQ=255)
                                         tatGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGc  >  1:552153/1‑69 (MQ=255)
                                         tatGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGc  >  1:1663353/1‑69 (MQ=255)
                                           tGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTaaa          <  1:630564/59‑1 (MQ=255)
                                               |                                                              
GTTGCCTGGGCTAAACAGCACGGCATCCATATTACTTCCTATATGACGCTGGCGTATGGTAAGGCCCTGAAAGATGAGGTTATTGCTCGTATCGCAGCTAAACACAATGC  >  minE/208483‑208592

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: