Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,055,820 A→T 36.3% G138G (GGA→GGT ygiH → conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,055,8200AT36.3% 55.1 / 52.1 66G138G (GGA→GGTygiHconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (24/18);  new base T (12/12);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

CATTGGCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAGTGCACTGATTGCTCCGTTTTATGTCTGGTGGTTTAAGCCAC  >  minE/2055757‑2055887
                                                               |                                                                   
cATTGGCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACtc                                                                >  1:3400322/1‑69 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                                     >  1:249900/1‑60 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                                     >  1:1401601/1‑60 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                                     >  1:1405559/1‑60 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                                     >  1:2528028/1‑60 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                                     >  1:2251928/1‑60 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                                     >  1:3332377/1‑60 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                                     >  1:475520/1‑60 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                                     >  1:797030/1‑60 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACtcgtcg                                                            >  1:1828692/1‑69 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACtcgtcg                                                            >  1:318449/1‑69 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACtcgtcg                                                            >  1:2629734/1‑69 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACtcgtcg                                                            >  1:1823934/1‑69 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACtcgtcg                                                            >  1:2754235/1‑69 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACtcgtcg                                                            >  1:142943/1‑69 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACtcgtcg                                                            >  1:3071646/1‑69 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACacgtcg                                                            >  1:2638422/1‑69 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                                          >  1:174272/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                                          >  1:2447108/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                                          >  1:1423011/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                                          >  1:849409/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                                          >  1:2273603/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                                          >  1:2189688/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                                          >  1:832350/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                                          >  1:2077511/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                                          >  1:3135082/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                                          >  1:710416/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                                          >  1:749122/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                                          >  1:795715/1‑68 (MQ=255)
        gggATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                                     <  1:310409/56‑1 (MQ=255)
        gggATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                                     <  1:2094080/56‑1 (MQ=255)
        gggATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                                     <  1:3329730/56‑1 (MQ=255)
        gggATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                                     <  1:238781/56‑1 (MQ=255)
        gggATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                                     <  1:814160/56‑1 (MQ=255)
        gggATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                                     <  1:1441060/56‑1 (MQ=255)
                  cGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCg                                                   <  1:1177922/64‑1 (MQ=255)
                  cGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGt                                              <  1:836159/69‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                                     <  1:3318899/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                                     <  1:3187731/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                                     <  1:2740448/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                                     <  1:530114/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                                     <  1:3158805/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                                     <  1:2039528/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                                     <  1:145582/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                                     <  1:212609/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                                     <  1:1897561/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                                     <  1:2733523/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                                     <  1:2561117/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                                     <  1:2645916/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCg                                                   <  1:2840428/62‑1 (MQ=255)
                     aGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAg                                           <  1:2056211/69‑1 (MQ=255)
                     aGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAg                                           <  1:2531208/69‑1 (MQ=255)
                     aGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAg                                           <  1:947564/69‑1 (MQ=255)
                       tAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAGTg                                         <  1:1680962/69‑1 (MQ=255)
                       tAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAGTg                                         <  1:1031254/69‑1 (MQ=255)
                       tAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAGTg                                         <  1:574679/69‑1 (MQ=255)
                         aTGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAGTGCa                                       >  1:77883/1‑69 (MQ=255)
                              gggAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGc                                                    <  1:759022/51‑1 (MQ=255)
                              gggAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGc                                                    <  1:2412437/51‑1 (MQ=255)
                              gggAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGc                                                    <  1:448528/51‑1 (MQ=255)
                                     tGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAGTGCACTGATTGc                               >  1:3402614/1‑65 (MQ=255)
                                      ggTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTc                                             <  1:17576/50‑1 (MQ=255)
                                      ggTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTc                                             <  1:1469088/50‑1 (MQ=255)
                                      ggTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTc                                             <  1:2793245/50‑1 (MQ=255)
                                      ggTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTc                                             <  1:2662297/50‑1 (MQ=255)
                                      ggTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTc                                             <  1:2772425/50‑1 (MQ=255)
                                         tACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAGTGCACTGATTGCTCCGttt                        <  1:1845432/68‑1 (MQ=255)
                                                  tGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAGTGCACTGATTGCTCCGTTTTATGTCtggt              >  1:2147265/1‑69 (MQ=255)
                                                  tGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAGTGCACTGATTGCTCCGTTTTATGTCtggt              >  1:2961927/1‑69 (MQ=255)
                                                              gATACTCGTGGCTGGGAGCGATTGTCAGTGCACTGATTGCTCCGTTTTATGTCTGGTGGTTTAAGCcac  >  1:3120152/1‑69 (MQ=255)
                                                              gATACTCGTCGCTGGGAGCGATTGTCAGTGCACTTATTGCTCCGTTTTATGTCTGGTGGTTTAAGCcac  >  1:2937576/1‑69 (MQ=255)
                                                              gATACTCGTCGCTGGGAGCGATTGTCAGTGCACTGATTGCTCCGTTTTATGTCTGGTGGTTTAAGCcac  >  1:513265/1‑69 (MQ=255)
                                                              gATACTCGTCGCTGGGAGCGATTGTCAGTGCACTGATTGCTCCGTTTTATGTCTGGTGGTTTAAGCcac  >  1:452767/1‑69 (MQ=255)
                                                              gATACTCGTCGCTGGGAGCGATTGTCAGTGCACTGATTGCTCCGTTTTATGTCTGGTGGTTTAAGCcac  >  1:3381168/1‑69 (MQ=255)
                                                              gATACTCGTCGCTGGGAGCGATTGTCAGTGCACTGATTGCTCCGTTTTATGTCTGGTGGTTTAAGCcac  >  1:3283461/1‑69 (MQ=255)
                                                              gATACTCGTCGCTGGGAGCGATTGTCAGTGCACTGATTGCTCCGTTTTATGTCTGGTGGTTTAAGCcac  >  1:2599281/1‑69 (MQ=255)
                                                              gATACTCGTCGCTGGGAGCGATTGTCAGTGCACTGATTGCTCCGTTTTATGTCTGGTGGTTTAAGCcac  >  1:1725253/1‑69 (MQ=255)
                                                              gATACTCGTCGCTGGGAGCGATTGTCAGTGCACTGATTGCTCCGTTTTATGTCTGGTGGTTTAAGCcac  >  1:833678/1‑69 (MQ=255)
                                                              gATACTCGTCGCTGGGAGCGATTGTCAGTGCACTGATTGCTCCGTTTTATGTCTGGTGGTTTAAGCcac  >  1:1428454/1‑69 (MQ=255)
                                                              gATACTCGTCGCTGGGAGCGATTGTCAGTGCACTGATTGCTCCGTTTTATGTCTGGTGGTTTAAGCcac  >  1:963288/1‑69 (MQ=255)
                                                               |                                                                   
CATTGGCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAGTGCACTGATTGCTCCGTTTTATGTCTGGTGGTTTAAGCCAC  >  minE/2055757‑2055887

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