Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,447,246 A→G 42.7% pseudogene (810/1647 nt) ilvG ← Mutation in protein compared to W3110: p.Glu327_Ilefs*223 Restoration of IlvG found in MG1655 (548 amino acids), fusion of N‑ter and C‑ter fragments ECK3760:JW3740+JW3741:b4488; large subunit of acetolactate synthase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,447,2460AG42.7% 23.0 / 68.2 89pseudogene (810/1647 nt)ilvGMutation in protein compared to W3110: p.Glu327_Ilefs*223 Restoration of IlvG found in MG1655 (548 amino acids), fusion of N‑ter and C‑ter fragments ECK3760:JW3740+JW3741:b4488; large subunit of acetolactate synthase II
Reads supporting (aligned to +/- strand):  ref base A (38/13);  new base G (25/13);  total (63/26)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.12e-01

TATGGATAACACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTCGCCGCTTTGGTGCCATGCATTCCC  >  minE/2447179‑2447312
                                                                   |                                                                  
tatGGATAACACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCa                                                                    <  1:1097760/68‑1 (MQ=255)
tatGGATAACACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCa                                                                    <  1:1475799/68‑1 (MQ=255)
  tGGATAACACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCCCCCGGTCATCAAAACGTGCAcc                                                                  >  1:1839755/1‑68 (MQ=255)
  tGGATAACACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCAccc                                                                 >  1:2416502/1‑69 (MQ=255)
  tGGATAACACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCAccc                                                                 >  1:1151336/1‑69 (MQ=255)
  tGGATAACACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCAccc                                                                 >  1:1462919/1‑69 (MQ=255)
    gATAACACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCACCCa                                                                <  1:3307126/68‑1 (MQ=255)
    gATAACACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCACCCa                                                                <  1:2505710/68‑1 (MQ=255)
       aaCACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCCCCCGGTCATCAAAACGTGCACCCACGGc                                                            >  1:2215951/1‑69 (MQ=255)
       aaCACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCACCCACGGc                                                            >  1:1030015/1‑69 (MQ=255)
       aaCACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCACCCACGGc                                                            >  1:511476/1‑69 (MQ=255)
       aaCACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCACCCACGGc                                                            >  1:3417582/1‑69 (MQ=255)
       aaCACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCACCCACGGc                                                            >  1:2720970/1‑69 (MQ=255)
       aaCACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCACCCACGGc                                                            >  1:2615687/1‑69 (MQ=255)
       aaCACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCACCCACGGc                                                            >  1:1632383/1‑69 (MQ=255)
       aaCACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCACCCACGGc                                                            >  1:1475579/1‑69 (MQ=255)
       aaCACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCACCCACGGc                                                            >  1:1063832/1‑69 (MQ=255)
       aaCACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCACCCACGGc                                                            >  1:1048079/1‑69 (MQ=255)
       aaCACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCACCCAAGGc                                                            >  1:2190983/1‑69 (MQ=255)
       aaCACTGGCGTGTGGTGCGAAGGTGTTAAGTTTGACGGTCACCCGGTCATCAAAACGTGCAccc                                                                 >  1:1903081/1‑64 (MQ=39)
                        cGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCACCCACGGCGATCAGCAAGTCGCACt                                           >  1:3180236/1‑69 (MQ=255)
                        cGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCACCCACGGCGATCAGCAAGTCGCACt                                           >  1:1301261/1‑69 (MQ=255)
                             gtgtTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCACCCa                                                                >  1:1468409/1‑43 (MQ=255)
                             gtgtTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGc                                      >  1:509229/1‑69 (MQ=255)
                             gtgtTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGc                                      <  1:1141188/69‑1 (MQ=255)
                             gtgtTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGc                                      >  1:841321/1‑69 (MQ=255)
                                  cAGTTTGCCAGTCACCCGGTCATCAAAACGTGCGCCCGCGGCGATCAGTAGATCGCACTCCTGCActg                                  >  1:3284917/1‑66 (MQ=25)
                                  cAGTTTGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCActg                                  >  1:3297696/1‑66 (MQ=37)
                                  cAGTTTGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCActg                                  >  1:628462/1‑66 (MQ=37)
                                  cAGTTTGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCActg                                  >  1:3278031/1‑66 (MQ=37)
                                  cAGTTTGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCActg                                  >  1:759418/1‑66 (MQ=37)
                                  cAGTTTGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCActg                                  >  1:3092447/1‑66 (MQ=37)
                                  cAGTTTGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCActg                                  >  1:2988825/1‑66 (MQ=37)
                                  cAGTTTGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCActg                                  >  1:1452978/1‑66 (MQ=37)
                                  cAGTTTGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCActg                                  >  1:2726628/1‑66 (MQ=37)
                                  cAGTTTGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCActg                                  >  1:283948/1‑66 (MQ=37)
                                  cAGTTTGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCActg                                  >  1:2371885/1‑66 (MQ=37)
                                  cAGTTTGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCActg                                  >  1:2354101/1‑66 (MQ=37)
                                  cAGTTTGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCActg                                  >  1:1604107/1‑66 (MQ=37)
                                  cAGTTTGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCActg                                  >  1:1731173/1‑66 (MQ=37)
                                  cAGTTTGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCActg                                  >  1:1101481/1‑66 (MQ=37)
                                  cAGTTTGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCActg                                  >  1:1917465/1‑66 (MQ=37)
                                  cAGTTTGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCActg                                  >  1:2153053/1‑66 (MQ=37)
                                  cAGTTTGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCActg                                  >  1:2077265/1‑66 (MQ=37)
                                  cAGTTTGCCAGTAACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCActg                                  >  1:3426814/1‑66 (MQ=255)
                                       tGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTTCACtgcg                                <  1:2315502/65‑1 (MQ=25)
                                       tGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCACtgcg                                <  1:2331725/65‑1 (MQ=37)
                                       tGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCACtgcg                                <  1:1047772/65‑1 (MQ=37)
                                       tGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCACtgcg                                <  1:179519/65‑1 (MQ=37)
                                       tGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCACtgcg                                <  1:3225984/65‑1 (MQ=37)
                                       tGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCACtgcg                                <  1:479703/65‑1 (MQ=37)
                                       tGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCACtgcg                                <  1:3074879/65‑1 (MQ=37)
                                       tGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCACtgcg                                <  1:1528856/65‑1 (MQ=37)
                                       tGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCACtgcg                                <  1:101782/65‑1 (MQ=37)
                                       tGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCACtgcg                                <  1:2379793/65‑1 (MQ=37)
                                       tGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCACtgcg                                <  1:1051650/65‑1 (MQ=37)
                                       tGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCACtgcg                                <  1:531515/65‑1 (MQ=37)
                                       tGCCAGTCACCCGGTCATCAAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCACtgcg                                <  1:520747/65‑1 (MQ=37)
                                              cacaCGGTCATCAAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTcgccgc                     >  1:764155/1‑69 (MQ=255)
                                              cACCCGGTCATCAAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTcgccgc                     >  1:1059379/1‑69 (MQ=255)
                                              cACCCGGTCATCAAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTcgccgc                     >  1:2381462/1‑69 (MQ=255)
                                              cACCCGGTCATCAAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTcgccgc                     >  1:1363081/1‑69 (MQ=255)
                                              cACCCGGTCATCAAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTcgccgc                     >  1:328175/1‑69 (MQ=255)
                                              cACCCGGTCATCAAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTcgccgc                     >  1:539459/1‑69 (MQ=255)
                                              cACCCGGTCATCAAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTcgccgc                     >  1:2342289/1‑69 (MQ=255)
                                              cACCCGGTCATCAAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTcgccgc                     >  1:427607/1‑69 (MQ=255)
                                              cACCCGGTCATCAAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTcgccgc                     >  1:1983789/1‑69 (MQ=255)
                                              cACCCGGTCATCAAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTcgccgc                     >  1:2707859/1‑69 (MQ=255)
                                              cACCCGGTCATCAAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTcgccgc                     >  1:2785814/1‑69 (MQ=255)
                                              cACCCGGTCATCAAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTcgccgc                     >  1:2886710/1‑69 (MQ=255)
                                                      catcaAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGTGAAGTTCGCCGCTTTGGTGc             <  1:351597/69‑1 (MQ=255)
                                                      catcaAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTCGCCGCTTTGGTGc             <  1:3254901/69‑1 (MQ=255)
                                                      catcaAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTCGCCGCTTTGGTGc             <  1:333618/69‑1 (MQ=255)
                                                      catcaAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTCGCCGCTTTGGTGc             <  1:1968187/69‑1 (MQ=255)
                                                      catcaAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTCGCCGCTTTGGTGc             <  1:3174420/69‑1 (MQ=255)
                                                      catcaAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTCGCCGCTTTGGTGc             <  1:1083168/69‑1 (MQ=255)
                                                      catcaAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTCGCCGCTTTGGTGc             <  1:3072130/69‑1 (MQ=255)
                                                         caAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCACTGCGAAGTTCGCCGCTTTGGTGccg           >  1:872374/1‑67 (MQ=255)
                                                         caAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCACTGCGAAGTTCGCCGCTTTGGTGccg           >  1:2844618/1‑67 (MQ=255)
                                                         caAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCACTGCGAAGTTCGCCGCTTTGGTGccg           >  1:1483585/1‑67 (MQ=255)
                                                         caAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCACTGCGAAGTTCGCCGCTTTGGTGccg           >  1:2693802/1‑67 (MQ=255)
                                                         caAAACGTGCGCCCACGGCGATCAGTAGATCGCACTCCTGCACTGCGAAGTTCGCCGCTTTGGTGccg           >  1:573745/1‑67 (MQ=255)
                                                         caAAACGTGCGCCCAAGGCGATCAGTAGATCGCACTCCTGCACTGCGAAGTTCGCCGCTTTGGTGccg           >  1:887406/1‑67 (MQ=255)
                                                         caAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGtt                           >  1:1360454/1‑52 (MQ=255)
                                                         caAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTCGCCGCTTTGGTGCCa           >  1:2553596/1‑68 (MQ=255)
                                                         caAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTCGCCGCTTTGGTGCCa           >  1:2109375/1‑68 (MQ=255)
                                                         caAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTCGCCGCTTTGGTGCCa           >  1:478196/1‑68 (MQ=255)
                                                         caAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTCGCCGCTTTGGTGCCa           >  1:3241463/1‑68 (MQ=255)
                                                          aaaaCGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTcgccgc                     <  1:695869/57‑1 (MQ=255)
                                                          aaaaCGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTCGCCGCTTTGGTGCCATg         <  1:482346/69‑1 (MQ=255)
                                                                 gCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTCGCCGCTTTGGTGCCATGCATTccc  <  1:24943/69‑1 (MQ=255)
                                                                   |                                                                  
TATGGATAACACTGGCGTGTGGTGCGAAGGTGTTCAGTTTGCCGGTCACCCGGTCATCAAAACGTGCACCCACGGCGATCAGCAAGTCGCACTCCTGCACCGCGAAGTTCGCCGCTTTGGTGCCATGCATTCCC  >  minE/2447179‑2447312

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: