Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,891,538 T→C 40.8% T15T (ACA→ACG gudP ← predicted D‑glucarate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,891,5380TC40.8% 35.4 / 67.1 93T15T (ACA→ACGgudPpredicted D‑glucarate transporter
Reads supporting (aligned to +/- strand):  ref base T (45/10);  new base C (10/28);  total (55/38)
Fisher's exact test for biased strand distribution p-value = 1.11e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGAAGGATGTGACGATAAACAACATCACCACTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTGAAAA  >  minE/1891487‑1891599
                                                   |                                                             
gTTGAAGGATGTGACGATAAACAACATCACCACTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACt                                              >  1:2481553/1‑69 (MQ=255)
gTTGAAGGATGTGACGATAAACAACATCACCACTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACt                                              >  1:2013314/1‑69 (MQ=255)
gTTGAAGGATGTGACGATAAACAACATCACCACTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACt                                              >  1:187996/1‑69 (MQ=255)
gTTGAAGGATGTGACGATAAACAACATCACCACTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACt                                              >  1:1821374/1‑69 (MQ=255)
gTTGAAGGATGTGACGATAAACAACATCACCACTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACt                                              >  1:506275/1‑69 (MQ=255)
gTTGAAGGATGTGACGATAAACAACATCACCACTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACt                                              >  1:2141515/1‑69 (MQ=255)
gTTGAAGGATGTGACGATAAACAACATCACCACTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACt                                              >  1:2684778/1‑69 (MQ=255)
gTTGAAGGATGTGACGATAAACAACATCACCACTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACt                                              >  1:289318/1‑69 (MQ=255)
gTTGAAGGATGTGACGATAAACAACATCACCACTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACt                                              >  1:1026455/1‑69 (MQ=255)
gTTGAAGGATGTGACGATAAACAACATCACCACTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACt                                              >  1:3054099/1‑69 (MQ=255)
           tgACGATAAACAACATCACCACTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCt                                   >  1:1126065/1‑69 (MQ=255)
           tgACGATAAACAACATCACCACTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCt                                   >  1:1486367/1‑69 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACTAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:2199803/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:377242/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:2963826/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:3110697/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:2045571/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:3085249/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:2059180/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:2082149/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:325877/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:413436/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:2704330/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:3278719/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:2976523/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:2812695/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:231981/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:3345793/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:3186608/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:2745162/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:935645/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:831398/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:749913/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:966143/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:112680/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:1634925/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:1687832/69‑1 (MQ=255)
                 tAAACAACATCACCACTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:1966109/69‑1 (MQ=255)
                          tCACCACTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCa                    >  1:1584833/1‑69 (MQ=255)
                           caccacTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAt                   <  1:1647076/69‑1 (MQ=255)
                           caccacTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAt                   <  1:1347576/69‑1 (MQ=255)
                           caccacTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAt                   <  1:2256248/69‑1 (MQ=255)
                           caccacTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAt                   <  1:2635491/69‑1 (MQ=255)
                           caccacTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAt                   <  1:567023/69‑1 (MQ=255)
                            accacTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAtt                  >  1:2960771/1‑69 (MQ=255)
                            accacTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAtt                  >  1:1018276/1‑69 (MQ=255)
                            accacTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAtt                  >  1:1192988/1‑69 (MQ=255)
                            accacTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAtt                  >  1:2000346/1‑69 (MQ=255)
                            accacTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAtt                  >  1:251764/1‑69 (MQ=255)
                            accacTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAtt                  >  1:2492623/1‑69 (MQ=255)
                            accacTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAtt                  >  1:2365843/1‑69 (MQ=255)
                            accacTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAtt                  >  1:2340724/1‑69 (MQ=255)
                            accacTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAtt                  >  1:1056409/1‑69 (MQ=255)
                            accacTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAtt                  >  1:3425250/1‑69 (MQ=255)
                            accacTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAtt                  >  1:3436194/1‑69 (MQ=255)
                            accacTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAtt                  >  1:663509/1‑69 (MQ=255)
                            accacTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAtt                  >  1:1817630/1‑69 (MQ=255)
                            accacTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAt                   <  1:1375665/68‑1 (MQ=255)
                            accacTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAt                   <  1:1229165/68‑1 (MQ=255)
                            accacTATCCAGTAACGAGCATTTGTGCGTTTTACcac                                                 >  1:2931710/1‑38 (MQ=255)
                                cTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAttgttg              <  1:3175142/69‑1 (MQ=255)
                                cTATCCAGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:2624348/54‑1 (MQ=255)
                                   tCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAttgt                <  1:1333742/64‑1 (MQ=255)
                                   tCCAGTAACGAGCATTTGTGAGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAttgt                <  1:1178660/64‑1 (MQ=255)
                                     cAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTc         >  1:2545031/1‑69 (MQ=255)
                                     cAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTc         >  1:1185933/1‑69 (MQ=255)
                                     cAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTc         >  1:646739/1‑69 (MQ=255)
                                      aGTAACGAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTa                             <  1:975184/48‑1 (MQ=255)
                                         aaCGCGCATTTGTGCGTTTTTCCTCACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTg      >  1:2934679/1‑68 (MQ=255)
                                         aaCGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTg      >  1:3345188/1‑68 (MQ=255)
                                         aaCGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTg      >  1:681128/1‑68 (MQ=255)
                                         aaCGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTg      >  1:2474863/1‑68 (MQ=255)
                                         aaCGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTg      >  1:3033334/1‑68 (MQ=255)
                                         aaCGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTg      >  1:2690632/1‑68 (MQ=255)
                                         aaCGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTg      >  1:2442989/1‑68 (MQ=255)
                                         aaCGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTg      >  1:1050264/1‑68 (MQ=255)
                                         aaCGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTg      >  1:2252991/1‑68 (MQ=255)
                                         aaCGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTg      >  1:2153281/1‑68 (MQ=255)
                                         aaCGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTg      >  1:2121477/1‑68 (MQ=255)
                                         aaCGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTg      >  1:2120142/1‑68 (MQ=255)
                                         aaCGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTg      >  1:211172/1‑68 (MQ=255)
                                         aaCGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTg      >  1:1836943/1‑68 (MQ=255)
                                         aaCGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTg      >  1:1250409/1‑68 (MQ=255)
                                            gAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGc             >  1:3253065/1‑58 (MQ=255)
                                            gAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTGaaaa  >  1:994793/1‑69 (MQ=255)
                                            gAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTGaaaa  >  1:1100954/1‑69 (MQ=255)
                                            gAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTGaaaa  >  1:1414858/1‑69 (MQ=255)
                                            gAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTGaaaa  >  1:146776/1‑69 (MQ=255)
                                            gAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTGaaaa  >  1:697487/1‑69 (MQ=255)
                                            gAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTGaaaa  >  1:1593275/1‑69 (MQ=255)
                                            gAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTGaaaa  >  1:397523/1‑69 (MQ=255)
                                            gAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTGaaaa  >  1:3375/1‑69 (MQ=255)
                                            gAGCATTCGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTGaaaa  >  1:3082169/1‑69 (MQ=255)
                                                   tgtgCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAttgttg              <  1:1469007/50‑1 (MQ=255)
                                                   tgtgCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCAttgttg              <  1:1614287/50‑1 (MQ=255)
                                                   |                                                             
GTTGAAGGATGTGACGATAAACAACATCACCACTATCCAGTAACGAGCATTTGTGCGTTTTTCCACACTGCTCGCAGCCTGACTTAAAGAACTCATTGTTGCACTCCTGAAAA  >  minE/1891487‑1891599

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