Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,670,042 C→A 54.5% Q474K (CAA→AAA)  malP → maltodextrin phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,670,0420CA54.5% 3.3 / 86.0 77Q474K (CAA→AAA) malPmaltodextrin phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (27/8);  new base A (36/6);  total (63/14)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCCACGTCGCTGGATCAAACAGTGCAACCCGGCACTGGCGGCTCTGTTGGATAAATCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAAT  >  minE/2669979‑2670108
                                                               |                                                                  
aCCCCACGTCGCTGGATCAAACAGTGCAACCCGGCACTGGCGGCTCTGTTGGATAAATCACTGCaaa                                                                 >  1:2384613/1‑67 (MQ=255)
aCCCCACGTCGCTGGATCAAACAGTGCAACCCGGCACTGGCGGCTCTGTTGGATAAATCACTGCaaa                                                                 >  1:182841/1‑67 (MQ=255)
aCCCCACGTCGCTGGATCAAACAGTGCAACCCGGCACTGGCGGCTCTGTTGGATAAATCACTGCaaa                                                                 >  1:1080327/1‑67 (MQ=255)
                       gTGCAACCCGGCACTGGCGGCTCTGTTGGATAAATCACTGCAAAAAGAGTGGGCTAACGATCTCGATCa                                        <  1:969385/69‑1 (MQ=255)
                       gTGCAACCCGGCACTGGCGGCTCTGTTGGATAAATCACTGCAAAAAGAGTGGGCTAACGATCTCGATCa                                        <  1:2405236/69‑1 (MQ=255)
                                ggCACTGGCGGCTCTGTTGGATAAATCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCaa                               <  1:3074922/69‑1 (MQ=255)
                                ggCACTGGCGGCTCTGTTGGATAAATCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCaa                               <  1:635487/69‑1 (MQ=255)
                                    cTGGCGGCTCTGTTGGATAAATCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAAtat                            >  1:1052800/1‑66 (MQ=255)
                                                ttGGATAAATCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGaaaaa                     <  1:1079934/63‑1 (MQ=255)
                                                ttGGATAAATCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGaaaaa                     <  1:1721910/63‑1 (MQ=255)
                                                ttGGATAAATCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGaaaaa                     <  1:2870701/63‑1 (MQ=255)
                                                ttGGATAAATCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGaaaaa                     <  1:389496/63‑1 (MQ=255)
                                                ttGGATAAATCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGaaaaa                     <  1:968692/63‑1 (MQ=255)
                                                 ttgATAAATCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGaaaaa                     <  1:2927205/60‑1 (MQ=255)
                                                    aTAAATCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCtgat            <  1:2192440/68‑1 (MQ=255)
                                                    aTAAATCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCtgat            <  1:2522555/68‑1 (MQ=255)
                                                    aTAAATCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCtgat            <  1:2772844/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:1952942/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:950826/1‑69 (MQ=255)
                                                      aaaTCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTgatgat         >  1:2323410/1‑69 (MQ=255)
                                                      aaaTCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTgatgat         >  1:750010/1‑69 (MQ=255)
                                                      aaaTCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTgatgat         >  1:2999630/1‑69 (MQ=255)
                                                      aaaTCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTgatgat         >  1:2825944/1‑69 (MQ=255)
                                                      aaaTCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTgatgat         >  1:112086/1‑69 (MQ=255)
                                                      aaaTCACTGCAAAAAGAGTGGGCTAACGATCACGATCAGCTGAGCAATCTGGaaaa                      >  1:930625/1‑56 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:1043337/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:355056/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:2742749/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:2784360/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:878009/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:849794/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:76264/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:3097427/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:3170136/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:463110/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:329216/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:2018250/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:1523850/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:1472339/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:1539425/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:2733235/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:1802426/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:2234844/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:235889/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:2365139/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:106832/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:1014538/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:1057164/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:2435693/1‑58 (MQ=255)
                                                        aTCACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTc                  >  1:2379719/1‑58 (MQ=255)
                                                          cACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTCGCTGATGATGCGa     >  1:1914001/1‑69 (MQ=255)
                                                          cACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTCGCTGATGATGCGa     >  1:1195716/1‑69 (MQ=255)
                                                          cACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTCGCTGATGATGCGa     >  1:2537978/1‑69 (MQ=255)
                                                          cACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTCGCTGATGATGCGa     >  1:2965694/1‑69 (MQ=255)
                                                          cACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTCGCTGATGATGCGa     >  1:2605103/1‑69 (MQ=255)
                                                          cACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTCGCTGATGATGCGa     >  1:1291795/1‑69 (MQ=255)
                                                          cACTGAAAAAAGAGTGGGCTAACGATCTCGATCAACTGATCAATCTGGAAAAATTCGCTGATGATGCGa     >  1:1402699/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:2256384/1‑69 (MQ=255)
                                                           aCTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGaa    >  1:1664573/1‑69 (MQ=255)
                                                           aCTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGaa    >  1:1505732/1‑69 (MQ=255)
                                                           aCTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGaa    >  1:634628/1‑69 (MQ=255)
                                                           aCTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGaa    >  1:587923/1‑69 (MQ=255)
                                                           aCTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGaa    >  1:481061/1‑69 (MQ=255)
                                                           aCTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGaa    >  1:3426043/1‑69 (MQ=255)
                                                           aCTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGaa    >  1:2792969/1‑69 (MQ=255)
                                                           aCTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGaa    >  1:2662469/1‑69 (MQ=255)
                                                           aCTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGaa    >  1:2666556/1‑69 (MQ=255)
                                                             tGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAAt  <  1:141997/69‑1 (MQ=255)
                                                               |                                                                  
ACCCCACGTCGCTGGATCAAACAGTGCAACCCGGCACTGGCGGCTCTGTTGGATAAATCACTGCAAAAAGAGTGGGCTAACGATCTCGATCAGCTGATCAATCTGGAAAAATTCGCTGATGATGCGAAAT  >  minE/2669979‑2670108

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: