Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,279,256 C→G 15.8% A345A (GCC→GCG) ‡ yeeJ → adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,279,2560CG15.8% 196.9 / 25.8 82A345A (GCC→GCG) ‡yeeJadhesin
Reads supporting (aligned to +/- strand):  ref base C (41/28);  new base G (9/4);  total (50/32)
Fisher's exact test for biased strand distribution p-value = 5.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.06e-01

GCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATACCCCC  >  minE/1279186‑1279322
                                                                      |                                                                  
gCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGcc                                                                    >  1:4672137/1‑71 (MQ=255)
  tGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCt                                                                  >  1:5298725/1‑71 (MQ=255)
      tACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGtt                                                               <  1:5211454/70‑1 (MQ=255)
      tACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTc                                                              >  1:789441/1‑71 (MQ=255)
       aCCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTc                                                              >  1:3480601/1‑70 (MQ=255)
       aCCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTc                                                              >  1:276759/1‑70 (MQ=255)
         ccGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGtttg                                                             >  1:930567/1‑67 (MQ=255)
          cGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGtttgg                                                            >  1:582751/1‑66 (MQ=255)
          cGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGAt                                                           >  1:399332/1‑70 (MQ=255)
           gCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaa                                                         <  1:4646674/71‑1 (MQ=255)
            ccTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATaaa                                                        <  1:3471370/71‑1 (MQ=255)
               ggCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGAc                                                     >  1:622087/1‑71 (MQ=255)
                gCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaacat                                                     <  1:2447591/70‑4 (MQ=255)
                gCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGa                                                      <  1:5044865/69‑1 (MQ=255)
                 ccGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGa                                                   <  1:212857/71‑1 (MQ=255)
                   gCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATc                                                 >  1:5469915/1‑71 (MQ=255)
                   gCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATc                                                 <  1:1670750/71‑1 (MQ=255)
                     aCCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCgg                                               >  1:3774442/1‑71 (MQ=255)
                      ccTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCgg                                               >  1:1661023/1‑70 (MQ=255)
                       cTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGccc                                                                   >  1:2790648/1‑49 (MQ=255)
                       cTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATa                                                          <  1:5462014/58‑1 (MQ=255)
                       cTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATa                                                          <  1:4725793/58‑1 (MQ=255)
                           gCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAg                                          <  1:4415364/70‑1 (MQ=255)
                           gCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAg                                          <  1:5163512/70‑1 (MQ=255)
                            cGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTa                                        <  1:4158708/71‑1 (MQ=255)
                              gTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGt                                                                <  1:74015/45‑1 (MQ=255)
                              gTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGt                                                                <  1:3053862/45‑1 (MQ=255)
                              gTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAAt                                      >  1:4577140/1‑71 (MQ=255)
                               tAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTg                                                                 <  1:4469042/43‑1 (MQ=255)
                                aaaCTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTc                                                              >  1:4641011/1‑45 (MQ=255)
                                aaaCTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCaaa                                           >  1:3612501/1‑64 (MQ=255)
                                aaaCTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATcc                                    >  1:3678873/1‑71 (MQ=255)
                                aaaCTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATcc                                    >  1:4000099/1‑71 (MQ=255)
                                 aaCTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGtt                                                               >  1:1547429/1‑43 (MQ=255)
                                     ggTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATg                               >  1:5237960/1‑71 (MQ=255)
                                      gTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGc                              <  1:713950/71‑1 (MQ=255)
                                       tCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCa                                             <  1:5048629/55‑1 (MQ=255)
                                       tCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCa                                             <  1:4801790/55‑1 (MQ=255)
                                        cTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGc                              >  1:2481538/1‑69 (MQ=255)
                                         tATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGa                                                   <  1:4473728/47‑1 (MQ=255)
                                         tATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGa                                                   <  1:1009565/47‑1 (MQ=255)
                                         tATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTa                                        <  1:1886378/58‑1 (MQ=255)
                                           tGAACAGTATTATGGCTATGAAGTCGGGCTGTTTGGTaaa                                                        >  1:947842/1‑40 (MQ=16)
                                           tGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGc                                              >  1:132388/1‑50 (MQ=255)
                                           tGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATaa                         >  1:1902824/1‑71 (MQ=255)
                                           tGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATa                          <  1:4411527/70‑1 (MQ=255)
                                           tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                                                        >  1:247319/1‑40 (MQ=16)
                                           tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                                                        >  1:5562060/1‑40 (MQ=16)
                                           tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                                                        >  1:4522852/1‑40 (MQ=16)
                                           tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                                                        >  1:644916/1‑40 (MQ=16)
                                           tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGa                                                      >  1:547210/1‑42 (MQ=19)
                                           tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGa                                                      >  1:3504138/1‑42 (MQ=19)
                                            gAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACgct                                                  >  1:1865255/1‑43 (MQ=255)
                                            gAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCaaa                                           >  1:4982896/1‑52 (MQ=255)
                                            gAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCaaa                                           >  1:5495788/1‑52 (MQ=255)
                                            gAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCaaa                                           >  1:1861234/1‑52 (MQ=255)
                                             aaCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGa                                                   >  1:1065708/1‑43 (MQ=255)
                                             aaCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGAt                                                  >  1:1696726/1‑44 (MQ=255)
                                             aaCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGAt                                                  >  1:390284/1‑44 (MQ=255)
                                             aaCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAAcc                       >  1:3649832/1‑71 (MQ=255)
                                              aCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCg                      <  1:2886515/71‑1 (MQ=255)
                                                 gTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCa                                 <  1:3520479/57‑1 (MQ=255)
                                                   attatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCg                                                >  1:3319100/1‑40 (MQ=255)
                                                   attatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCg                                                >  1:1046219/1‑40 (MQ=255)
                                                    ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCa                                 >  1:1897872/1‑54 (MQ=255)
                                                    ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACGGCTgg                  <  1:2634890/69‑1 (MQ=255)
                                                    ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTgg                  <  1:690727/69‑1 (MQ=255)
                                                    ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTgg                  <  1:3714117/69‑1 (MQ=255)
                                                    ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGAc                <  1:409899/71‑1 (MQ=255)
                                                    ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGAc                >  1:3419483/1‑71 (MQ=255)
                                                     tatGGCCATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATaa                         >  1:4468066/1‑61 (MQ=255)
                                                         gCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAAc           <  1:795407/71‑1 (MQ=255)
                                                          cGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAAc           >  1:3659153/1‑70 (MQ=255)
                                                           gATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACta         >  1:260368/1‑71 (MQ=255)
                                                            aTGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTaa                                       <  1:4151459/40‑1 (MQ=255)
                                                             tGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACtata       <  1:1159844/71‑1 (MQ=255)
                                                              gAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCa                            >  1:5288081/1‑49 (MQ=255)
                                                              gAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCa                            >  1:1669836/1‑49 (MQ=255)
                                                              gAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCa                            >  1:773287/1‑49 (MQ=255)
                                                              gAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATAc      >  1:4278948/1‑71 (MQ=255)
                                                                 gTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATAcccc   >  1:3522912/1‑71 (MQ=255)
                                                                  tGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTc                                  >  1:1573677/1‑39 (MQ=255)
                                                                  tGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTg                   <  1:1941665/54‑1 (MQ=255)
                                                                  tGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATAccccc  >  1:5304496/1‑71 (MQ=255)
                                                                     cccTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACta         >  1:3542256/1‑61 (MQ=255)
                                                                      |                                                                  
GCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATACCCCC  >  minE/1279186‑1279322

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: