Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,872,275 T→C 77.9% C176R (TGC→CGC)  ulaF → L‑ribulose 5‑phosphate 4‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,872,2750TC77.9% 110.5 / 22.0 81C176R (TGC→CGC) ulaFL‑ribulose 5‑phosphate 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (16/2);  new base C (30/33);  total (46/35)
Fisher's exact test for biased strand distribution p-value = 2.38e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATT  >  minE/2872226‑2872338
                                                 |                                                               
ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGATAGATGCCCACGACGCGGTGCATAACGCggtg                                               >  1:74650/1‑68 (MQ=255)
ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                                               >  1:2533167/1‑68 (MQ=255)
ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                                               >  1:984315/1‑68 (MQ=255)
ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                                               >  1:681360/1‑68 (MQ=255)
ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                                               >  1:520097/1‑68 (MQ=255)
ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                                               >  1:506326/1‑68 (MQ=255)
ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                                               >  1:340682/1‑68 (MQ=255)
ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                                               >  1:3059951/1‑68 (MQ=255)
ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                                               >  1:3034215/1‑68 (MQ=255)
ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                                               >  1:2885186/1‑68 (MQ=255)
ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                                               >  1:2370091/1‑68 (MQ=255)
ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                                               >  1:2156777/1‑68 (MQ=255)
ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                                               >  1:1881764/1‑68 (MQ=255)
ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                                               >  1:1708479/1‑68 (MQ=255)
ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                                               >  1:1614997/1‑68 (MQ=255)
 tgtgGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCggtggt                                             <  1:3303790/69‑1 (MQ=255)
  gtggtgTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCggtgg                                              >  1:912962/1‑67 (MQ=255)
  gtggtgTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGgtggtg                                            >  1:2115442/1‑69 (MQ=255)
          tCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaaga                                    >  1:1289280/1‑69 (MQ=255)
          tCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaaga                                    >  1:3089497/1‑69 (MQ=255)
          tCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaaga                                    >  1:2184702/1‑69 (MQ=255)
          tCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaaga                                    >  1:2710180/1‑69 (MQ=255)
          tCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaaga                                    >  1:2077280/1‑69 (MQ=255)
            aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGga                                       >  1:1180736/1‑64 (MQ=255)
            aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGaagaag                                  >  1:540956/1‑69 (MQ=255)
            aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGaagaag                                  >  1:2089168/1‑69 (MQ=255)
            aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGaagaag                                  >  1:2721466/1‑69 (MQ=255)
               aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                                               >  1:1647915/1‑37 (MQ=25)
               aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                                               >  1:3197161/1‑37 (MQ=25)
               aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                                               >  1:3186756/1‑37 (MQ=25)
               aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                                               >  1:357796/1‑37 (MQ=25)
               aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                                               >  1:2956298/1‑37 (MQ=25)
               aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                                               >  1:3020342/1‑37 (MQ=25)
               aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                                               >  1:2468780/1‑37 (MQ=25)
               aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                                               >  1:1326801/1‑37 (MQ=25)
               aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                                               >  1:1401308/1‑37 (MQ=25)
               aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                                               >  1:554411/1‑37 (MQ=25)
               aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                                               >  1:1790361/1‑37 (MQ=25)
               aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                                               >  1:745873/1‑37 (MQ=25)
               aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                                               >  1:1838338/1‑37 (MQ=25)
               aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                                               >  1:2177427/1‑37 (MQ=25)
               aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                                               >  1:2061916/1‑37 (MQ=25)
                   gCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaa                           <  1:902955/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGGCGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:3265477/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:2735877/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:830459/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:77709/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:747051/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:328907/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:1258623/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:3397796/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:1379368/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:627055/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:383387/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:1327710/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:2976774/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:172296/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:1758086/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:1782504/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:1633490/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:2049078/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:2217703/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:2252857/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:1476252/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:2375466/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:244404/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:2601116/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:26156/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:279014/69‑1 (MQ=255)
                     cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:1404523/69‑1 (MQ=255)
                      gTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:284224/68‑1 (MQ=255)
                             tGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:114505/61‑1 (MQ=255)
                              gggtgAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGTAAGAAGTGGCGaaaa                         <  1:2198662/60‑1 (MQ=255)
                                          cccACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:265950/46‑1 (MQ=38)
                                          cccACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:2160215/46‑1 (MQ=38)
                                            cACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCAtt  >  1:2332336/1‑69 (MQ=255)
                                            cACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCAtt  >  1:2180039/1‑69 (MQ=255)
                                            cACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCAtt  >  1:2730504/1‑69 (MQ=255)
                                            cACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCAtt  >  1:75472/1‑69 (MQ=255)
                                            cACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCAtt  >  1:1406114/1‑69 (MQ=255)
                                              cGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                         <  1:2024953/44‑1 (MQ=255)
                                                 |                                                               
TTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATT  >  minE/2872226‑2872338

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: