Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,669,963 C→T 36.2% N447N (AAC→AAT malP → maltodextrin phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,669,9630CT36.2% 60.5 / 53.5 91N447N (AAC→AATmalPmaltodextrin phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (4/54);  new base T (9/24);  total (13/78)
Fisher's exact test for biased strand distribution p-value = 1.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

GGTGAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAACAGTGCAACCCGGCACTGGCGGCTCTGTTGG  >  minE/2669909‑2670030
                                                      |                                                                   
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACgg                                                          <  1:994992/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACgg                                                          <  1:2362719/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACgg                                                          <  1:2654918/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACgg                                                          <  1:1910151/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACgg                                                          <  1:1909871/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACgg                                                          <  1:1891700/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACgg                                                          <  1:2898169/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACgg                                                          <  1:1720891/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACgg                                                          <  1:3263959/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACgg                                                          <  1:350799/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACgg                                                          <  1:373200/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACgg                                                          <  1:726395/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACgg                                                          <  1:1257060/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACgg                                                          <  1:784114/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACgg                                                          <  1:1221557/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACgg                                                          <  1:1141057/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACgg                                                          <  1:752749/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACgg                                                          <  1:1362787/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACgg                                                          <  1:1204452/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:3131060/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:2184058/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:2167519/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:835481/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:1988558/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:2669164/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:2967355/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:3227574/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:1671419/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:370148/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:1512101/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:593534/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:674264/66‑1 (MQ=255)
ggtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:1271789/66‑1 (MQ=255)
 gtgAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:2074687/65‑1 (MQ=255)
   gAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACGGTATTa                                                     >  1:801642/1‑68 (MQ=255)
   gAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACGGTATTa                                                     >  1:2570120/1‑68 (MQ=255)
   gAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACGGTATTAc                                                    >  1:3306644/1‑69 (MQ=255)
   gAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACGGTATTAc                                                    >  1:82977/1‑69 (MQ=255)
   gAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACGGTATTAc                                                    >  1:560271/1‑69 (MQ=255)
   gAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACGGTATTAc                                                    >  1:1060802/1‑69 (MQ=255)
   gAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACGGTATTAc                                                    >  1:1070110/1‑69 (MQ=255)
   gAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACGGTATTAc                                                    >  1:1837314/1‑69 (MQ=255)
   gAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACGGTATTAc                                                    >  1:1212268/1‑69 (MQ=255)
      aGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACGGTATTAccc                                                  <  1:1530810/68‑1 (MQ=255)
      aGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACGGTATTAccc                                                  <  1:890409/68‑1 (MQ=255)
      aGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACGGTATTAccc                                                  <  1:2588136/68‑1 (MQ=255)
      aGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAATGTCACCAACGGTATTAccc                                                  <  1:2610328/68‑1 (MQ=255)
           tGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGt                                             >  1:952935/1‑68 (MQ=255)
           tGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTc                                            >  1:1940984/1‑69 (MQ=255)
               cccGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:2859528/51‑1 (MQ=255)
               cccGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:361553/51‑1 (MQ=255)
               cccGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:1519438/51‑1 (MQ=255)
               cccGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:1654046/51‑1 (MQ=255)
               cccGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:523620/51‑1 (MQ=255)
               cccGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:1724389/51‑1 (MQ=255)
               cccGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:742155/51‑1 (MQ=255)
               cccGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:2788821/51‑1 (MQ=255)
               cccGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:1135337/51‑1 (MQ=255)
               cccGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:2164944/51‑1 (MQ=255)
               cccGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:2517669/51‑1 (MQ=255)
               cccGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:848745/51‑1 (MQ=255)
               cccGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:2186278/51‑1 (MQ=255)
               cccGGAATATCACCAGCTATGGCAGAACAAATTCCATAACGTCACCAACgg                                                          <  1:1132867/51‑1 (MQ=255)
                     atatCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCa                                   >  1:2178392/1‑68 (MQ=255)
                       atCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAAc                                <  1:2264209/69‑1 (MQ=255)
                       atCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAAc                                <  1:1403871/69‑1 (MQ=255)
                       atCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAAc                                <  1:1428477/69‑1 (MQ=255)
                       atCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAAc                                <  1:650464/69‑1 (MQ=255)
                       atCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAAc                                <  1:380413/69‑1 (MQ=255)
                       atCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAAc                                <  1:1490512/69‑1 (MQ=255)
                       atCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAAc                                <  1:1581784/69‑1 (MQ=255)
                       atCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAAc                                <  1:2296410/69‑1 (MQ=255)
                       atCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAAc                                <  1:2379056/69‑1 (MQ=255)
                       atCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAAc                                <  1:320222/69‑1 (MQ=255)
                        tCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCa                                   <  1:2469502/65‑1 (MQ=255)
                             aGCTATGGCCGATCTAATTCCATAACGTCACCAACTGTATTACCCCACGTCGCTGGATCATAc                                <  1:98030/63‑1 (MQ=38)
                             aGCTATGGCCGAACAAATTCCATAACGTCACCAACgg                                                          <  1:2893139/37‑1 (MQ=255)
                                tATGGCCGAACAAATTCCATAATGTCACCAACGGTATTAccc                                                  <  1:3294823/42‑1 (MQ=255)
                                tATGGCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAACAGTGCAAcc                       >  1:1441688/1‑69 (MQ=255)
                                    gCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAACAGTGCAACCCgg                    <  1:1456142/68‑1 (MQ=255)
                                    gCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAACAGTGCAACCCgg                    <  1:1689453/68‑1 (MQ=255)
                                    gCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAACAGTGCAACCCgg                    <  1:3038885/68‑1 (MQ=255)
                                    gCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAACAGTGCAACCCgg                    <  1:3353649/68‑1 (MQ=255)
                                    gCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAACAGTGCAACCCgg                    <  1:2893197/68‑1 (MQ=255)
                                    gCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAACAGTGCAACCCgg                    <  1:1997646/68‑1 (MQ=255)
                                    gCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAACAGTGCAACCCgg                    <  1:2625883/68‑1 (MQ=255)
                                    gCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAACAGTGCAACCCgg                    <  1:1101230/68‑1 (MQ=255)
                                    gCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAACAGTGCAACCCgg                    <  1:2395014/68‑1 (MQ=255)
                                    gCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAACAGTGCAACCCgg                    <  1:947378/68‑1 (MQ=255)
                                    gCCGAACAAATTCCATAACGTCACCAACGGTATTAACCCACGTCGCTGGATCAAACAGTGCAACCCgg                    <  1:1005888/68‑1 (MQ=255)
                                                     aCGTCACCAACGGTATTACCCCACGTCGCTGGATCAAACAGTGCAACCCGGCACTGGCGGCTCTGTTgg  <  1:1846636/69‑1 (MQ=255)
                                                      cGTCACCAACGGTATTACCCCACGTCGCTGGATCAAACAGTGCAACCCgg                    <  1:768636/50‑1 (MQ=255)
                                                      |                                                                   
GGTGAAAGATCTGTTCCCGGAATATCACCAGCTATGGCCGAACAAATTCCATAACGTCACCAACGGTATTACCCCACGTCGCTGGATCAAACAGTGCAACCCGGCACTGGCGGCTCTGTTGG  >  minE/2669909‑2670030

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: