Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,246,542 T→C 20.9% I85I (ATT→ATC thiC → thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,246,5420TC20.9% 133.5 / 25.8 95I85I (ATT→ATCthiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base T (54/21);  new base C (8/12);  total (62/33)
Fisher's exact test for biased strand distribution p-value = 1.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCGCGGCG  >  minE/2246476‑2246603
                                                                  |                                                             
gCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTg                                                              >  1:3337803/1‑68 (MQ=255)
gCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTg                                                              >  1:2683070/1‑68 (MQ=255)
gCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTg                                                              >  1:2248704/1‑68 (MQ=255)
gCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTg                                                              >  1:1569090/1‑68 (MQ=255)
     aCGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa                                                        <  1:1079652/69‑1 (MQ=255)
     aCGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa                                                        <  1:3398720/69‑1 (MQ=255)
     aCGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa                                                        <  1:3140486/69‑1 (MQ=255)
     aCGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa                                                        <  1:2627979/69‑1 (MQ=255)
     aCGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa                                                        <  1:2141788/69‑1 (MQ=255)
     aCGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa                                                        <  1:207332/69‑1 (MQ=255)
     aCGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa                                                        <  1:1928998/69‑1 (MQ=255)
     aCGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa                                                        <  1:171432/69‑1 (MQ=255)
     aCGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa                                                        <  1:1291207/69‑1 (MQ=255)
     aCGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa                                                        <  1:1330964/69‑1 (MQ=255)
      cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTaa                                                       >  1:2671465/1‑69 (MQ=255)
      cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTaa                                                       >  1:232537/1‑69 (MQ=255)
      cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTaa                                                       >  1:1773763/1‑69 (MQ=255)
      cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTaa                                                       >  1:1771395/1‑69 (MQ=255)
      cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTaa                                                       >  1:3079204/1‑69 (MQ=255)
      cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTaa                                                       >  1:1211756/1‑69 (MQ=255)
      cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTaa                                                       >  1:1510760/1‑69 (MQ=255)
      cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTaa                                                       >  1:1535707/1‑69 (MQ=255)
      cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa                                                        <  1:534856/68‑1 (MQ=255)
      cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGTCCGTATGGCGATCCGCAGATCGCCATTa                                                        <  1:1177248/68‑1 (MQ=255)
      cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTaa                                                       >  1:837904/1‑69 (MQ=255)
      cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTaa                                                       >  1:280636/1‑69 (MQ=255)
      cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTaa                                                       >  1:2810329/1‑69 (MQ=255)
      cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTaa                                                       >  1:298021/1‑69 (MQ=255)
      cGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTaa                                                       >  1:1939589/1‑69 (MQ=255)
                   gCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAggg                                          >  1:1319477/1‑69 (MQ=255)
                   gCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAggg                                          >  1:2567655/1‑69 (MQ=255)
                   gCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAggg                                          >  1:2448481/1‑69 (MQ=255)
                   gCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAggg                                          >  1:1391099/1‑69 (MQ=255)
                   gCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAggg                                          >  1:1841543/1‑69 (MQ=255)
                   gCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAggg                                          >  1:858596/1‑69 (MQ=255)
                               tACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTa                              >  1:2690262/1‑69 (MQ=255)
                               tACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTa                              >  1:262835/1‑69 (MQ=255)
                               tACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTa                              >  1:3349629/1‑69 (MQ=255)
                               tACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTa                              >  1:2982035/1‑69 (MQ=255)
                               tACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTa                              >  1:1633548/1‑69 (MQ=255)
                               tACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTa                              >  1:3097816/1‑69 (MQ=255)
                               tACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTa                              >  1:3099152/1‑69 (MQ=255)
                                                  aTGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCg              >  1:2693108/1‑66 (MQ=255)
                                                  aTGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATg           >  1:2080443/1‑69 (MQ=255)
                                                  aTGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATg           >  1:492246/1‑69 (MQ=255)
                                                  aTGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATg           >  1:488356/1‑69 (MQ=255)
                                                   tGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATgc          >  1:623070/1‑69 (MQ=255)
                                                     gTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATgcgc        >  1:2819381/1‑69 (MQ=255)
                                                     gTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATgcgc        >  1:3030886/1‑69 (MQ=255)
                                                     gTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATgcgc        >  1:998535/1‑69 (MQ=255)
                                                     gTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATgcgc        >  1:687174/1‑69 (MQ=255)
                                                     gTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATgcgc        >  1:671020/1‑69 (MQ=255)
                                                     gTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATgcgc        >  1:2043356/1‑69 (MQ=255)
                                                     gTGATCCGCAGATTGCCATTAACGTGCAGCAACGGCTGGCAAAACTACGCCAGCCGTGGATCGATgcgc        >  1:255138/1‑69 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTAc                             >  1:2637647/1‑43 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTAc                             >  1:1038308/1‑43 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTAc                             >  1:3151838/1‑43 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTAc                             >  1:1478604/1‑43 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTAc                             >  1:1619353/1‑43 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTAc                             >  1:1907748/1‑43 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTAc                             >  1:2058120/1‑43 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTAc                             >  1:2384555/1‑43 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTAc                             >  1:2567266/1‑43 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTAc                             >  1:264654/1‑43 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   <  1:631664/69‑1 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   >  1:59691/1‑69 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   <  1:51168/69‑1 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   <  1:3418263/69‑1 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   <  1:694235/69‑1 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   <  1:582385/69‑1 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   >  1:937645/1‑69 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   <  1:574381/69‑1 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   <  1:2688427/69‑1 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   >  1:337963/1‑69 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   >  1:2323508/1‑69 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   <  1:1491478/69‑1 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   >  1:1608338/1‑69 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   <  1:1645268/69‑1 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   >  1:1724805/1‑69 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   <  1:1753777/69‑1 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   <  1:1974448/69‑1 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   >  1:2093444/1‑69 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   >  1:3362761/1‑69 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   >  1:1488098/1‑69 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   <  1:2715934/69‑1 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   <  1:2758134/69‑1 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   >  1:2974527/1‑69 (MQ=255)
                                                          ccGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   <  1:303544/69‑1 (MQ=255)
                                                           cGCAGATTGCCATTAACGTGCAGCAAGGGCTGTCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   <  1:599599/68‑1 (MQ=255)
                                                           cGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggcg  <  1:2393621/69‑1 (MQ=255)
                                                           cGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggcg  <  1:1266866/69‑1 (MQ=255)
                                                           cGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggcg  <  1:525374/69‑1 (MQ=255)
                                                           cGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggcg  <  1:1898773/69‑1 (MQ=255)
                                                           cGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggcg  <  1:1087837/69‑1 (MQ=255)
                                                                aTTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCgcggc   <  1:1435735/63‑1 (MQ=255)
                                                                  |                                                             
GCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTCCGGCCCGTATGGTGATCCGCAGATTGCCATTAACGTGCAGCAAGGGCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCGCGGCG  >  minE/2246476‑2246603

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