Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 253,106 C→G 6.7% P94P (CCG→CCC yajB ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE253,1060CG 6.7% 212.9 / 2.9 74P94P (CCG→CCCyajBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/38);  new base G (2/3);  total (33/41)
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.98e-01

CGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAAC  >  minE/253037‑253174
                                                                     |                                                                    
cGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCg                                                                     <  1:3205937/71‑1 (MQ=255)
   tGAGTCCGGCAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGc                                                                   >  1:785259/1‑70 (MQ=255)
    gAGTCCGGCAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGa                                                                 >  1:963467/1‑71 (MQ=255)
        ccGGCAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGc                                                             >  1:1116653/1‑71 (MQ=255)
         cGGCAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAg                                                              >  1:595818/1‑69 (MQ=255)
         cGGCAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCt                                                            >  1:4465088/1‑71 (MQ=255)
          ggCAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCt                                                            <  1:3148291/70‑1 (MQ=255)
           gCAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGc                                                          >  1:2770646/1‑71 (MQ=255)
           gCAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGc                                                          <  1:2015330/71‑1 (MQ=255)
            cAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGa                                                                 >  1:1567038/1‑63 (MQ=255)
             aaaaTCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGa                                                        >  1:4203690/1‑71 (MQ=255)
             aaaaTCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGa                                                        >  1:351072/1‑71 (MQ=255)
               aaTCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGTAAAGTCCGGCGACAGCTGCGAcc                                                      <  1:1412009/71‑1 (MQ=255)
               aaTCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTcc                                                                      <  1:875340/55‑1 (MQ=255)
                aTCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGAc                                                                <  1:3090839/60‑1 (MQ=255)
                   gTCATCACTTGCTCGCGGGCATAACAGACGAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGa                                                        >  1:1961871/1‑65 (MQ=255)
                      atcaCTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCg                                               >  1:2715597/1‑71 (MQ=255)
                      atcaCTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCg                                               >  1:924703/1‑71 (MQ=255)
                      atcaCTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCg                                               >  1:1220390/1‑71 (MQ=255)
                         aCTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACa                                                               >  1:4064456/1‑52 (MQ=255)
                         aCTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCggg                                             >  1:1631821/1‑70 (MQ=255)
                          cTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGc                                                                   <  1:771913/47‑1 (MQ=255)
                            tGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGa                                                                 >  1:519618/1‑47 (MQ=255)
                             gCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGtaaag                                         <  1:719261/70‑1 (MQ=255)
                             gCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGgaaag                                         >  1:3901398/1‑70 (MQ=255)
                             gCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGgaaag                                         >  1:2865832/1‑70 (MQ=255)
                             gCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGgaaag                                         >  1:4614777/1‑70 (MQ=255)
                              cTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGgaaaga                                        <  1:2626579/70‑1 (MQ=255)
                              cTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGgaaaga                                        <  1:4108210/70‑1 (MQ=255)
                              cTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGgaaaga                                        <  1:3089030/70‑1 (MQ=255)
                                 gcgGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGaaagaaag                                     <  1:866931/70‑1 (MQ=255)
                                  cgGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTg                                   <  1:4687668/71‑1 (MQ=255)
                                     gCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCggg                                             <  1:704097/58‑1 (MQ=255)
                                        tAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGgaaagaa                                       <  1:2747707/61‑1 (MQ=255)
                                          aCAGACAAATTCCTGTAGCGGAAAGTCCGGCGACa                                                               >  1:434620/1‑35 (MQ=255)
                                          aCAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCa                                                     >  1:350263/1‑45 (MQ=255)
                                          aCAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCcaca                           >  1:2702795/1‑71 (MQ=255)
                                          aCAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCcaca                           >  1:430510/1‑71 (MQ=255)
                                          aCAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCcaca                           >  1:212789/1‑71 (MQ=255)
                                          aCAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCcac                            <  1:315668/70‑1 (MQ=255)
                                          aCAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCcac                            <  1:1484874/70‑1 (MQ=255)
                                          aCAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCcac                            <  1:1033690/70‑1 (MQ=255)
                                           cAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACAt                          <  1:99663/71‑1 (MQ=255)
                                            aGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATc                                >  1:1689288/1‑64 (MQ=255)
                                             gACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGAc                                                       <  1:1697726/40‑1 (MQ=255)
                                             gACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCca                             >  1:4552861/1‑66 (MQ=255)
                                             gACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGa                        >  1:4604117/1‑71 (MQ=255)
                                               cAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACa                      >  1:4022855/1‑71 (MQ=255)
                                                aaaTTCCTGTAGCGGAAAGTCCGGCGACAGCTGAGAcc                                                      <  1:2095636/38‑1 (MQ=255)
                                                 aaTTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGgaaagaaa                                      <  1:1869637/53‑1 (MQ=255)
                                                         tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGaaagaaag                                     <  1:3250545/46‑1 (MQ=255)
                                                         tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt              <  1:1151928/69‑1 (MQ=255)
                                                         tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt              <  1:1391312/69‑1 (MQ=255)
                                                         tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt              <  1:1820171/69‑1 (MQ=255)
                                                         tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt              <  1:1903508/69‑1 (MQ=255)
                                                         tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt              <  1:4225871/69‑1 (MQ=255)
                                                         tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt              <  1:2980141/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:1471233/71‑1 (MQ=255)
                                                           gCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATa          <  1:3161041/71‑1 (MQ=255)
                                                            cGGAAAGTCCGGCGACAGCTGCGACCAGTGGCggg                                             <  1:2486253/35‑1 (MQ=255)
                                                             ggAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAgg        >  1:3452352/1‑71 (MQ=255)
                                                              gAAAATCGGGCGACAGCTGCGACCAGTGGCGGGACAGAAAGTGATCCCACATCACATCCAGCg               <  1:53627/63‑1 (MQ=255)
                                                              gAAAATCGGGCGACAGCTGCGACCAGTGGCGGGACAGAAAGTGATCCCACATCACATCCAGCg               <  1:3590256/63‑1 (MQ=255)
                                                              gAAAATCGGGCGACAGCTGCGACCAGTGGCGGGACAGAAAGTGATCCCACATCACATCCAGCg               <  1:840047/63‑1 (MQ=255)
                                                               aaaaTCGGGCGACAGCTGCGACCAGTGGCGGGACAGAAAGTGATCCCACATCACATCCAGCGTAATAGgcg      >  1:3221462/1‑71 (MQ=255)
                                                               aaaaTCGGGCGACAGCTGCGACCAGTGGCGGGACAGAAAGTGATCCCACATCACATCCAGCGTAATAGgcg      >  1:1193590/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)
                                                                    ccGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACAt                          >  1:336821/1‑46 (MQ=255)
                                                                     cGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATcc                               >  1:1792563/1‑40 (MQ=255)
                                                                     |                                                                    
CGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAAC  >  minE/253037‑253174

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: