Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,530,970 G→T 48.9% T216T (ACG→ACT mutM → formamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,530,9700GT48.9% 15.3 / 105.8 96T216T (ACG→ACTmutMformamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base G (23/26);  new base T (15/32);  total (38/58)
Fisher's exact test for biased strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

GTGTGAATTGTTAGCTCGGGTGATTAAAGCGGTGTTGCTGCGTTCGATTGAGCAGGGTGGTACAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAGGAATTGCAGGTTTACGG  >  minE/2530904‑2531035
                                                                  |                                                                 
gtgtGAATTGTTAGCTCGGGTGATTAAAGCGGTGTTGCTGCGTTCGATTGAGCAGGGTGGTACAACGc                                                                  >  1:2083713/1‑68 (MQ=255)
gtgtGAATTGTTAGCTCGGGTGATTAAAGCGGTGTTGCTGCGTTCGATTGAGCAGGGTGGTACAACGc                                                                  >  1:1344232/1‑68 (MQ=255)
gtgtGAATTGTTAGCTCGGGTGATTAAAGCGGTGTTGCTGCGTTCGATTGAGCAGGGTGGTACAACGc                                                                  >  1:1988767/1‑68 (MQ=255)
gtgtGAATTGTTAGCTCGGGTGATTAAAGCGGTGTTGCTGCGTTCGATTGAGCAGGGTGGTACAACGc                                                                  >  1:254001/1‑68 (MQ=255)
gtgtGAATTGTTAGCTCGGGTGATTAAAGCGGTGTTGCTGCGTTCGATTGAGCAGGGTGGTACAACGc                                                                  >  1:3059657/1‑68 (MQ=255)
gtgtGAATTGTTAGCTCGGGTGATTAAAGCGGTGTTGCTGCGTTCGATTGAGCAGGGTGGTACAACGc                                                                  >  1:658900/1‑68 (MQ=255)
gtgtGAATTGTTAGCTCGGGTGATTAAAGCGGTGTTGCTGCGATCGATTGAGCAGGGTGGTACAACGc                                                                  >  1:1583907/1‑68 (MQ=255)
     aaTTGTTAGCTCGGGTGATTAAAGCGGTGTTGCTGCGTTCGATTGAGCAGGGTGGTACAACGCTGAAAg                                                            <  1:3377444/69‑1 (MQ=255)
     aaTTGTTAGCTCGGGTGATTAAAGCGGTGTTGCTGCGTTCGATTGAGCAGGGTGGTACAACGCTGAAAg                                                            <  1:756111/69‑1 (MQ=255)
                         aaaGCGGTGTTGCTGCGTTCGATTGAGCAGGGTGGTACAACGCTGAAAGATTTTCTGCAAAGTGa                                            <  1:2423779/65‑1 (MQ=255)
                         aaaGCGGTGTTGCTGCGTTCGATTGAGCAGGGTGGTACAACGCTGAAAGATTTTCTGCAAAGTGa                                            <  1:231059/65‑1 (MQ=255)
                         aaaGCGGTGTTGCTGCGTTCGATTGAGCAGGGTGGTACAACGCTGAAAGATTTTCTGCAAAGTGa                                            <  1:2105587/65‑1 (MQ=255)
                         aaaGCGGTGTTGCTGCGTTCGATTGAGCAGGGTGGTACAACGCTGAAAGATTTTCTGCAAAGTGa                                            <  1:1887259/65‑1 (MQ=255)
                         aaaGCGGTGTTGCTGCGTTCGATTGAGCAGGGTGGTACAACGCTGAAAGATTTTCTGCAAAGTGa                                            <  1:1576936/65‑1 (MQ=255)
                                 gtTGCTGCGTTCGATTGAGCAGGGTGGTACAACGCTGAAAGATTTTCt                                                     <  1:1617788/48‑1 (MQ=255)
                                 gtTGCTGCGTTCGATTGAGCAGGGTGGTACAACGCTGAAAGATTTTCt                                                     <  1:1316944/48‑1 (MQ=255)
                                        cGTTCGATTGAGCAGGGTGGTACAACGCTGAAAGa                                                           >  1:999095/1‑35 (MQ=255)
                                        cGTTCGATTGAGCAGGGTGGTACAACGCTGAAAGa                                                           >  1:1832174/1‑35 (MQ=255)
                                        cGTTCGATTGAGCAGGGTGGTACAACGCTGAAAGa                                                           >  1:1547002/1‑35 (MQ=255)
                                        cGTTCGATTGAGCAGGGTGGTACAACGCTGAAAGa                                                           >  1:1311891/1‑35 (MQ=255)
                                        cGTTCGATTGAGCAGGGTGGTACAACGCTGAAAGa                                                           >  1:1343879/1‑35 (MQ=255)
                                        cGTTCGATTGAGCAGGGTGGTACAACGCTGAAAGa                                                           >  1:1272727/1‑35 (MQ=255)
                                           tCGATTGAGCAGGGTGGTACAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCt                      <  1:3322442/69‑1 (MQ=255)
                                               ttGAGCAGGGTGGTACAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAgg                  <  1:1244324/69‑1 (MQ=255)
                                               ttGAGCAGGGTGGTACAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAgg                  <  1:2723068/69‑1 (MQ=255)
                                               ttGAGCAGGGTGGTACAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAgg                  <  1:2719131/69‑1 (MQ=255)
                                               ttGAGCAGGGTGGTACAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAgg                  <  1:2573358/69‑1 (MQ=255)
                                               ttGAGCAGGGTGGTACAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAgg                  <  1:3106473/69‑1 (MQ=255)
                                               ttGAGCAGGGTGGTACAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAgg                  <  1:2628568/69‑1 (MQ=255)
                                               ttGAGCAGGGTGGTACAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAgg                  <  1:1711104/69‑1 (MQ=255)
                                                tGAGCAGGGTGGTACAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCTGGCTATTTCGCTCAGGa                 <  1:2428402/69‑1 (MQ=255)
                                                tGAGCAGGGTGGTACAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAGGa                 <  1:2908289/69‑1 (MQ=255)
                                                tGAGCAGGGTGGTACAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAGGa                 <  1:2608540/69‑1 (MQ=255)
                                                tGAGCAGGGTGGTACAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAGGa                 <  1:2678461/69‑1 (MQ=255)
                                                tGAGCAGGGTGGTACAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAGGa                 <  1:2903853/69‑1 (MQ=255)
                                                     aGGGTGGTACAACTCTGAAAGATTTTCTGCAAAGTGATgg                                         >  1:61467/1‑40 (MQ=255)
                                                     aGGGTGGTACAACTCTGAAAGATTTTCTGCAAAGTGATg                                          >  1:757066/1‑39 (MQ=255)
                                                     aGGGTGGTACAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCGCAGGAAt               >  1:2513195/1‑66 (MQ=255)
                                                     aGGGTGGTACAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCGCAGGAAt               >  1:2600378/1‑66 (MQ=255)
                                                     aGGGTGGTACAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCGCAGGAAt               >  1:3232435/1‑66 (MQ=255)
                                                     aGGGTGGTACAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCGCAGGAAt               >  1:1080721/1‑66 (MQ=255)
                                                     aGGGTGGTACAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCGCAGGAAt               >  1:2864991/1‑66 (MQ=255)
                                                     aGGGTGGTACAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCGCAGGAAt               >  1:176845/1‑66 (MQ=255)
                                                     aGGGTGGTACAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCGCAGGAAt               >  1:1496096/1‑66 (MQ=255)
                                                     aGGGTGGTACAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCGCAGGAAt               >  1:1520892/1‑66 (MQ=255)
                                                     aGGGTGGTACAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCGCAGGAAt               >  1:1678569/1‑66 (MQ=255)
                                                     aGGGTGGTACAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCGCAGGAAt               >  1:640029/1‑66 (MQ=255)
                                                     aGGGTGGTACAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCGCAGGAAt               >  1:1771942/1‑66 (MQ=255)
                                                     aGGGTGGTACAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCGCAGGAAt               >  1:556863/1‑66 (MQ=255)
                                                     aGGGTGGTACAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCGCAGGAAt               >  1:1253285/1‑66 (MQ=255)
                                                         tggtACAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAGGAATTGCAgg         >  1:2637180/1‑68 (MQ=255)
                                                           gtACAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCg                                  >  1:1208607/1‑41 (MQ=255)
                                                           gtACAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAGGAATTGCAGGttt      >  1:1128159/1‑69 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:2948372/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:289478/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:382763/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:956120/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:954002/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:849975/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:554481/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:1100859/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:1058302/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:1028915/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:3099339/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:3053457/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:585819/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:27013/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:1310792/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:15833/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:1732763/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:1738055/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:1807408/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:1253475/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:1883352/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:189307/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:1898386/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:2147867/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:1224161/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:1193331/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:2335497/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:2337897/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:2566433/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:277150/39‑1 (MQ=255)
                                                              cAACTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:2632844/39‑1 (MQ=255)
                                                              cAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAGGTATTGCAGGTTTACg   <  1:1751000/69‑1 (MQ=255)
                                                              cAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAGGAATTGCAGGTTTACg   <  1:647989/69‑1 (MQ=255)
                                                              cAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAGGAATTGCAGGTTTACg   <  1:2675243/69‑1 (MQ=255)
                                                               aaCTCTGAAAGATTTTCTGCAAAGTGATGGTAAACCgg                                 <  1:1258430/38‑1 (MQ=255)
                                                                aCGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAGGAATTGCAGGtt       >  1:1557815/1‑63 (MQ=255)
                                                                aCGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAGGAATTGCAGGtt       >  1:957235/1‑63 (MQ=255)
                                                                aCGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAGGAATTGCAGGtt       >  1:94167/1‑63 (MQ=255)
                                                                aCGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAGGAATTGCAGGtt       >  1:902962/1‑63 (MQ=255)
                                                                aCGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAGGAATTGCAGGtt       >  1:1036469/1‑63 (MQ=255)
                                                                aCGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAGGAATTGCAGGtt       >  1:549626/1‑63 (MQ=255)
                                                                aCGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAGGAATTGCAGGTTTAcag  >  1:715113/1‑66 (MQ=255)
                                                                aCGCTGAAAGATTATCTGCAAAGTGATGGTAAACCGGGCTATTTCGCt                      <  1:1155309/48‑1 (MQ=255)
                                                                  |                                                                 
GTGTGAATTGTTAGCTCGGGTGATTAAAGCGGTGTTGCTGCGTTCGATTGAGCAGGGTGGTACAACGCTGAAAGATTTTCTGCAAAGTGATGGTAAACCGGGCTATTTCGCTCAGGAATTGCAGGTTTACGG  >  minE/2530904‑2531035

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: