Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 467,413 G→T 49.3% P751P (CCG→CCT sucA → 2‑oxoglutarate decarboxylase, thiamin‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE467,4130GT49.3% 11.3 / 78.5 73P751P (CCG→CCTsucA2‑oxoglutarate decarboxylase, thiamin‑requiring
Reads supporting (aligned to +/- strand):  ref base G (8/29);  new base T (16/20);  total (24/49)
Fisher's exact test for biased strand distribution p-value = 4.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCCTCCGCGCGTCTGGAACGTTATCTGCAACTTTGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGCTGCGCGGGATGCGTCG  >  minE/467347‑467472
                                                                  |                                                           
ctcctcCGCGCGTCTGGAACGTTATCTGCAACTTTGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTc                                                           <  1:2780230/69‑1 (MQ=255)
ctcctcCGCGCGTCTGGAACGTTATCTGCAACTTTGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTc                                                           <  1:2762540/69‑1 (MQ=255)
ctcctcCGCGCGTCTGGAACGTTATCTGCAACTTTGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTc                                                           <  1:342633/69‑1 (MQ=255)
ctcctcCGCGCGTCTGGAACGTTATCTGCAACTTTGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTc                                                           <  1:69727/69‑1 (MQ=255)
ctcctcCGCGCGTCTGGAACGTTATCTGCAACTTTGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTc                                                           <  1:1835276/69‑1 (MQ=255)
ctcctcCGCGCGTCTGGAACGTTATCTGCAACTTTGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTc                                                           <  1:144043/69‑1 (MQ=255)
ctcctcCGCGCGTCTGGAACGTTATCTGCAACTTTGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTc                                                           <  1:913913/69‑1 (MQ=255)
ctcctcCGCGCGTCTGGAACGTTATCTGCAACTTTGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTc                                                           <  1:825322/69‑1 (MQ=255)
ctcctcCGCGCGTCTGGAACGTTATCTGCAACTTTGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTc                                                           <  1:1797549/69‑1 (MQ=255)
ctcctcCGCGCGTCTGGAACGTTATCTGCAACTTTGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTc                                                           <  1:1720052/69‑1 (MQ=255)
  cctccGCGCGTCTGGAACGTTATCTGCAACTTTGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTCTa                                                         >  1:2682858/1‑69 (MQ=255)
               ggAACGTTATCTGCAACTTTGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTCTACCCCGGCACAgg                                             >  1:1124056/1‑68 (MQ=255)
                  aCGTTATCTGCAACTTTGTGCTTAGCAAAACATGCAGGTTTGCGTACCGTCTACCCCGGCACAGGttt                                          <  1:1573714/68‑1 (MQ=255)
                  aCGTTATCTGCAACTTTGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTCTACCCCGGCACAGGttt                                          <  1:2452352/68‑1 (MQ=255)
                  aCGTTATCTGCAACTTTGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTCTACCCCGGCACAGGttt                                          <  1:2092798/68‑1 (MQ=255)
                  aCGTTATCTGCAACTTTGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTCTACCCCGGCACAGGttt                                          <  1:606839/68‑1 (MQ=255)
                  aCGTTATCTGCAACTTTGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTCTACCCCGGCACAGGttt                                          <  1:2396183/68‑1 (MQ=255)
                  aCGTTATCTGCAACTTTGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTCTACCCCGGCACAGGttt                                          <  1:2615630/68‑1 (MQ=255)
                             aaCTTTGCGCTGAGCAAAACATGCAGGTTTGCGTACCTTCt                                                          >  1:1371379/1‑41 (MQ=38)
                             aaCTTTGCGCTGAGCAAAACATGCAGGTTTGCGTACCTTCTACCCCGGca                                                 >  1:1557774/1‑50 (MQ=255)
                             aaCTTTGCGCTGAGCAAAACATGCAGGTTTGCGTACCTTCTACCCCGGCACAGGTTTACCACAtgctgc                              >  1:944237/1‑69 (MQ=255)
                             aaCTTTGCGCTGAGCAAAACATGCAGGTTTGCGTACCTTCTACCCCGGCACAGGTTTACCACAtgctgc                              >  1:542368/1‑69 (MQ=255)
                             aaCTTTGCGCTGAGCAAAACATGCAGGTTTGCGTACCTTCTACCCCGGCACAGGTTTACCACAtgctgc                              >  1:381503/1‑69 (MQ=255)
                             aaCTTTGCGCTGAGCAAAACATGCAGGTTTGCGTACCTTCTACCCCGGCACAGGTTTACCACAtgctgc                              >  1:2208707/1‑69 (MQ=255)
                             aaCTTTGCGCTGAGCAAAACATGCAGGTTTGCGTACCTTCTACCCCGGCACAGGTTTACCACAtgctgc                              >  1:2190531/1‑69 (MQ=255)
                             aaCTTTGCGCTGAGCAAAACATGCAGGTTTGCGTACCTTCTACCCCGGCACAGGTTTACCACAtgctgc                              >  1:2132761/1‑69 (MQ=255)
                             aaCTTTGCGCTGAGCAAAACATGCAGGTTTGCGTACCTTCTACCCCGGCACAGGTTTACCACAtgctgc                              >  1:281042/1‑69 (MQ=255)
                             aaCTTTGCGCTGAGCAAAACATGCAGGTTTGCGTACCTTCTACCCCGGCACAGGTTTACCACAtgctgc                              >  1:1191346/1‑69 (MQ=255)
                             aaCTTTGCGCTGAGCAAAACATGCAGGTTTGCGTACCTTCTACCCCGGCACAGGTTTACCACAtgctgc                              >  1:1296432/1‑69 (MQ=255)
                             aaCTTTGCGCTGAGCAAAACATGCAGGTTTGCGTACCTTCTACCCCGGCACAGGTTTACCACAtgctgc                              >  1:82370/1‑69 (MQ=255)
                             aaCTTTGCGCTGAGCAAAACATGCAGGTTTGCGTACCTTCTACCCCGGCACAGGTTTACCACAtgctgc                              >  1:1849485/1‑69 (MQ=255)
                             aaCTTTGCGCTGAGCAAAACATGCAGGTTTGCGTACCTTCTACCCCGGCACAGGTTTACCACAtgctgc                              >  1:860076/1‑69 (MQ=255)
                             aaCTTTGCGCTGAGCAAAACATGCAGGTTTGCGTACCTTCTACCCCGGCACAGGTTTACCACAtgctgc                              >  1:1404438/1‑69 (MQ=255)
                             aaCTTTGCGCTGAGCAAAACATGCAGGTTTGCGTACCTTCTACCCCGGCACAGGTTTACCACAtgatgc                              >  1:2996520/1‑69 (MQ=255)
                                 ttGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTCTACCCCGGCACAGGttt                                          <  1:1502340/53‑1 (MQ=255)
                                 ttGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTCTACCCCGGCACAGGttt                                          <  1:923159/53‑1 (MQ=255)
                                 ttGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTCTACCCCGGCACAGGttt                                          <  1:2843440/53‑1 (MQ=255)
                                 ttGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTCTACCCCGGCACAGGttt                                          <  1:1310602/53‑1 (MQ=255)
                                 ttGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTCTACCCCGGCACAGGttt                                          <  1:3357685/53‑1 (MQ=255)
                                 ttGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTCTACCCCGGCACAGGttt                                          <  1:2534878/53‑1 (MQ=255)
                                 ttGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTCTACCCCGGCACAGGttt                                          <  1:3031221/53‑1 (MQ=255)
                                 ttGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTCTACCCCGGCACAGGttt                                          <  1:2206317/53‑1 (MQ=255)
                                 ttGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTCTACCCCGGCACAGGttt                                          <  1:2174866/53‑1 (MQ=255)
                                 ttGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTCTACCCCGGCACAGGttt                                          <  1:21450/53‑1 (MQ=255)
                                        gAGCAAAACATGCAGGTTTGCGTACCTTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGgcg                   <  1:1325109/69‑1 (MQ=255)
                                        gAGCAAAACATGCAGGTTTGCGTACCTTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGgcg                   <  1:3261021/69‑1 (MQ=255)
                                        gAGCAAAACATGCAGGTTTGCGTACCTTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGgcg                   <  1:1960649/69‑1 (MQ=255)
                                            aaaaCATGCAGGTTTGCGTACCGTCTACCCCGGCACAGGttt                                          <  1:3341679/42‑1 (MQ=255)
                                             aaaCATGCAGGTTTGCGTACCGTCTACCCCGGCACAGGttt                                          <  1:241196/41‑1 (MQ=255)
                                                  tGCAGGTTTGCGTACCGTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGCTGCGCGGGa         <  1:2007309/69‑1 (MQ=255)
                                                        tttGCGTACCTTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGTTGCGCGGGATGcgtc   <  1:1098447/69‑1 (MQ=255)
                                                        tttGCGTACCTTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGTTGCGCGGGATGcgtc   <  1:3364312/69‑1 (MQ=255)
                                                        tttGCGTACCTTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGTTGCGCGGGATGcgtc   <  1:3368240/69‑1 (MQ=255)
                                                        tttGCGTACCTTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGTTGCGCGGGATGcgtc   <  1:895601/69‑1 (MQ=255)
                                                        tttGCGTACCTTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGTTGCGCGGGATGcgtc   <  1:2994439/69‑1 (MQ=255)
                                                        tttGCGTACCTTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGTTGCGCGGGATGcgtc   <  1:1044961/69‑1 (MQ=255)
                                                        tttGCGTACCTTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGTTGCGCGGGATGcgtc   <  1:1270005/69‑1 (MQ=255)
                                                        tttGCGTACCTTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGTTGCGCGGGATGcgtc   <  1:1451415/69‑1 (MQ=255)
                                                        tttGCGTACCTTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGTTGCGCGGGATGcgtc   <  1:1630346/69‑1 (MQ=255)
                                                        tttGCGTACCTTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGTTGCGCGGGATGcgtc   <  1:1905786/69‑1 (MQ=255)
                                                        tttGCGTACCTTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGTTGCGCGGGATGcgtc   <  1:2279083/69‑1 (MQ=255)
                                                        tttGCGTACCTTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGTTGCGCGGGATGcgtc   <  1:2511353/69‑1 (MQ=255)
                                                        tttGCGTACCTTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGTTGCGCGGGATGcgtc   <  1:2599914/69‑1 (MQ=255)
                                                        tttGCGTACCTTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGTTGCGCGGGATGcgtc   <  1:3288783/69‑1 (MQ=255)
                                                        tttGCGTACCTTCTACCCCGGCACAGGTTTACCACATCCTGCGTCGTCAGGCGTTGCGCGGGATGcgtc   <  1:1449506/69‑1 (MQ=255)
                                                        tttGCGTACCGTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGCTGCGCGGGATGcgtc   >  1:521502/1‑69 (MQ=255)
                                                        tttGCGTACCGTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGCTGCGCGGGATGcgtc   >  1:1823615/1‑69 (MQ=255)
                                                        tttGCGTACCGTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGCTGCGCGGGATGcgtc   >  1:3362786/1‑69 (MQ=255)
                                                        tttGCGTACCGTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGCTGCGCGGGATGcgt    >  1:1915958/1‑68 (MQ=255)
                                                         ttGCGTACCTTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGTTGCGCGGGATGcgtc   <  1:264877/68‑1 (MQ=255)
                                                         ttGCGTACCTTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGTTGCGCGGGATGcgtc   <  1:268792/68‑1 (MQ=255)
                                                         ttGCGTACCGTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGCTGCGCGGGATGcgtcg  >  1:949974/1‑69 (MQ=255)
                                                         atGCGTACCGTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGCTGCGCGGGATGcgtcg  >  1:281254/2‑69 (MQ=255)
                                                                  |                                                           
CTCCTCCGCGCGTCTGGAACGTTATCTGCAACTTTGTGCTGAGCAAAACATGCAGGTTTGCGTACCGTCTACCCCGGCACAGGTTTACCACATGCTGCGTCGTCAGGCGCTGCGCGGGATGCGTCG  >  minE/467347‑467472

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: