Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 636,129 C→T 35.8% L1035L (CTC→CTT mukB → fused chromosome partitioning proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE636,1290CT35.8% 63.9 / 53.2 92L1035L (CTC→CTTmukBfused chromosome partitioning proteins
Reads supporting (aligned to +/- strand):  ref base C (47/12);  new base T (11/22);  total (58/34)
Fisher's exact test for biased strand distribution p-value = 1.52e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAGTACAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCTGATAGCGGGGCAGAAGAGCGGGCGCG  >  minE/636068‑636197
                                                             |                                                                    
gTCAGTACAACCAGGTGCTGGCTTCGCTGAAAAGCTCTTACGACACCAAAAAAGAGCTACTCAACGAt                                                                <  1:1527094/68‑1 (MQ=255)
gTCAGTACAACCAGGTGCTGGCTTCGCTGAAAAGCTCTTACGACACCAAAAAAGAGCTACTCAACGAt                                                                <  1:2218674/68‑1 (MQ=255)
      aCAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAAc                                                         >  1:2551610/1‑69 (MQ=255)
       cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCaa                                                          >  1:1318573/1‑67 (MQ=255)
       cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAAc                                                         >  1:50104/1‑68 (MQ=255)
       cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAAc                                                         >  1:1535514/1‑68 (MQ=255)
       cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAAc                                                         >  1:2013777/1‑68 (MQ=255)
       cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAAc                                                         >  1:2077272/1‑68 (MQ=255)
       cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAAc                                                         >  1:136581/1‑68 (MQ=255)
       cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAAc                                                         >  1:2542749/1‑68 (MQ=255)
       cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAAc                                                         >  1:301817/1‑68 (MQ=255)
       cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAAc                                                         >  1:3324272/1‑68 (MQ=255)
       cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAAc                                                         >  1:540733/1‑68 (MQ=255)
              gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGGTACTCAACGATCTGCAACGTGAAt                                                   >  1:408785/1‑67 (MQ=255)
              gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                                                   >  1:2862297/1‑67 (MQ=255)
              gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                                                   >  1:2838938/1‑67 (MQ=255)
              gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                                                   >  1:1036820/1‑67 (MQ=255)
              gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                                                   >  1:699264/1‑67 (MQ=255)
              gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                                                   >  1:3025441/1‑67 (MQ=255)
              gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                                                   >  1:743956/1‑67 (MQ=255)
              gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                                                   >  1:1777692/1‑67 (MQ=255)
              gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                                                   >  1:161815/1‑67 (MQ=255)
              gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                                                   >  1:1579557/1‑67 (MQ=255)
              gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                                                   >  1:800704/1‑67 (MQ=255)
              gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                                                   >  1:1501654/1‑67 (MQ=255)
                gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                                                      >  1:898674/1‑62 (MQ=255)
                gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                                                      >  1:2881893/1‑62 (MQ=255)
                gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                                                      >  1:284986/1‑62 (MQ=255)
                gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                                                      >  1:1315268/1‑62 (MQ=255)
                gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                                                      >  1:3295507/1‑62 (MQ=255)
                gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                                                      >  1:1627840/1‑62 (MQ=255)
                gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                                                      >  1:361032/1‑62 (MQ=255)
                gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                                                      >  1:2593319/1‑62 (MQ=255)
                gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                                                      >  1:1334843/1‑62 (MQ=255)
                gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                                                      >  1:1346331/1‑62 (MQ=255)
                gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                                                      >  1:1440977/1‑62 (MQ=255)
                gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                                                      >  1:1990609/1‑62 (MQ=255)
                gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAtag                                                 >  1:1399631/1‑65 (MQ=255)
                   ggCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGa                                            <  1:459351/69‑1 (MQ=255)
                   ggCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGa                                            <  1:1581685/69‑1 (MQ=255)
                   ggCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGa                                            <  1:1593343/69‑1 (MQ=255)
                   ggCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGa                                            <  1:3311121/69‑1 (MQ=255)
                   ggCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGa                                            <  1:2978201/69‑1 (MQ=255)
                   ggCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGa                                            <  1:120886/69‑1 (MQ=255)
                    gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCATGat                                           >  1:1095005/1‑69 (MQ=255)
                    gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGat                                           >  1:1312211/1‑69 (MQ=255)
                    gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGat                                           >  1:1512032/1‑69 (MQ=255)
                    gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGat                                           >  1:2607440/1‑69 (MQ=255)
                    gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGat                                           >  1:839064/1‑69 (MQ=255)
                    gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGat                                           >  1:2984542/1‑69 (MQ=255)
                    gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGat                                           >  1:479009/1‑69 (MQ=255)
                    gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGat                                           >  1:52364/1‑69 (MQ=255)
                    gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGat                                           >  1:1294995/1‑69 (MQ=255)
                    gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGat                                           >  1:2011681/1‑69 (MQ=255)
                    gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTACAGGat                                           >  1:1475444/1‑69 (MQ=255)
                         gCTGAAAAGCTCTTACGACACCAAAAAAGAGCTACTCAACGAt                                                                <  1:882936/43‑1 (MQ=255)
                          cTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCgg                                      <  1:2957303/68‑1 (MQ=255)
                          cTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCgg                                      <  1:3023847/68‑1 (MQ=255)
                            tctaaGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCATGATATCGgc                                     <  1:1873267/64‑1 (MQ=255)
                                  ttCTTACGCCACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGcgtgcgtg                              >  1:2171435/1‑68 (MQ=255)
                                  ttCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGcgtgcgtg                              >  1:3173977/1‑68 (MQ=255)
                                  ttCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGcgtgcgtg                              >  1:2983051/1‑68 (MQ=255)
                                  ttCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGcgtgcgtg                              >  1:2908984/1‑68 (MQ=255)
                                  ttCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGcgtgcgtg                              >  1:703296/1‑68 (MQ=255)
                                  ttCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGcgtgcgtg                              >  1:1603783/1‑68 (MQ=255)
                                  ttCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGcgtgcgtg                              >  1:2646841/1‑68 (MQ=255)
                                  ttCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGcgtgcgtg                              >  1:2608053/1‑68 (MQ=255)
                                  ttCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGcgtgcgtg                              >  1:2555252/1‑68 (MQ=255)
                                  ttCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGcgtgcgtg                              >  1:679262/1‑68 (MQ=255)
                                  ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc                                     <  1:1025855/61‑1 (MQ=255)
                                  ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc                                     <  1:1041518/61‑1 (MQ=255)
                                  ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc                                     <  1:2837595/61‑1 (MQ=255)
                                  ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc                                     <  1:1181821/61‑1 (MQ=255)
                                  ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc                                     <  1:3115649/61‑1 (MQ=255)
                                  ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc                                     <  1:2347757/61‑1 (MQ=255)
                                  ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc                                     <  1:1775682/61‑1 (MQ=255)
                                  ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc                                     <  1:192172/61‑1 (MQ=255)
                                  ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc                                     <  1:203745/61‑1 (MQ=255)
                                  ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc                                     <  1:2107206/61‑1 (MQ=255)
                                  ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc                                     <  1:2195483/61‑1 (MQ=255)
                                  ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc                                     <  1:2407312/61‑1 (MQ=255)
                                  ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc                                     <  1:2804967/61‑1 (MQ=255)
                                  ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc                                     <  1:1281798/61‑1 (MQ=255)
                                   tCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCt                            <  1:2795506/69‑1 (MQ=255)
                                   tCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCt                            <  1:1637681/69‑1 (MQ=255)
                                   tCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCt                            <  1:1656660/69‑1 (MQ=255)
                                   tCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCt                            <  1:1521702/69‑1 (MQ=255)
                                   tCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCt                            <  1:3231450/69‑1 (MQ=255)
                                   tCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCt                            <  1:1900215/69‑1 (MQ=255)
                                               aaaaaaGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCTGATAGCGGGGc                 >  1:1259117/1‑68 (MQ=255)
                                                       gCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc                                     <  1:1539031/40‑1 (MQ=255)
                                                        cTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc                                     <  1:1190665/39‑1 (MQ=255)
                                                             cAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCTGATAGCGGGGCAGAAGAGCGGGcgcg  <  1:1777376/69‑1 (MQ=255)
                                                             cAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCTGATAGCGGGGCAGAAGAGCGGGcgcg  <  1:717634/69‑1 (MQ=255)
                                                             cAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCTGATAGCGGGGCAGAAGAGCGGGcgcg  <  1:1380591/69‑1 (MQ=255)
                                                             cAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCTGATAGCGGGGCAGAAGAGCGGGcgcg  <  1:2990980/69‑1 (MQ=255)
                                                             |                                                                    
GTCAGTACAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCTGATAGCGGGGCAGAAGAGCGGGCGCG  >  minE/636068‑636197

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: