Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,055,817 C→T 34.7% S137S (AGC→AGT ygiH → conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,055,8170CT34.7% 49.8 / 37.1 69S137S (AGC→AGTygiHconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (21/24);  new base T (12/12);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 8.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTGGCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAGTGCACTGATTGCTCCGTTTTATGTCTGGT  >  minE/2055757‑2055875
                                                            |                                                          
cATTGGCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACtc                                                    >  1:3400322/1‑69 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                         >  1:1401601/1‑60 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                         >  1:1405559/1‑60 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                         >  1:3332377/1‑60 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                         >  1:2528028/1‑60 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                         >  1:249900/1‑60 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                         >  1:2251928/1‑60 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                         >  1:475520/1‑60 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                         >  1:797030/1‑60 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACtcgtcg                                                >  1:3071646/1‑69 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACtcgtcg                                                >  1:2754235/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:142943/1‑69 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACtcgtcg                                                >  1:1828692/1‑69 (MQ=255)
    ggCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACacgtcg                                                >  1:2638422/1‑69 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                              >  1:2447108/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                              >  1:3135082/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                              >  1:1423011/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                              >  1:849409/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                              >  1:832350/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                              >  1:2273603/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                              >  1:2189688/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                              >  1:795715/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                              >  1:749122/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                              >  1:710416/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                              >  1:2077511/1‑68 (MQ=255)
       tGGGATCTCACTGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCt                                              >  1:174272/1‑68 (MQ=255)
        gggATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                         <  1:2094080/56‑1 (MQ=255)
        gggATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                         <  1:1441060/56‑1 (MQ=255)
        gggATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                         <  1:3329730/56‑1 (MQ=255)
        gggATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                         <  1:814160/56‑1 (MQ=255)
        gggATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                         <  1:238781/56‑1 (MQ=255)
        gggATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGa                                                         <  1:310409/56‑1 (MQ=255)
                  cGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCg                                       <  1:1177922/64‑1 (MQ=255)
                  cGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGt                                  <  1:836159/69‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                         <  1:530114/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                         <  1:3318899/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                         <  1:3187731/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                         <  1:3158805/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                         <  1:2645916/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                         <  1:145582/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                         <  1:1897561/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                         <  1:2740448/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                         <  1:2733523/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                         <  1:2039528/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                         <  1:2561117/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGTGGTTACTCGTCGCTGGGAg                                         <  1:212609/60‑1 (MQ=255)
                    gAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCg                                       <  1:2840428/62‑1 (MQ=255)
                     aGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAg                               <  1:947564/69‑1 (MQ=255)
                     aGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAg                               <  1:2056211/69‑1 (MQ=255)
                     aGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAg                               <  1:2531208/69‑1 (MQ=255)
                       tAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAGTg                             <  1:1680962/69‑1 (MQ=255)
                       tAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAGTg                             <  1:574679/69‑1 (MQ=255)
                       tAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAGTg                             <  1:1031254/69‑1 (MQ=255)
                         aTGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAGTGCa                           >  1:77883/1‑69 (MQ=255)
                              gggAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGc                                        <  1:2412437/51‑1 (MQ=255)
                              gggAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGc                                        <  1:448528/51‑1 (MQ=255)
                              gggAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGc                                        <  1:759022/51‑1 (MQ=255)
                                     tGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAGTGCACTGATTGc                   >  1:3402614/1‑65 (MQ=255)
                                      ggTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTc                                 <  1:2662297/50‑1 (MQ=255)
                                      ggTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTc                                 <  1:2772425/50‑1 (MQ=255)
                                      ggTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTc                                 <  1:17576/50‑1 (MQ=255)
                                      ggTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTc                                 <  1:1469088/50‑1 (MQ=255)
                                      ggTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTc                                 <  1:2793245/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)
                                                            |                                                          
CATTGGCTGGGATCTCACCGGAGTAATGGCGGGAACCTGGTTACTGACCGTGCTATTGAGCGGATACTCGTCGCTGGGAGCGATTGTCAGTGCACTGATTGCTCCGTTTTATGTCTGGT  >  minE/2055757‑2055875

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: