Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,723,207 G→C 43.7% L225L (CTC→CTG prkB ← predicted phosphoribulokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,723,2070GC43.7% 40.7 / 120.6 112L225L (CTC→CTGprkBpredicted phosphoribulokinase
Reads supporting (aligned to +/- strand):  ref base G (23/40);  new base C (14/35);  total (37/75)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

CAGCCAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGAAGTGTCGACGGTGGGAACGCGCTGGAAGTTAA  >  minE/2723147‑2723274
                                                            |                                                                   
cagcCAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACgg                                                             <  1:772318/69‑1 (MQ=255)
cagcCAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACgg                                                             <  1:1079620/69‑1 (MQ=255)
cagcCAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACgg                                                             <  1:2210596/69‑1 (MQ=255)
cagcCAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACgg                                                             <  1:1362923/69‑1 (MQ=255)
cagcCAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACgg                                                             <  1:437209/69‑1 (MQ=255)
cagcCAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACgg                                                             <  1:2253937/69‑1 (MQ=255)
cagcCAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACgg                                                             <  1:1319194/69‑1 (MQ=255)
cagcCAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACgg                                                             <  1:3269434/69‑1 (MQ=255)
cagcCAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACgg                                                             <  1:130729/69‑1 (MQ=255)
   ccAGGGGAAATCGATCCCTTCCAGATTACTAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGAt                                                          <  1:2942277/69‑1 (MQ=255)
   ccAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGAt                                                          <  1:572789/69‑1 (MQ=255)
   ccAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGAt                                                          <  1:1028924/69‑1 (MQ=255)
   ccAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGAt                                                          <  1:397302/69‑1 (MQ=255)
   ccAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGAt                                                          <  1:670292/69‑1 (MQ=255)
   ccAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGAt                                                          <  1:2632415/69‑1 (MQ=255)
   ccAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGAt                                                          <  1:2367433/69‑1 (MQ=255)
   ccAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGAt                                                          <  1:708513/69‑1 (MQ=255)
   ccAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGAt                                                          <  1:2814662/69‑1 (MQ=255)
   ccAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGAt                                                          <  1:2383713/69‑1 (MQ=255)
   ccAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGAt                                                          <  1:720767/69‑1 (MQ=255)
   ccAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGAt                                                          <  1:1317289/69‑1 (MQ=255)
   ccAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGAt                                                          <  1:858134/69‑1 (MQ=255)
   ccAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGAt                                                          <  1:98316/69‑1 (MQ=255)
   ccAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGAt                                                          <  1:1136814/69‑1 (MQ=255)
   ccAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGAt                                                          <  1:1109340/69‑1 (MQ=255)
    cAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGAt                                                          <  1:3184596/68‑1 (MQ=255)
        ggAAATCGATCCCTTCCAGAATACCAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCtt                                                     >  1:120348/1‑69 (MQ=255)
         gAAATCGATCCCTTCCAGATTACGAAAATGGATCAGCACAAAGCTTTCATCGAGCGACGGGATGCCttt                                                    >  1:805472/1‑69 (MQ=255)
         gAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACgg                                                             <  1:3208383/60‑1 (MQ=255)
         gAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGAt                                                          <  1:2665853/63‑1 (MQ=255)
          aaaTCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACgg                                                             <  1:2537134/59‑1 (MQ=255)
          aaaTCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACgg                                                             <  1:44824/59‑1 (MQ=255)
          aaaTCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACgg                                                             <  1:733399/59‑1 (MQ=255)
                 tCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCcgcg                                             <  1:2143125/68‑1 (MQ=255)
                 tCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCcgcg                                             <  1:1944953/68‑1 (MQ=255)
                 tCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCcgcg                                             <  1:362751/68‑1 (MQ=255)
                 tCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCcgcg                                             <  1:465971/68‑1 (MQ=255)
                 tCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCcgcg                                             <  1:2575101/68‑1 (MQ=255)
                 tCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCcgcg                                             <  1:1691995/68‑1 (MQ=255)
                 tCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCcgcg                                             <  1:713330/68‑1 (MQ=255)
                 tCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCcgcg                                             <  1:1647037/68‑1 (MQ=255)
                 tCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCcgcg                                             <  1:2846750/68‑1 (MQ=255)
                  cccTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCcgcg                                             <  1:29300/67‑1 (MQ=255)
                     ttCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGAc                                                               >  1:2553374/1‑46 (MQ=255)
                     ttCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCTTTTGCCGCGAACgg                                        >  1:2639362/1‑69 (MQ=255)
                     ttCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCTTTTGCCGCGAACgg                                        >  1:290715/1‑69 (MQ=255)
                     ttCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCTTTTGCCGCGAACgg                                        >  1:406409/1‑69 (MQ=255)
                     ttCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCTTTTGCCGCGAACgg                                        >  1:2486892/1‑69 (MQ=255)
                     ttCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCTTTTGCCGCGAACgg                                        >  1:2443627/1‑69 (MQ=255)
                     ttCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCTTTTGCCGCGAACgg                                        >  1:3269858/1‑69 (MQ=255)
                     ttCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCTTTTGCCGCGAACgg                                        >  1:216679/1‑69 (MQ=255)
                     ttCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCTTTTGCCGCGAACgg                                        >  1:1976104/1‑69 (MQ=255)
                     ttCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCTTTTGCCGCGAACgg                                        >  1:410404/1‑69 (MQ=255)
                     ttCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCTTTTGCCGCGAACgg                                        >  1:1837215/1‑69 (MQ=255)
                     ttCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCTTTTGCCGCGAACgg                                        >  1:1070102/1‑69 (MQ=255)
                      tCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACggg                                       <  1:3343104/69‑1 (MQ=255)
                           aTTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGa                                  <  1:2618427/69‑1 (MQ=255)
                           aTTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGa                                  <  1:964147/69‑1 (MQ=255)
                           aTTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGa                                  <  1:1957226/69‑1 (MQ=255)
                           aTTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGa                                  <  1:2631252/69‑1 (MQ=255)
                           aTTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCttt                                                    <  1:2777296/51‑1 (MQ=255)
                           aTTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCttt                                                    <  1:987572/51‑1 (MQ=255)
                           aTTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCttt                                                    <  1:1469491/51‑1 (MQ=255)
                           aTTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCttt                                                    <  1:1926581/51‑1 (MQ=255)
                           aTTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCttt                                                    <  1:1938793/51‑1 (MQ=255)
                           aTTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCttt                                                    <  1:1941749/51‑1 (MQ=255)
                           aTTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCttt                                                    <  1:2267362/51‑1 (MQ=255)
                           aTTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCttt                                                    <  1:2361824/51‑1 (MQ=255)
                           aTTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCttt                                                    <  1:236368/51‑1 (MQ=255)
                           aTTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCttt                                                    <  1:2882492/51‑1 (MQ=255)
                           aTTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCttt                                                    <  1:3039919/51‑1 (MQ=255)
                           aTTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCttt                                                    <  1:2553617/51‑1 (MQ=255)
                           aTTACGAAAATGGATCACCACAAAGCTTTCATCCAGCGACGGGATACCttt                                                    <  1:2710751/51‑1 (MQ=255)
                             tACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGAAg                                >  1:1745070/1‑69 (MQ=255)
                             tACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGAAg                                >  1:1699584/1‑69 (MQ=255)
                             tACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGAAg                                >  1:1968742/1‑69 (MQ=255)
                                      ggATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACg                                         >  1:3426952/1‑51 (MQ=255)
                                      ggATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACg                                         >  1:631433/1‑51 (MQ=255)
                                      ggATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACg                                         >  1:1664765/1‑51 (MQ=255)
                                      ggATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACg                                         >  1:2264454/1‑51 (MQ=255)
                                      ggATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACg                                         >  1:1573320/1‑51 (MQ=255)
                                      ggATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACg                                         >  1:2506528/1‑51 (MQ=255)
                                      ggATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACg                                         >  1:3183110/1‑51 (MQ=255)
                                             cacaAAGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGAAGTGTCGACGGTGGGAAc                >  1:1457529/1‑69 (MQ=255)
                                              acaAAGCTTTCATCGAGCGACGGGATACCTTTTGCcgcg                                             <  1:2454311/39‑1 (MQ=255)
                                              acaAAGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACggg                                       <  1:2852693/45‑1 (MQ=255)
                                                 aaGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGAAGTGTCGACGGTGGGAACgcgc            >  1:2250466/1‑69 (MQ=255)
                                                 aaGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGAAGTGTCGACGGTGGGAACgcgc            >  1:3347956/1‑69 (MQ=255)
                                                 aaGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGAAGTGTCGACGGTGGGAACgcgc            >  1:1196342/1‑69 (MQ=255)
                                                   gCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGAAGTGTCGACGGTGGGAACGCGCt           >  1:2896754/1‑68 (MQ=255)
                                                   gCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGAAGTGTCGACGGTGGGAACGCGCt           >  1:1648015/1‑68 (MQ=255)
                                                   gCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGAAGTGTCGACGGTGGGAACGCGCt           >  1:2734678/1‑68 (MQ=255)
                                                   gCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGAAGTGTCGACGGTGGGAACGCGCt           >  1:3100843/1‑68 (MQ=255)
                                                   gCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGAAGTGTCGACGGTGGGAACGCGCt           >  1:95371/1‑68 (MQ=255)
                                                   gCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGAAGTGTCGACGGTGGGAACGCGCt           >  1:2783751/1‑68 (MQ=255)
                                                    cTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACggg                                       <  1:1212459/39‑1 (MQ=255)
                                                           cGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGAAGTGTCGACGGTGGGAACGCGCTGGAAGTTaa  <  1:3152998/69‑1 (MQ=255)
                                                           cGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGAAGTGTCGACGGTGGGAACGCGCTGGAAGTTaa  <  1:155887/69‑1 (MQ=255)
                                                           cGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGAAGTGTCGACGGTGGGAACGCGCTGGAAGTTaa  <  1:922644/69‑1 (MQ=255)
                                                           cGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGAAGTGTCGACGGTGGGAACGCGCTGGAAGTTaa  <  1:114950/69‑1 (MQ=255)
                                                            |                                                                   
CAGCCAGGGGAAATCGATCCCTTCCAGATTACGAAAATGGATCACCACAAAGCTTTCATCGAGCGACGGGATACCTTTTGCCGCGAACGGGTTTGAAGTGTCGACGGTGGGAACGCGCTGGAAGTTAA  >  minE/2723147‑2723274

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: