Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,670,071 G→A 43.2% Q483Q (CAG→CAA malP → maltodextrin phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,670,0710GA43.2% 32.4 / 71.4 97Q483Q (CAG→CAAmalPmaltodextrin phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (33/22);  new base A (36/6);  total (69/28)
Fisher's exact test for biased strand distribution p-value = 6.68e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCACTGGCGGCTCTGTTGGATAAATCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGCAATATCGCGAGA  >  minE/2670011‑2670129
                                                            |                                                          
ggCACTGGCGGCTCTGTTGGATAAATCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCaa                                                    <  1:3074922/69‑1 (MQ=255)
ggCACTGGCGGCTCTGTTGGATAAATCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCaa                                                    <  1:635487/69‑1 (MQ=255)
    cTGGCGGCTCTGTTGGATAAATCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAAtat                                                 >  1:1052800/1‑66 (MQ=255)
                ttGGATAAATCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGaaaaa                                          <  1:968692/63‑1 (MQ=255)
                ttGGATAAATCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGaaaaa                                          <  1:2870701/63‑1 (MQ=255)
                ttGGATAAATCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGaaaaa                                          <  1:1721910/63‑1 (MQ=255)
                ttGGATAAATCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGaaaaa                                          <  1:1079934/63‑1 (MQ=255)
                ttGGATAAATCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGaaaaa                                          <  1:389496/63‑1 (MQ=255)
                 ttgATAAATCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGaaaaa                                          <  1:2927205/60‑1 (MQ=255)
                    aTAAATCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCtgat                                 <  1:2772844/68‑1 (MQ=255)
                    aTAAATCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCtgat                                 <  1:2192440/68‑1 (MQ=255)
                    aTAAATCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCtgat                                 <  1:2522555/68‑1 (MQ=255)
                      aaaTCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCtgatgac                              >  1:688535/1‑68 (MQ=255)
                      aaaTCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCtgatga                               >  1:1956939/1‑68 (MQ=255)
                      aaaTCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCtgatga                               >  1:2691415/1‑68 (MQ=255)
                      aaaTCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTgatgat                              >  1:750010/1‑69 (MQ=255)
                      aaaTCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTgatgat                              >  1:2323410/1‑69 (MQ=255)
                      aaaTCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTgatgat                              >  1:1551003/1‑69 (MQ=255)
                      aaaTCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTgatgat                              >  1:3243569/1‑69 (MQ=255)
                      aaaTCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTgatgat                              >  1:2024135/1‑69 (MQ=255)
                      aaaTCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTgatgat                              >  1:112086/1‑69 (MQ=255)
                      aaaTCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTgatgat                              >  1:2999630/1‑69 (MQ=255)
                      aaaTCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTgatgat                              >  1:1952942/1‑69 (MQ=255)
                      aaaTCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTgatgat                              >  1:950826/1‑69 (MQ=255)
                      aaaTCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTgatgat                              >  1:2825944/1‑69 (MQ=255)
                      aaaTCACTGCAAAAAGAGTGGGCTAACGATCACGATCAGCTGAGCAATCTGGaaaa                                           >  1:930625/1‑56 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:2365139/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:2733235/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:1014538/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:2379719/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:2435693/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:2742749/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:2784360/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:3097427/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:3170136/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:329216/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:355056/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:463110/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:76264/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:849794/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:878009/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:235889/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:1523850/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:1539425/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:1472339/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:1802426/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:2018250/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:106832/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:1043337/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:2234844/1‑58 (MQ=255)
                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                                       >  1:1057164/1‑58 (MQ=255)
                          cACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTCGCTGATGATGCGa                          >  1:2605103/1‑69 (MQ=255)
                          cACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTCGCTGATGATGCGa                          >  1:632710/1‑69 (MQ=255)
                          cACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTCGCTGATGATGCGa                          >  1:633831/1‑69 (MQ=255)
                          cACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTCGCTGATGATGCGa                          >  1:633875/1‑69 (MQ=255)
                          cACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTCGCTGATGATGCGa                          >  1:2965694/1‑69 (MQ=255)
                          cACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTCGCTGATGATGCGa                          >  1:2537978/1‑69 (MQ=255)
                          cACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTCGCTGATGATGCGa                          >  1:2256384/1‑69 (MQ=255)
                          cACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTCGCTGATGATGCGa                          >  1:1402699/1‑69 (MQ=255)
                          cACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTCGCTGATGATGCGa                          >  1:1914001/1‑69 (MQ=255)
                          cACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTCGCTGATGATGCGa                          >  1:1291795/1‑69 (MQ=255)
                          cACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTCGCTGATGATGCGa                          >  1:1195716/1‑69 (MQ=255)
                           aCTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGaa                         >  1:2792969/1‑69 (MQ=255)
                           aCTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGaa                         >  1:1505732/1‑69 (MQ=255)
                           aCTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGaa                         >  1:634628/1‑69 (MQ=255)
                           aCTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGaa                         >  1:3426043/1‑69 (MQ=255)
                           aCTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGaa                         >  1:2666556/1‑69 (MQ=255)
                           aCTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGaa                         >  1:1664573/1‑69 (MQ=255)
                           aCTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGaa                         >  1:2662469/1‑69 (MQ=255)
                           aCTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGaa                         >  1:587923/1‑69 (MQ=255)
                           aCTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGaa                         >  1:481061/1‑69 (MQ=255)
                             tGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAAt                       <  1:141997/69‑1 (MQ=255)
                                      agTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              <  1:569463/69‑1 (MQ=255)
                                      agTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              <  1:1000716/69‑1 (MQ=255)
                                      agTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              >  1:1195831/1‑69 (MQ=255)
                                      agTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              <  1:915280/69‑1 (MQ=255)
                                      agTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              >  1:2344881/1‑69 (MQ=255)
                                      agTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              <  1:1879969/69‑1 (MQ=255)
                                      agTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              >  1:2277431/1‑69 (MQ=255)
                                      agTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              >  1:2419905/1‑69 (MQ=255)
                                      agTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              >  1:2467031/1‑69 (MQ=255)
                                      agTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              <  1:2050577/69‑1 (MQ=255)
                                      agTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              <  1:2704532/69‑1 (MQ=255)
                                      agTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              <  1:270859/69‑1 (MQ=255)
                                      agTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              <  1:1989349/69‑1 (MQ=255)
                                      agTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              >  1:3372133/1‑69 (MQ=255)
                                      agTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              >  1:2898478/1‑69 (MQ=255)
                                      agTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              >  1:1690972/1‑69 (MQ=255)
                                      agTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              <  1:3012221/69‑1 (MQ=255)
                                      agTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              <  1:1658553/69‑1 (MQ=255)
                                      agTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              <  1:327664/69‑1 (MQ=255)
                                      agTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              <  1:3287720/69‑1 (MQ=255)
                                        tGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGCaa            <  1:1726621/69‑1 (MQ=255)
                                        tGGGCTAACGATCTCGATCAGCTGAGCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              <  1:156480/67‑1 (MQ=255)
                                                  aTCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGCAATATCGCgaga  <  1:1621636/69‑1 (MQ=255)
                                                    ctcGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGCAATATCGCgaga  >  1:71820/1‑67 (MQ=255)
                                                         tCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGc              <  1:2022957/50‑1 (MQ=255)
                                                            |                                                          
GGCACTGGCGGCTCTGTTGGATAAATCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAATTCCGTCAGCAATATCGCGAGA  >  minE/2670011‑2670129

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: