Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,914,374 G→T 73.1% I366I (ATC→ATA ptsP ← fused PEP‑protein phosphotransferase (enzyme I) of PTS system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,914,3740GT73.1% 114.2 / 53.8 93I366I (ATC→ATAptsPfused PEP‑protein phosphotransferase (enzyme I) of PTS system
Reads supporting (aligned to +/- strand):  ref base G (14/11);  new base T (39/29);  total (53/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

GCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACA  >  minE/1914308‑1914437
                                                                  |                                                               
gCCATCTACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGAt                                                               >  1:440339/1‑69 (MQ=255)
gCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGAt                                                               >  1:2855097/1‑69 (MQ=255)
gCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGAt                                                               >  1:3391775/1‑69 (MQ=255)
gCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGAt                                                               >  1:1482446/1‑69 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATcccc                                                           >  1:1016943/3‑68 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATcccc                                                           >  1:2359381/3‑68 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATcccc                                                           >  1:3259727/3‑68 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATcccc                                                           >  1:3315943/3‑68 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATcccc                                                           >  1:1799261/3‑68 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCa                                                          >  1:1405008/3‑69 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCa                                                          >  1:2927404/3‑69 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCa                                                          >  1:2315435/3‑69 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCa                                                          >  1:2412979/3‑69 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCa                                                          >  1:2221339/3‑69 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCa                                                          >  1:517612/3‑69 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCa                                                          >  1:3328164/3‑69 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCa                                                          >  1:3370822/3‑69 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCa                                                          >  1:2438543/3‑69 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCa                                                          >  1:1961518/3‑69 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCa                                                          >  1:1916396/3‑69 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCa                                                          >  1:1081603/3‑69 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCa                                                          >  1:1892668/3‑69 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCa                                                          >  1:343435/3‑69 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCa                                                          >  1:1478724/3‑69 (MQ=255)
     cgaCGATCAGCGTCCGACGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCa                                                          >  1:2512023/3‑69 (MQ=255)
                gTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTAcc                                               <  1:2022238/69‑1 (MQ=255)
                gTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTAcc                                               <  1:2245738/69‑1 (MQ=255)
                gTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTAcc                                               <  1:2663795/69‑1 (MQ=255)
                 tCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTAcc                                               >  1:1526735/1‑68 (MQ=255)
                 tCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTAcc                                               >  1:216206/1‑68 (MQ=255)
                 tCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTAcc                                               >  1:1226538/1‑68 (MQ=255)
                   cTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATg                                            <  1:2631243/69‑1 (MQ=255)
                   cTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATg                                            <  1:1807866/69‑1 (MQ=255)
                   cTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATg                                            <  1:1217682/69‑1 (MQ=255)
                   cTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATg                                            <  1:487107/69‑1 (MQ=255)
                   cTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATg                                            <  1:1721575/69‑1 (MQ=255)
                    tACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATg                                            <  1:3365562/68‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGGGCACGTACCATGAt                                          <  1:2509936/69‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:2698171/69‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:240377/69‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:2809585/69‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:2760007/69‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:254747/69‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:2461818/69‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:2461710/69‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:300669/69‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:3016910/69‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:3062489/69‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:3372417/69‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:365240/69‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:383692/69‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:1775170/69‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:1736567/69‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:1714123/69‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:1565692/69‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:1257841/69‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:675578/69‑1 (MQ=255)
                     aCGATGCAGTACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:714774/69‑1 (MQ=255)
                          gcagtaCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGAt                                          <  1:1022768/64‑1 (MQ=255)
                              taCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTa                                                 >  1:1390038/2‑53 (MQ=255)
                              taCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCgcag                                     >  1:2324544/2‑65 (MQ=255)
                              taCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCgcag                                     >  1:1860074/2‑65 (MQ=255)
                              taCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCgcag                                     >  1:1054877/2‑65 (MQ=255)
                              taCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCgcag                                     >  1:1856136/2‑65 (MQ=255)
                              taCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCgcag                                     >  1:1832666/2‑65 (MQ=255)
                              taCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCgcag                                     >  1:406454/2‑65 (MQ=255)
                              taCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCgcag                                     >  1:2038591/2‑65 (MQ=255)
                              taCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCgcag                                     >  1:419021/2‑65 (MQ=255)
                              taCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCgcag                                     >  1:22533/2‑65 (MQ=255)
                              taCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCgcag                                     >  1:2341553/2‑65 (MQ=255)
                              taCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCgcag                                     >  1:921839/2‑65 (MQ=255)
                              taCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCgcag                                     >  1:1411752/2‑65 (MQ=255)
                              taCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCgcag                                     >  1:2355007/2‑65 (MQ=255)
                              taCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCgcag                                     >  1:2379712/2‑65 (MQ=255)
                              taCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCgcag                                     >  1:663231/2‑65 (MQ=255)
                              taCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCgcag                                     >  1:2698214/2‑65 (MQ=255)
                              taCCGAAGGCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCgca                                      >  1:2279207/2‑64 (MQ=255)
                                      gCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTgg                         <  1:2009079/69‑1 (MQ=255)
                                      gCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTgg                         <  1:3416747/69‑1 (MQ=255)
                                      gCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTgg                         <  1:1905059/69‑1 (MQ=255)
                                      gCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTgg                         <  1:3364403/69‑1 (MQ=255)
                                      gCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTgg                         <  1:2448876/69‑1 (MQ=255)
                                      gCTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTgg                         <  1:2774255/69‑1 (MQ=255)
                                       cTGAATATCCGCGCCCATCACGGTAGGTATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTgg                         <  1:203022/68‑1 (MQ=255)
                                            tatCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTgcgc                   <  1:2938617/69‑1 (MQ=255)
                                                       aTCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACg        <  1:557699/69‑1 (MQ=255)
                                                          aCGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGac      >  1:2141494/1‑68 (MQ=255)
                                                          aCGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGac      >  1:1069477/1‑68 (MQ=255)
                                                          aCGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGac      >  1:1065481/1‑68 (MQ=255)
                                                             gTATGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGacaaca  >  1:3428839/1‑69 (MQ=255)
                                                             gTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGacaaca  >  1:2652290/1‑69 (MQ=255)
                                                             gTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGacaaca  >  1:2951275/1‑69 (MQ=255)
                                                               aGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTgg                         >  1:2008862/1‑44 (MQ=255)
                                                                  |                                                               
GCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACA  >  minE/1914308‑1914437

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: