Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 253,151 G→C 8.5% V79V (GTC→GTG yajB ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE253,1510GC 8.5% 193.4 / 6.7 70V79V (GTC→GTGyajBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/39);  new base C (2/4);  total (27/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAACCACTCCCGTGCTTCGCGGACTTCC  >  minE/253082‑253220
                                                                     |                                                                     
gACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGa                                                                      >  1:4604117/1‑71 (MQ=255)
  cAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACa                                                                    >  1:4022855/1‑71 (MQ=255)
            tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt                                                            <  1:1151928/69‑1 (MQ=255)
            tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt                                                            <  1:1820171/69‑1 (MQ=255)
            tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt                                                            <  1:1903508/69‑1 (MQ=255)
            tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt                                                            <  1:2980141/69‑1 (MQ=255)
            tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt                                                            <  1:1391312/69‑1 (MQ=255)
            tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt                                                            <  1:4225871/69‑1 (MQ=255)
            tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                                          <  1:3739492/71‑1 (MQ=255)
            tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                                          >  1:1838916/1‑71 (MQ=255)
              gCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                                          >  1:25562/1‑69 (MQ=255)
              gCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATa                                                        <  1:3161041/71‑1 (MQ=255)
              gCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATa                                                        <  1:1471233/71‑1 (MQ=255)
                ggAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAgg                                                      >  1:3452352/1‑71 (MQ=255)
                 gAAAATCGGGCGACAGCTGCGACCAGTGGCGGGACAGAAAGTGATCCCACATCACATCCAGCg                                                             <  1:3590256/63‑1 (MQ=255)
                 gAAAATCGGGCGACAGCTGCGACCAGTGGCGGGACAGAAAGTGATCCCACATCACATCCAGCg                                                             <  1:53627/63‑1 (MQ=255)
                 gAAAATCGGGCGACAGCTGCGACCAGTGGCGGGACAGAAAGTGATCCCACATCACATCCAGCg                                                             <  1:840047/63‑1 (MQ=255)
                  aaaaTCGGGCGACAGCTGCGACCAGTGGCGGGACAGAAAGTGATCCCACATCACATCCAGCGTAATAGgcg                                                    >  1:1193590/1‑71 (MQ=255)
                  aaaaTCGGGCGACAGCTGCGACCAGTGGCGGGACAGAAAGTGATCCCACATCACATCCAGCGTAATAGgcg                                                    >  1:3221462/1‑71 (MQ=255)
                  aaaGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt                                                            <  1:2964105/63‑1 (MQ=255)
                   aaGTCCGTCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGgcgc                                                   <  1:452907/71‑1 (MQ=255)
                   aaGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGgcgc                                                   <  1:3661874/71‑1 (MQ=255)
                     gTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                                          <  1:1387258/62‑1 (MQ=255)
                      tCCGGCGACAGCTGCGTCCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAAc                                                <  1:4210695/71‑1 (MQ=255)
                      tCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                                          <  1:2923251/61‑1 (MQ=255)
                        gggCGACAGCTGCGACCAGTGGCGGGACAGAAAGTGATCCCACATCACATCCAGCGTAATAGGCGCGAcac                                              <  1:538874/70‑3 (MQ=255)
                               aGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCaa                                                 <  1:4159342/61‑1 (MQ=255)
                                gCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt                                                            <  1:4645886/49‑1 (MQ=255)
                                 cTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGttt                                     <  1:4619550/71‑1 (MQ=255)
                                 cTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGttt                                     <  1:1920091/71‑1 (MQ=255)
                                   gCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTc                                    >  1:1516869/1‑70 (MQ=255)
                                     gACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCAc                                  >  1:753007/1‑70 (MQ=255)
                                     gACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACt                                 <  1:3199801/71‑1 (MQ=255)
                                      aCCAGTGGCGGGAAAGAAAGTGATCCCACATGACATc                                                                  <  1:802504/37‑1 (MQ=255)
                                        cAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt                                                            <  1:3180278/41‑1 (MQ=255)
                                         aGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCTTAATAGGCGCAACGCGGCGCGTTTCACTACGa                             <  1:2282900/71‑1 (MQ=255)
                                         aGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                                          >  1:818307/1‑42 (MQ=255)
                                         aGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCaa                                                 <  1:1233666/51‑1 (MQ=255)
                                           tGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                                          <  1:4292858/40‑1 (MQ=255)
                                           tGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                                          <  1:422384/40‑1 (MQ=255)
                                           tGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACt                                 <  1:3164160/65‑1 (MQ=255)
                                             gCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                                          >  1:1575056/1‑38 (MQ=255)
                                             gCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAAc                          >  1:655327/1‑70 (MQ=255)
                                               gggAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                                          >  1:4333926/1‑36 (MQ=255)
                                                ggAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                                          <  1:1415144/35‑1 (MQ=255)
                                                ggAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                                          <  1:606373/35‑1 (MQ=255)
                                                ggAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAAc                          <  1:452597/67‑1 (MQ=255)
                                                  aaagaaagTGATCCCACATGACATCCAGCGTAATAGGCGCa                                                  >  1:2959334/1‑41 (MQ=255)
                                                   aagaaagTGATCCCACATGACATCCAGCGTAATAGGCGCaa                                                 >  1:1868362/1‑41 (MQ=255)
                                                   aagaaagTGATCCCACATGACATCCAGCGTAATAGGCGCaa                                                 >  1:1396796/1‑41 (MQ=255)
                                                      aaagTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAACCACTCCCGTg                >  1:4596237/1‑71 (MQ=255)
                                                       aagTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAACCACTCCCGTGc               >  1:1506839/1‑71 (MQ=255)
                                                          tGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGcgcg                                        <  1:3150599/43‑1 (MQ=255)
                                                          tGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAACCACTCCCGTGCTTc            <  1:3481779/71‑1 (MQ=255)
                                                           gATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAACCACTCCCGTGCTTcg           >  1:4550452/1‑71 (MQ=255)
                                                           gATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAACCACTCCCGTGCTTcg           >  1:3776019/1‑71 (MQ=255)
                                                            aTCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAACCACTCCCGTGCTTcg           <  1:120503/70‑1 (MQ=255)
                                                              cccACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGa                             >  1:1156791/1‑50 (MQ=255)
                                                                cacaTGACATCCAGCGTAATAGGCGCAACGCGGCGCGttt                                     <  1:2514798/40‑1 (MQ=255)
                                                                cacaTGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGa                             <  1:1426850/48‑1 (MQ=255)
                                                                 acaTGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAACCAc                       >  1:2768526/1‑53 (MQ=255)
                                                                 acaTGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAACCACt                      <  1:1434039/54‑1 (MQ=255)
                                                                  caTGACATCCAGCGTAATAGGCGCAACGCGGCGCGttt                                     <  1:1190372/38‑1 (MQ=255)
                                                                  caTGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAACCACTCCCGTGCTTCGCGGACtt    <  1:1162786/71‑1 (MQ=255)
                                                                  caTGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAACCACTCCCGTGCTTCGCGGACtt    <  1:4151564/71‑1 (MQ=255)
                                                                   aTGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCa                                   >  1:3428453/1‑39 (MQ=255)
                                                                   aTGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGaaa                           >  1:4413632/1‑47 (MQ=255)
                                                                   aTGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAACCACTCCCGTGCTTCGCGGACTTc   >  1:3142482/1‑71 (MQ=255)
                                                                   aTGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAACCACTCCCGTGCTTCGCGGACTTc   >  1:1081955/1‑71 (MQ=255)
                                                                    tGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAACCACTCCCGTGCTTCGCGGACTTcc  >  1:4663208/1‑71 (MQ=255)
                                                                     |                                                                     
GACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAACCACTCCCGTGCTTCGCGGACTTCC  >  minE/253082‑253220

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: