Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,447,032 C→A 27.7% pseudogene (1024/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
*Exported2,447,0320CA27.7% 148.0 / 42.3 72pseudogene (1024/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 C (21/31);  new base A (12/8);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

ACAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTAACGGCTGCTGTAATGCT  >  Exported/2446972‑2447099
                                                            |                                                                   
aCAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTtcat                                                             <  1:5173303/69‑1 (MQ=255)
aCAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTtcat                                                             <  1:356536/69‑1 (MQ=255)
aCAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTtcat                                                             <  1:4348284/69‑1 (MQ=255)
aCAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTtcat                                                             <  1:1376178/69‑1 (MQ=255)
aCAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTtcat                                                             <  1:957482/69‑1 (MQ=255)
       ttttAACAACAACGGCGCGTAGATAGCTTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGc                                                     <  1:4385217/70‑1 (MQ=255)
       ttttAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGc                                                     <  1:1128517/70‑1 (MQ=255)
       ttttAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGc                                                     <  1:3908789/70‑1 (MQ=255)
       ttttAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCt                                                    >  1:4060920/1‑71 (MQ=255)
        tttAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTtcatca                                                           <  1:4691927/63‑1 (MQ=255)
        tttAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTtcatca                                                           <  1:4164695/63‑1 (MQ=255)
              atcaacGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCa                                             <  1:4025463/69‑1 (MQ=255)
              aacaacGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCa                                             <  1:47395/71‑1 (MQ=255)
              aacaacGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCa                                             <  1:4019423/71‑1 (MQ=255)
              aacaacGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCa                                             <  1:2457290/71‑1 (MQ=255)
                  acGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAgtgt                                         <  1:2806725/71‑1 (MQ=255)
                  acGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAgtgt                                         <  1:2917133/71‑1 (MQ=255)
                   cGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTa                                        <  1:934855/71‑1 (MQ=255)
                   cGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTa                                        <  1:913713/71‑1 (MQ=255)
                   cGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTa                                        <  1:3987740/71‑1 (MQ=255)
                        cgTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGc                                   <  1:3102548/71‑1 (MQ=255)
                         gTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGcc                                  <  1:3046489/71‑1 (MQ=255)
                         gTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGc                                   <  1:2282211/70‑1 (MQ=255)
                            gATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGt                               >  1:2723095/1‑71 (MQ=255)
                            gATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGt                               >  1:1227460/1‑71 (MQ=255)
                               aGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCt                                                    >  1:2472992/1‑47 (MQ=255)
                               aGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCAt                            >  1:3202408/1‑71 (MQ=255)
                               aGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCAt                            >  1:1156933/1‑71 (MQ=255)
                               aGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCAt                            >  1:3341492/1‑71 (MQ=255)
                               aGCGTCACCAGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCAt                            >  1:1204739/1‑71 (MQ=255)
                                gCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCAtt                           <  1:3775739/71‑1 (MQ=255)
                                gCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCAtt                           <  1:990658/71‑1 (MQ=255)
                                gCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCAtt                           <  1:3203704/71‑1 (MQ=255)
                                gCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCAtt                           <  1:3504149/71‑1 (MQ=255)
                                gCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCAtt                           <  1:4791185/71‑1 (MQ=255)
                                   tCACCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTgcgcgc                                              <  1:5040860/49‑1 (MQ=255)
                                   tCACCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTgcgcgc                                              <  1:4477580/49‑1 (MQ=255)
                                   tCACCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTgcgcgc                                              <  1:3824236/49‑1 (MQ=255)
                                   tCACCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTgcgcgc                                              <  1:3301418/49‑1 (MQ=255)
                                   tCACCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTgcgcgc                                              <  1:1565021/49‑1 (MQ=255)
                                   tCACCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTgcgcgc                                              <  1:2600424/49‑1 (MQ=255)
                                   tCACCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTgcgcgc                                              <  1:581126/49‑1 (MQ=255)
                                    cACCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTgcgcgc                                              <  1:291868/48‑1 (MQ=255)
                                     aCCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGata                       <  1:3454894/70‑1 (MQ=255)
                                     aCCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTGCGCGCAg                                            >  1:2103439/1‑49 (MQ=255)
                                     aCCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTGCGCGCAg                                            >  1:1210790/1‑49 (MQ=255)
                                     aCCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTGCGCGCAg                                            >  1:98786/1‑49 (MQ=255)
                                     aCCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTGCGCGCAg                                            >  1:4290368/1‑49 (MQ=255)
                                     aCCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTGCGCGCAg                                            >  1:177476/1‑49 (MQ=255)
                                     aCCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTGCGCGCAg                                            >  1:1797936/1‑49 (MQ=255)
                                     aCCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTGCGCGCAg                                            >  1:1850842/1‑49 (MQ=255)
                                     aCCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTGCGCGCAg                                            >  1:1964121/1‑49 (MQ=255)
                                     aCCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTGCGCGCAg                                            >  1:4977332/1‑49 (MQ=255)
                                     aCCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTGCGCGCAg                                            >  1:266362/1‑49 (MQ=255)
                                     aCCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTGCGCGCAg                                            >  1:3561136/1‑49 (MQ=255)
                                     aCCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTGCGCGCAg                                            >  1:3381303/1‑49 (MQ=255)
                                       cGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCTCAGTGTAGCTGCCAGTCATTGATAttt                    <  1:467162/71‑1 (MQ=255)
                                        gggATGGTCGTAACGCCAGGCATGTTCATCACGCACCTGCGCGCAGTGTAGCTGCCAGTCAt                            >  1:4628820/1‑62 (MQ=255)
                                         ggATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTaa                  >  1:1132825/1‑71 (MQ=255)
                                         ggATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTa                   >  1:3410924/1‑70 (MQ=255)
                                              gTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTAACGgct             >  1:2952270/1‑71 (MQ=255)
                                              gTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTAACGgct             >  1:224739/1‑71 (MQ=255)
                                                 gTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTAACGgctgct          >  1:4193882/1‑71 (MQ=255)
                                                 gTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTAACGgctgct          >  1:1386039/1‑71 (MQ=255)
                                                 gTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTAACGgctgct          >  1:1201697/1‑71 (MQ=255)
                                                  tAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTAACGGctgctg         <  1:998092/71‑1 (MQ=255)
                                                    aCGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTAACGGCTGCTGt        <  1:1737256/70‑1 (MQ=255)
                                                        cAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTAACGGCTGCTGTAATGc   >  1:4038312/1‑71 (MQ=255)
                                                        cAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTAACGGCTGCTGTAATGc   >  1:1357983/1‑71 (MQ=255)
                                                         aGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTAACGGCTGCTGTAATGCt  >  1:1344415/1‑71 (MQ=255)
                                                           gCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGat                        >  1:3450071/1‑47 (MQ=255)
                                                           gCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGat                        >  1:2625337/1‑47 (MQ=255)
                                                            |                                                                   
ACAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTAACGGCTGCTGTAATGCT  >  Exported/2446972‑2447099

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: