Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 636,108 C→T 33.2% Y1028Y (TAC→TAT mukB → fused chromosome partitioning proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE636,1080CT33.2% 77.2 / 48.4 93Y1028Y (TAC→TATmukBfused chromosome partitioning proteins
Reads supporting (aligned to +/- strand):  ref base C (47/15);  new base T (11/20);  total (58/35)
Fisher's exact test for biased strand distribution p-value = 2.43e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

GCAGCGCAGTTGAGTCAGTACAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCT  >  minE/636055‑636171
                                                     |                                                               
gcagcgcagTTGAGTCAGTACAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAagag                                                  <  1:2268510/69‑1 (MQ=255)
gcagcgcagTTGAGTCAGTACAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAagag                                                  <  1:1124719/69‑1 (MQ=255)
gcagcgcagTTGAGTCAGTACAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAagag                                                  <  1:671912/69‑1 (MQ=255)
             gTCAGTACAACCAGGTGCTGGCTTCGCTGAAAAGCTCTTACGACACCAAAAAAGAGCTACTCAACGAt                                      <  1:2218674/68‑1 (MQ=255)
             gTCAGTACAACCAGGTGCTGGCTTCGCTGAAAAGCTCTTACGACACCAAAAAAGAGCTACTCAACGAt                                      <  1:1527094/68‑1 (MQ=255)
                   aCAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAAc                               >  1:2551610/1‑69 (MQ=255)
                    cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGa                                                               >  1:2290889/1‑36 (MQ=255)
                    cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCaa                                >  1:1318573/1‑67 (MQ=255)
                    cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAAc                               >  1:2077272/1‑68 (MQ=255)
                    cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAAc                               >  1:2542749/1‑68 (MQ=255)
                    cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAAc                               >  1:136581/1‑68 (MQ=255)
                    cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAAc                               >  1:1535514/1‑68 (MQ=255)
                    cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAAc                               >  1:301817/1‑68 (MQ=255)
                    cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAAc                               >  1:3324272/1‑68 (MQ=255)
                    cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAAc                               >  1:50104/1‑68 (MQ=255)
                    cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAAc                               >  1:540733/1‑68 (MQ=255)
                    cAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAAc                               >  1:2013777/1‑68 (MQ=255)
                           gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGGTACTCAACGATCTGCAACGTGAAt                         >  1:408785/1‑67 (MQ=255)
                           gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                         >  1:2862297/1‑67 (MQ=255)
                           gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                         >  1:1777692/1‑67 (MQ=255)
                           gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                         >  1:161815/1‑67 (MQ=255)
                           gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                         >  1:2838938/1‑67 (MQ=255)
                           gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                         >  1:1036820/1‑67 (MQ=255)
                           gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                         >  1:3025441/1‑67 (MQ=255)
                           gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                         >  1:699264/1‑67 (MQ=255)
                           gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                         >  1:743956/1‑67 (MQ=255)
                           gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                         >  1:800704/1‑67 (MQ=255)
                           gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                         >  1:1501654/1‑67 (MQ=255)
                           gTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAAt                         >  1:1579557/1‑67 (MQ=255)
                             gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                            >  1:284986/1‑62 (MQ=255)
                             gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                            >  1:898674/1‑62 (MQ=255)
                             gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                            >  1:361032/1‑62 (MQ=255)
                             gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                            >  1:3295507/1‑62 (MQ=255)
                             gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                            >  1:1315268/1‑62 (MQ=255)
                             gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                            >  1:1334843/1‑62 (MQ=255)
                             gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                            >  1:2881893/1‑62 (MQ=255)
                             gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                            >  1:1346331/1‑62 (MQ=255)
                             gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                            >  1:2593319/1‑62 (MQ=255)
                             gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                            >  1:1440977/1‑62 (MQ=255)
                             gCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTg                            >  1:1627840/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:120886/69‑1 (MQ=255)
                                ggCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGa                  <  1:3311121/69‑1 (MQ=255)
                                ggCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGa                  <  1:1593343/69‑1 (MQ=255)
                                ggCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGa                  <  1:2978201/69‑1 (MQ=255)
                                 gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCATGat                 >  1:1095005/1‑69 (MQ=255)
                                 gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGat                 >  1:1294995/1‑69 (MQ=255)
                                 gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGat                 >  1:479009/1‑69 (MQ=255)
                                 gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGat                 >  1:2607440/1‑69 (MQ=255)
                                 gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGat                 >  1:1312211/1‑69 (MQ=255)
                                 gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGat                 >  1:52364/1‑69 (MQ=255)
                                 gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGat                 >  1:839064/1‑69 (MQ=255)
                                 gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGat                 >  1:2011681/1‑69 (MQ=255)
                                 gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGat                 >  1:2984542/1‑69 (MQ=255)
                                 gCTTCGCTGAAAAGTTCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGat                 >  1:1512032/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:2908984/1‑68 (MQ=255)
                                               ttCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGcgtgcgtg    >  1:2555252/1‑68 (MQ=255)
                                               ttCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGcgtgcgtg    >  1:2608053/1‑68 (MQ=255)
                                               ttCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGcgtgcgtg    >  1:2646841/1‑68 (MQ=255)
                                               ttCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGcgtgcgtg    >  1:3173977/1‑68 (MQ=255)
                                               ttCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGcgtgcgtg    >  1:679262/1‑68 (MQ=255)
                                               ttCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGcgtgcgtg    >  1:1603783/1‑68 (MQ=255)
                                               ttCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGcgtgcgtg    >  1:2983051/1‑68 (MQ=255)
                                               ttCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGcgtgcgtg    >  1:703296/1‑68 (MQ=255)
                                               ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc           <  1:1025855/61‑1 (MQ=255)
                                               ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc           <  1:1181821/61‑1 (MQ=255)
                                               ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc           <  1:1041518/61‑1 (MQ=255)
                                               ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc           <  1:2107206/61‑1 (MQ=255)
                                               ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc           <  1:1281798/61‑1 (MQ=255)
                                               ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc           <  1:2407312/61‑1 (MQ=255)
                                               ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc           <  1:203745/61‑1 (MQ=255)
                                               ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc           <  1:2195483/61‑1 (MQ=255)
                                               ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc           <  1:3115649/61‑1 (MQ=255)
                                               ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc           <  1:1775682/61‑1 (MQ=255)
                                               ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc           <  1:2837595/61‑1 (MQ=255)
                                               ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc           <  1:2804967/61‑1 (MQ=255)
                                               ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc           <  1:192172/61‑1 (MQ=255)
                                               ttCGTATGACACCAAAAAAGAGCTACTTAACGATCTGCAACGTGAATTGCAGGATATCGgc           <  1:2347757/61‑1 (MQ=255)
                                                tCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCt  <  1:1637681/69‑1 (MQ=255)
                                                tCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCt  <  1:1656660/69‑1 (MQ=255)
                                                tCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCt  <  1:3231450/69‑1 (MQ=255)
                                                tCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCt  <  1:1521702/69‑1 (MQ=255)
                                                tCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCt  <  1:1900215/69‑1 (MQ=255)
                                                tCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCt  <  1:2795506/69‑1 (MQ=255)
                                                     |                                                               
GCAGCGCAGTTGAGTCAGTACAACCAGGTGCTGGCTTCGCTGAAAAGTTCTTACGACACCAAAAAAGAGCTACTCAACGATCTGCAACGTGAATTGCAGGATATCGGCGTGCGTGCT  >  minE/636055‑636171

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: