Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 697,873 C→T 28.1% R186R (CGC→CGT appA → phosphoanhydride phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE697,8730CT28.1% 83.7 / 31.2 78R186R (CGC→CGTappAphosphoanhydride phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (42/14);  new base T (15/7);  total (57/21)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGATCCTCAGCAGGGCAGGAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTGAGAAACAG  >  minE/697808‑697933
                                                                 |                                                            
gcgATCCTCAGCAGGGCAGGAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTcgcg                                                             <  1:1860416/67‑1 (MQ=255)
      cTCAGCAGGGCAGGAGGGTCAATTGCTGCCTTTACCGGGCATCGGCAAACGGCGTTTCGTGAGCTGGaa                                                     >  1:2617760/1‑69 (MQ=255)
      cTCAGCAGGGCAGGAGGGTCAATTGCTGACTTTACCGTGCATCGGCAAACGGCGTTTCGTGAACTGGaa                                                     >  1:1154129/1‑69 (MQ=255)
      cTCAGCAGGGCAGGAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGTGAACTGGaa                                                     >  1:1009985/1‑69 (MQ=255)
      cTCAGCAGGGCAGGAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGTGAACTGGaa                                                     >  1:3277844/1‑69 (MQ=255)
      cTCAGCAGGGCAGGAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGTGAACTGGaa                                                     >  1:3248346/1‑69 (MQ=255)
      cTCAGCAGGGCAGGAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGTGAACTGGaa                                                     >  1:1825873/1‑69 (MQ=255)
      cTCAGCAGGGCAGGAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGTGAACTGGaa                                                     >  1:1738859/1‑69 (MQ=255)
      cTCAGCAGGGCAGGAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGTGAACTGGaa                                                     >  1:3038045/1‑69 (MQ=255)
      cTCAGCAGGGCAGGAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGTGAACTGGaa                                                     >  1:2451728/1‑69 (MQ=255)
      cTCAGCAGGGCAGGAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGTGAACTGGaa                                                     >  1:552905/1‑69 (MQ=255)
      cTCAGCAGGGCAGGAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGTGAACTGGaa                                                     >  1:562904/1‑69 (MQ=255)
      cTCAGCAGGGCAGGAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGTGAACTGGaa                                                     >  1:1562131/1‑69 (MQ=255)
      cTCAGCAGGGCAGGAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGTGAACTGGaa                                                     >  1:27981/1‑69 (MQ=255)
      cTCAGCAGGGCAGGAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGTGAACTGGaa                                                     >  1:1477080/1‑69 (MQ=255)
      cTCAGCAGGGCAGGAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGTGAACTGGaa                                                     >  1:2728747/1‑69 (MQ=255)
            agggcaggAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTg                                               >  1:1562458/1‑69 (MQ=255)
            agggcaggAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTg                                               >  1:2737619/1‑69 (MQ=255)
             gggcaggAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGTGAACTGGAACGGGTGc                                              <  1:1635601/69‑1 (MQ=255)
             gggcaggAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGTGAACTGGAACGGGTGc                                              <  1:2612187/69‑1 (MQ=255)
             gggcaggAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGTGAACTGGAACGGGTGc                                              <  1:2546310/69‑1 (MQ=255)
             gggcaggAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGTGAACTGGAACGGGTGc                                              <  1:1525348/69‑1 (MQ=255)
             gggcaggAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGTGAACTGGAACGGGTGc                                              <  1:277904/69‑1 (MQ=255)
              ggcaggAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGTGAACTGGAACGGGTGc                                              <  1:2168122/68‑1 (MQ=255)
              ggcaggAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGTGAACTGGAACGGGTGc                                              <  1:2758656/68‑1 (MQ=255)
                caggAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTa                                           >  1:1367855/1‑69 (MQ=255)
                          aaTTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCa                                 >  1:1833067/1‑69 (MQ=255)
                          aaTTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCa                                 >  1:2015499/1‑69 (MQ=255)
                          aaTTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCa                                 >  1:2393372/1‑69 (MQ=255)
                          aaTTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCa                                 >  1:182909/1‑69 (MQ=255)
                          aaTTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCa                                 >  1:306764/1‑69 (MQ=255)
                          aaTTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCa                                 >  1:2248442/1‑69 (MQ=255)
                          aaTTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCa                                 >  1:598148/1‑69 (MQ=255)
                              gCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCa                             >  1:1687096/1‑69 (MQ=255)
                              gCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCa                             >  1:229482/1‑69 (MQ=255)
                              gCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCa                             >  1:2268021/1‑69 (MQ=255)
                              gCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCa                             >  1:930564/1‑69 (MQ=255)
                              gCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCa                             >  1:2027414/1‑69 (MQ=255)
                                             cATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTaaa              >  1:2238390/1‑69 (MQ=255)
                                             cATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTaaa              >  1:2990173/1‑69 (MQ=255)
                                             cATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTaaa              >  1:664666/1‑69 (MQ=255)
                                             cATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTaaa              >  1:802059/1‑69 (MQ=255)
                                                 ggaaaaCGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAAc             >  1:1954671/4‑66 (MQ=255)
                                                 ggCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTaaa              >  1:3026471/1‑65 (MQ=255)
                                                 ggCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAAc             >  1:2841262/1‑66 (MQ=255)
                                                 ggCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAAc             >  1:545352/1‑66 (MQ=255)
                                                 ggCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAAc             >  1:3029739/1‑66 (MQ=255)
                                                 ggCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAAc             >  1:462693/1‑66 (MQ=255)
                                                 ggCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAAc             >  1:3305688/1‑66 (MQ=255)
                                                 ggCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAAc             >  1:864610/1‑66 (MQ=255)
                                                 ggCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAAc             >  1:3157124/1‑66 (MQ=255)
                                                 ggCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAAc             >  1:108482/1‑66 (MQ=255)
                                                 ggCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAAc             >  1:1542124/1‑66 (MQ=255)
                                                 ggCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAAc             >  1:161187/1‑66 (MQ=255)
                                                 ggCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAAc             >  1:1650674/1‑66 (MQ=255)
                                                 ggCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAAc             >  1:16882/1‑66 (MQ=255)
                                                 ggCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAAc             >  1:1804185/1‑66 (MQ=255)
                                                 ggCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAAc             >  1:1876549/1‑66 (MQ=255)
                                                 ggCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAAc             >  1:1521531/1‑66 (MQ=255)
                                                 ggCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAAc             >  1:1486436/1‑66 (MQ=255)
                                                 ggCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAAc             >  1:1351451/1‑66 (MQ=255)
                                                 ggCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAAc             >  1:224419/1‑66 (MQ=255)
                                                    aaaCGGCGTTTCTCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTgaga       <  1:2462250/69‑1 (MQ=255)
                                                    aaaCGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTgaga       <  1:2088903/69‑1 (MQ=255)
                                                    aaaCGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTgaga       <  1:2173262/69‑1 (MQ=255)
                                                    aaaCGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTgaga       <  1:2217432/69‑1 (MQ=255)
                                                    aaaCGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTgaga       <  1:389924/69‑1 (MQ=255)
                                                    aaaCGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTgaga       <  1:390542/69‑1 (MQ=255)
                                                    aaaCGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTgaga       <  1:394520/69‑1 (MQ=255)
                                                    aaaCGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTgaga       <  1:1681848/69‑1 (MQ=255)
                                                    aaaCGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTgaga       <  1:892061/69‑1 (MQ=255)
                                                     aaCGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTgaga       <  1:959743/68‑1 (MQ=255)
                                                       cGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAAt                               >  1:2028109/1‑42 (MQ=255)
                                                       cGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTGAGAAAc    >  1:3211205/1‑69 (MQ=255)
                                                       cGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTGAGAAAc    >  1:1293687/1‑69 (MQ=255)
                                                         gCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTGAGAAACAg  <  1:2157583/69‑1 (MQ=255)
                                                         gCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTGAGAAACAg  <  1:2173438/69‑1 (MQ=255)
                                                         gCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTGAGAAACAg  <  1:2650501/69‑1 (MQ=255)
                                                            tttCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTgaga       <  1:2711023/61‑1 (MQ=255)
                                                                 |                                                            
GCGATCCTCAGCAGGGCAGGAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTGAGAAACAG  >  minE/697808‑697933

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