Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,447,498 A→G 57.6% pseudogene (558/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,4980AG57.6% 47.1 / 95.4 128pseudogene (558/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 (23/31);  new base G (16/57);  total (39/89)
Fisher's exact test for biased strand distribution p-value = 1.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

CCCACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCGGCATGTGGGAAAGTCACTTCGTTTTCAACGGTGGTGAAC  >  minE/2447433‑2447543
                                                                 |                                             
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAAc                                             >  1:3306578/1‑68 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACt                                            >  1:257941/1‑69 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACt                                            >  1:2809467/1‑69 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACt                                            >  1:1739046/1‑69 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACt                                            >  1:2456711/1‑69 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACt                                            >  1:3146068/1‑69 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACt                                            >  1:2274010/1‑69 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACt                                            >  1:1890697/1‑69 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTAAACt                                            >  1:1294227/1‑69 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGAc                                             >  1:3144788/1‑68 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            <  1:1621807/69‑1 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            <  1:1801145/69‑1 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            >  1:1648969/1‑69 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            >  1:1825465/1‑69 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            >  1:1873453/1‑69 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            <  1:1896812/69‑1 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            >  1:2037292/1‑69 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            <  1:2070871/69‑1 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            <  1:3240149/69‑1 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            >  1:2241332/1‑69 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            >  1:2242942/1‑69 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            <  1:2246697/69‑1 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            <  1:2958353/69‑1 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            <  1:278775/69‑1 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            <  1:257765/69‑1 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            >  1:1296279/1‑69 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            <  1:983304/69‑1 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            >  1:1023588/1‑69 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            <  1:96522/69‑1 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            <  1:840469/69‑1 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            <  1:1113485/69‑1 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            <  1:1117103/69‑1 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            <  1:783611/69‑1 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            <  1:1223350/69‑1 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            <  1:1289143/69‑1 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            >  1:40904/1‑69 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            <  1:1424474/69‑1 (MQ=255)
cccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            >  1:581177/1‑69 (MQ=255)
 ccACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACt                                            <  1:2415925/68‑1 (MQ=255)
    cGCCACCGCCAACTTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCgg                                        <  1:1987032/69‑1 (MQ=255)
    cGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCgg                                        <  1:3307924/69‑1 (MQ=255)
    cGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCgg                                        <  1:1514149/69‑1 (MQ=255)
    cGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCgg                                        <  1:410883/69‑1 (MQ=255)
    cGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCgg                                        <  1:1799879/69‑1 (MQ=255)
    cGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCgg                                        <  1:3068942/69‑1 (MQ=255)
    cGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCgg                                        <  1:31199/69‑1 (MQ=255)
    cGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCgg                                        <  1:1844807/69‑1 (MQ=255)
    cGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCgg                                        <  1:2318567/69‑1 (MQ=255)
    cGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCgg                                        <  1:1198574/69‑1 (MQ=255)
    cGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCgg                                        <  1:3413508/69‑1 (MQ=255)
    cGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCgg                                        <  1:52463/69‑1 (MQ=255)
    cGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCgg                                        <  1:2125384/69‑1 (MQ=255)
    cGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCgg                                        <  1:1353645/69‑1 (MQ=255)
                 gTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGt                           <  1:73947/69‑1 (MQ=255)
                 gTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGt                           <  1:729111/69‑1 (MQ=255)
                 gTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGt                           <  1:412732/69‑1 (MQ=255)
                 gTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGt                           <  1:1009982/69‑1 (MQ=255)
                 gTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGt                           <  1:366224/69‑1 (MQ=255)
                 gTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGt                           <  1:721634/69‑1 (MQ=255)
                 gTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGt                           <  1:631905/69‑1 (MQ=255)
                 gTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGt                           <  1:263371/69‑1 (MQ=255)
                 gTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGt                           <  1:2496031/69‑1 (MQ=255)
                 gTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGt                           <  1:1179786/69‑1 (MQ=255)
                 gTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGt                           <  1:1479242/69‑1 (MQ=255)
                 gTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGt                           <  1:2470746/69‑1 (MQ=255)
                 gTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGt                           <  1:1959682/69‑1 (MQ=255)
                 gTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGt                           <  1:2305186/69‑1 (MQ=255)
                    cAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCGGCATGTGGGAAAGTCa                         >  1:1590409/1‑68 (MQ=255)
                    cAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCGGCATGTGGGAAAGTCAc                        >  1:1005111/1‑69 (MQ=255)
                    cAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCGGCATGTGGGAAAGTCAc                        >  1:3174302/1‑69 (MQ=255)
                    cAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCGGCATGTGGGAAAGTCAc                        >  1:2355151/1‑69 (MQ=255)
                           ggTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCGGCAtg                                    >  1:505588/1‑50 (MQ=255)
                           ggTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCGGCAtg                                    >  1:1379150/1‑50 (MQ=255)
                               tttGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCGGCATGTGGGAAAGTCACTTCGTTTTCaa             >  1:1046948/1‑69 (MQ=255)
                               tttGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCGGCATGTGGGAAAGTCACTTCGTTTTCaa             >  1:859587/1‑69 (MQ=255)
                               tttGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGTCACTTCGTTTTCaa             <  1:1151428/69‑1 (MQ=255)
                               tttGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGTCACTTCGTTTTCaa             <  1:582538/69‑1 (MQ=255)
                               tttGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGTCACTTCGTTTTCaa             <  1:804708/69‑1 (MQ=255)
                               tttGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGTCACTTCGTTTTCaa             <  1:2323288/69‑1 (MQ=255)
                               tttGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGTCACTTCGTTTTCaa             <  1:514234/69‑1 (MQ=255)
                               tttGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGTCACTTCGTTTTCaa             <  1:1513518/69‑1 (MQ=255)
                               tttGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGTCACTTCGTTTTCaa             <  1:3012504/69‑1 (MQ=255)
                               tttGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGTCACTTCGTTTTCaa             <  1:1867540/69‑1 (MQ=255)
                               tttGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGTCACTTCGTTTTCaa             <  1:1940074/69‑1 (MQ=255)
                               tttGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGTCACTTCGTTTTCaa             <  1:960730/69‑1 (MQ=255)
                               tttGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGTCACTTCGTTTTCaa             <  1:3295315/69‑1 (MQ=255)
                               tttGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGTCACTTCGTTTTCaa             <  1:2083853/69‑1 (MQ=255)
                               tttGCGCTTTTGCCAGCATCTGACGCGCTTGCTCGACTTCGGCATGTGAGAAAGTCACTTCGTTTTCaa             <  1:137758/69‑1 (MQ=255)
                                    gcTTTTGCCAGCATCTGGCGCGCTTGCTCTACTTCGGCATGTGGGAAAGTCACTTCGTTTTCAACggt         <  1:1197628/68‑1 (MQ=255)
                                    gcTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCGGCATGTGGGAAAGTCACTTCGTTTTCAACggt         <  1:1103331/68‑1 (MQ=255)
                                    gcTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCGGCATGTGGGAAAGTCACTTCGTTTTCAACggt         <  1:2563811/68‑1 (MQ=255)
                                    gcTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCGGCATGTGGGAAAGTCACTTCGTTTTCAACggt         <  1:1451266/68‑1 (MQ=255)
                                    gcTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCGGCATGTGGGAAAGTCACTTCGTTTTCAACggt         <  1:6118/68‑1 (MQ=255)
                                    gcTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCGGCATGTGGGAAAGTCACTTCGTTTTCAACggt         <  1:1443027/68‑1 (MQ=255)
                                    gcTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCGGCATGTGGGAAAGTCACTTCGTTTTCAACggt         <  1:1759134/68‑1 (MQ=255)
                                    gcTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCGGCATGTGGGAAAGTCACTTCGTTTTCAACggt         <  1:3019533/68‑1 (MQ=255)
                                    gcTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCGGCATGTGGGAAAGTCACTTCGTTTTCAACggt         <  1:319570/68‑1 (MQ=255)
                                           ccAGCATCTGGCGCGCTTGCTCAACTTCGGCATGTGGGAAAGTCACTTCGTTTTCAACGGTGGTGAAc  <  1:2846924/68‑1 (MQ=255)
                                           ccAGCATCTGGCGCGCTTGCTCAACTTCGGCATGTGGGAAAGTCACTTCGTTTTCAACGGTGGTGAAc  <  1:2909139/68‑1 (MQ=255)
                                           ccAGCATCTGGCGCGCTTGCTCAACTTCGGCATGTGGGAAAGTCACTTCGTTTTCAACGGTGGTGAAc  <  1:2351657/68‑1 (MQ=255)
                                                                 |                                             
CCCACGCCACCGCCAACGTACAGCATCGGTTTTTGCGCTTTTGCCAGCATCTGGCGCGCTTGCTCAACTTCGGCATGTGGGAAAGTCACTTCGTTTTCAACGGTGGTGAAC  >  minE/2447433‑2447543

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: