Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,084,365 G→A 45.9% Q65Q (CAG→CAA yqjC → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,084,3650GA45.9% 15.9 / 70.3 87Q65Q (CAG→CAAyqjCconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/21);  new base A (20/20);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCGATCATCAGAAGAAAATCGCAAAGCAGAAAGATGAGGTGGCGGAACGCCAGCAAG  >  minE/2084300‑2084429
                                                                 |                                                                
cAGAATCGTATTGACTGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTg                                                               <  1:1198837/69‑1 (MQ=255)
cAGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTg                                                               <  1:458589/69‑1 (MQ=255)
cAGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTg                                                               <  1:954789/69‑1 (MQ=255)
 aGAATCGTATTGACGGTCTGATTAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGc                                                              <  1:233545/69‑1 (MQ=255)
 aGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATATCCAACTGc                                                              <  1:2446626/69‑1 (MQ=255)
 aGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGc                                                              <  1:2274373/69‑1 (MQ=255)
 aGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGc                                                              <  1:2978427/69‑1 (MQ=255)
 aGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGc                                                              <  1:444324/69‑1 (MQ=255)
 aGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGc                                                              <  1:197351/69‑1 (MQ=255)
 aGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGc                                                              <  1:1845880/69‑1 (MQ=255)
 aGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGc                                                              <  1:1975839/69‑1 (MQ=255)
 aGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGc                                                              <  1:340716/69‑1 (MQ=255)
 aGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGc                                                              <  1:1693881/69‑1 (MQ=255)
 aGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGc                                                              <  1:2230575/69‑1 (MQ=255)
 aGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGc                                                              <  1:1072628/69‑1 (MQ=255)
 aGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGc                                                              <  1:2328723/69‑1 (MQ=255)
 aGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGc                                                              <  1:1299831/69‑1 (MQ=255)
 aGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGc                                                              <  1:1262850/69‑1 (MQ=255)
 aGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGc                                                              <  1:3331708/69‑1 (MQ=255)
 aGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGc                                                              <  1:2930527/69‑1 (MQ=255)
 aGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGc                                                              <  1:3260569/69‑1 (MQ=255)
 aGAATCGTATTGACGGTCTGAATAAAGCCCCGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGc                                                              <  1:3105257/69‑1 (MQ=255)
  gAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGc                                                              <  1:2497469/68‑1 (MQ=255)
      cGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGcc                                                         <  1:828827/69‑1 (MQ=255)
      cGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGcc                                                         <  1:241475/69‑1 (MQ=255)
      cGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGcc                                                         <  1:2450940/69‑1 (MQ=255)
      cGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGcc                                                         <  1:3380478/69‑1 (MQ=255)
      cGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGcc                                                         <  1:729833/69‑1 (MQ=255)
      cGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGcc                                                         <  1:3206997/69‑1 (MQ=255)
      cGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGcc                                                         <  1:2647364/69‑1 (MQ=255)
      cGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGcc                                                         <  1:2624354/69‑1 (MQ=255)
      cGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGcc                                                         <  1:2586345/69‑1 (MQ=255)
       gTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCt                                                                >  1:1390091/1‑61 (MQ=255)
       gTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCg                                                        >  1:1236307/1‑69 (MQ=255)
       gTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCg                                                        >  1:1907184/1‑69 (MQ=255)
       gTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCg                                                        >  1:242780/1‑69 (MQ=255)
       gTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCg                                                        >  1:729940/1‑69 (MQ=255)
       gTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCg                                                        >  1:1063630/1‑69 (MQ=255)
       gTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCg                                                        >  1:3084471/1‑69 (MQ=255)
       gTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCg                                                        >  1:2033710/1‑69 (MQ=255)
       gTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCg                                                        >  1:491291/1‑69 (MQ=255)
       gTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCg                                                        >  1:1846234/1‑69 (MQ=255)
                     aaTAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCGATCATCagaagata                                          >  1:2431892/1‑67 (MQ=255)
                     aaTAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCGATCATCAGAAGaaa                                          >  1:3362194/1‑69 (MQ=255)
                     aaTAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCGATCATCAGAAGaaa                                          >  1:3364884/1‑69 (MQ=255)
                     aaTAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCGATCATCAGAAGaaa                                          >  1:3299692/1‑69 (MQ=255)
                     aaTAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCGATCATCAGAAGaaa                                          >  1:258267/1‑69 (MQ=255)
                     aaTAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCGATCATCAGAAGaaa                                          >  1:3369619/1‑69 (MQ=255)
                     aaTAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCGATCATCAGAAGaaa                                          >  1:1469735/1‑69 (MQ=255)
                     aaTAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCGATCATCAGAAGaaa                                          >  1:897562/1‑69 (MQ=255)
                     aaTAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCGATCATCAGAAGaaa                                          >  1:213995/1‑69 (MQ=255)
                     aaTAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCGATCATCAGAAGaaa                                          >  1:315984/1‑69 (MQ=255)
                     aaTAAAGCCCTGAGTGAAGTCCGGGCCAACTGTACAGATAGCCAGCTGCGTGCCGATCATCAGAAGaaa                                          >  1:1729647/1‑69 (MQ=255)
                         aaGCCCTGAGTGAAGTCCGGGCGAACTGTTCAGATAGCCAACTGCGTGCCGATCATCAGAAGaaaa                                         >  1:2315483/1‑66 (MQ=255)
                         aaGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGCTTGCCGATCATCAGAAGaaa                                          >  1:3372949/1‑65 (MQ=255)
                         aaGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGCGTGCCGATCATCAGAAGaaaa                                         >  1:3208301/1‑66 (MQ=255)
                         aaGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGCGTGCCGATCATCAGAAGaaaa                                         >  1:1034198/1‑66 (MQ=255)
                         aaGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGCGTGCCGATCATCAGAAGaaaa                                         >  1:3114978/1‑66 (MQ=255)
                         aaGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGCGTGCCGATCATCAGAAGaaaa                                         >  1:977163/1‑66 (MQ=255)
                         aaGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGCGTGCCGATCATCAGAAGaaaa                                         >  1:2989719/1‑66 (MQ=255)
                         aaGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGCGTGCCGATCATCAGAAGaaaa                                         >  1:1185075/1‑66 (MQ=255)
                         aaGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGCGTGCCGATCATCAGAAGaaaa                                         >  1:1692723/1‑66 (MQ=255)
                         aaGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGCGTGCCGATCATCAGAAGaaaa                                         >  1:1811315/1‑66 (MQ=255)
                         aaGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGCGTGCCGATCATCAGAAGaaaa                                         >  1:1940660/1‑66 (MQ=255)
                         aaGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGCGTGCCGATCATCAGAAGaaaa                                         >  1:2054961/1‑66 (MQ=255)
                         aaGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGCGTGCCGATCATCAGAAGaaaa                                         >  1:2196818/1‑66 (MQ=255)
                         aaGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGCGTGCCGATCATCAGAAGaaaa                                         >  1:2332543/1‑66 (MQ=255)
                         aaGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGCGTGCCGATCATCAGAAGaaaa                                         >  1:2351042/1‑66 (MQ=255)
                         aaGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGCGTGCCGATCATCAGAAGaaaa                                         >  1:2710749/1‑66 (MQ=255)
                         aaGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGCGTGCCGATCATCAGAAGaaaa                                         >  1:290388/1‑66 (MQ=255)
                         aaGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGCGTGCCGATCATCAGAAGaaaa                                         >  1:2931642/1‑66 (MQ=255)
                         aaGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGCGTGCCGATCATCAGAAGaaa                                          >  1:2354988/1‑65 (MQ=255)
                         aaGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAACTGCGTGCCGATCATCAGAAGaaa                                          >  1:2294917/1‑65 (MQ=255)
                                       gTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGcc                                                         <  1:1925922/36‑1 (MQ=255)
                                         ccGGGCCAACTTTTCAGATAGCCAGCTGCGTGCCGATCATCAAAAGAAAATCGCAAAGCAGAAAGATGa                      <  1:478287/69‑1 (MQ=255)
                                         ccGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCGATCATCAGAAGAAAATCGCAAAGCAGAAAGATGa                      <  1:2448540/69‑1 (MQ=255)
                                         ccGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCGATCATCAGAAGAAAATCGCAAAGCAGAAAGATGa                      <  1:841337/69‑1 (MQ=255)
                                         ccGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCGATCATCAGAAGAAAATCGCAAAGCAGAAAGATGa                      <  1:1287290/69‑1 (MQ=255)
                                         ccGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCGATCATCAGAAGAAAATCGCAAAGCAGAAAGATGa                      <  1:3140760/69‑1 (MQ=255)
                                            ggCCAACTGTTCAGATAGCCAGCTGCGTGCCGATCATCAGAAGAAAATCGCAAAGCAGAAAGATGAGGt                   <  1:681253/69‑1 (MQ=255)
                                            ggCCAACTGTTCAGATAGCCAGCTGCGTGCCGATCATCAGAAGAAAATCGCAAAGCAGAAAGATGAGGt                   <  1:104116/69‑1 (MQ=255)
                                              ccAACTGTTCAGATAGCCAGCTGCGTGCCGATCATCAGAAGAAAATCGCAAAGCAGAAAGATGa                      <  1:2566515/64‑1 (MQ=255)
                                                          aTAGCCAGCTGCGTGCCGATCATCAGAAGAAAATCGCAAAGCAGAAAGATGAGGTGGCGGAACGCCa       >  1:600618/1‑67 (MQ=255)
                                                           tAGCCAGCTGCGTGCCGATCATCAGAAGAAAATCGCAAAGCAGAAAGATGAGGTGGCGGAACGCCa       >  1:2686271/1‑66 (MQ=255)
                                                           tAGCCAGCTGCGTGCCGATCATCAGAAGAAAATCGCAAAGCAGAAAGATGAGGTGGCGGAACGCCa       >  1:3085968/1‑66 (MQ=255)
                                                           tAGCCAGCTGCGTGCCGATCATCAGAAGAAAATCGCAAAGCAGAAAGATGAGGTGGCGGAACGCCa       >  1:1099427/1‑66 (MQ=255)
                                                             gCCAGCTGCGTGCCGATCATCAGAAGAAAATCGCAAAGCAGAAAGATGAGGTGGCGGAACGCCAGCAAg  >  1:1599097/1‑69 (MQ=255)
                                                                 |                                                                
CAGAATCGTATTGACGGTCTGAATAAAGCCCTGAGTGAAGTCCGGGCCAACTGTTCAGATAGCCAGCTGCGTGCCGATCATCAGAAGAAAATCGCAAAGCAGAAAGATGAGGTGGCGGAACGCCAGCAAG  >  minE/2084300‑2084429

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: