Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 86,263 C→T S253S (AGC→AGT murC → UDP‑N‑acetylmuramate:L‑alanine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE86,2630CT98.2% 168.8 / ‑3.8 54S253S (AGC→AGTmurCUDP‑N‑acetylmuramate:L‑alanine ligase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (30/23);  total (31/23)
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.33e-01

TTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACGCTGCTGC  >  minE/86196‑86330
                                                                   |                                                                   
ttGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGTGaa                                                                  <  1:2392174/71‑1 (MQ=255)
  gatgatCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGCg                                                                    >  1:2159064/1‑67 (MQ=255)
       tCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGTGAAGATGcc                                                            >  1:1053882/1‑70 (MQ=255)
       tCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTCACGGCTTCAGTGAAGATGcc                                                            >  1:1840687/1‑70 (MQ=255)
            tGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGTGAAGATGCCGAcgtg                                                      >  1:1028461/1‑71 (MQ=255)
             gATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGTGAAGATGCCGAcgtgc                                                     >  1:1374400/1‑71 (MQ=255)
                 cgcgAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGTGAAGATGCCGAcgtgcgtg                                                  <  1:2003655/70‑1 (MQ=255)
                  tcgAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGTGAAGATGCCGa                                                          >  1:1357397/2‑61 (MQ=255)
                   cgAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGTGAAGATGCCGAcgtgcgtg                                                  >  1:1341795/1‑68 (MQ=255)
                   cgAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGTGAAGATGCCGAcgtgcgtg                                                  >  1:1471532/1‑68 (MQ=255)
                   cgAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTag                                               <  1:1380652/71‑1 (MQ=255)
                   cgAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTag                                               <  1:1416803/71‑1 (MQ=255)
                   cgAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTag                                               <  1:2295719/71‑1 (MQ=255)
                    gAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTaga                                              <  1:2119582/71‑1 (MQ=255)
                       ttGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGTGAAGATg                                                              >  1:2009247/1‑52 (MQ=255)
                        tGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGAt                                          >  1:1852196/1‑71 (MQ=255)
                           tACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTAt                                       >  1:1612817/1‑71 (MQ=255)
                            aCCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTAt                                       >  1:1862770/1‑70 (MQ=255)
                               gcgAGTGGGGCGTCAGACCACGACTTACGGCTTCAGTGaa                                                                  >  1:999039/1‑40 (MQ=255)
                               gcgAGTGGGGCGTCAGACCACGACTTACGGCTTCAGTGAAGATGc                                                             <  1:1854836/45‑1 (MQ=255)
                                  aGTGGGGCGTCAGACCACGACTTACGGCTTCAGTGAAGATGCCGAcgt                                                       >  1:1671803/1‑48 (MQ=255)
                                   gTGGGGCGTCAGACCACGACTTACGGCTTCAGTGaa                                                                  <  1:2148884/36‑1 (MQ=255)
                                    tGGGGCGTCAGACCACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGAtt                              <  1:2085287/71‑1 (MQ=255)
                                    tGGGGCGTCAGACCACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGAtt                              <  1:1515534/71‑1 (MQ=255)
                                     ggggCGTCAGACCACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGAtt                              <  1:2082439/70‑1 (MQ=255)
                                          gTCAGACCACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCg                        >  1:522264/1‑71 (MQ=255)
                                           tCAGACCACGACTTACGGCTTCAGTGAAGATGCCGAc                                                         <  1:417704/37‑1 (MQ=255)
                                           tCAGACCACGACTTACGGCTTCAGTGAAGATGCCGAc                                                         <  1:1244613/37‑1 (MQ=255)
                                           tCAGACCACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCg                        >  1:2132366/1‑70 (MQ=255)
                                            cAGACCACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTg                             >  1:788571/1‑64 (MQ=255)
                                             aGACCACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGa                                >  1:1225079/1‑60 (MQ=255)
                                               aCCACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAggg                   >  1:1562511/1‑71 (MQ=255)
                                                 cACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTagaaga                                           >  1:424817/1‑45 (MQ=255)
                                                 cACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTagaaga                                           >  1:1113394/1‑45 (MQ=255)
                                                 cACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTagaaga                                           >  1:2210773/1‑45 (MQ=255)
                                                 cACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTagaaga                                           >  1:93066/1‑45 (MQ=255)
                                                 cACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCa                      <  1:2416600/66‑1 (MQ=255)
                                                 cACGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGc                 <  1:1746847/71‑1 (MQ=255)
                                                  aCGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCa                >  1:419675/1‑71 (MQ=255)
                                                   cGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGAt                               >  1:67420/1‑55 (MQ=255)
                                                   cGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGccc                         >  1:1744946/1‑61 (MQ=255)
                                                   cGACTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCAc               >  1:2047910/1‑71 (MQ=255)
                                                     aCTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACtt             <  1:1878659/71‑1 (MQ=255)
                                                      cTTACGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACtt             >  1:1133328/1‑70 (MQ=255)
                                                       ttACGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAg                     <  1:272253/61‑1 (MQ=255)
                                                         aCGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGAt                               <  1:710084/49‑1 (MQ=255)
                                                         aCGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGAt                               <  1:1492442/49‑1 (MQ=255)
                                                         aCGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGAt                               <  1:2184720/49‑1 (MQ=255)
                                                         aCGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACtt             <  1:1568033/67‑1 (MQ=255)
                                                         aCGGCTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACg         >  1:1505765/1‑71 (MQ=255)
                                                             cTTCAGTGAAGATGCCGACGTGCGTGTAGAAGATTAt                                       <  1:2241406/37‑1 (MQ=255)
                                                              ttCAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACgctgct    >  1:69365/1‑71 (MQ=255)
                                                                cAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACgctgc     <  1:1931697/68‑1 (MQ=255)
                                                                cAGTGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACGCtgctgc  >  1:248492/1‑71 (MQ=255)
                                                                   |                                                                   
TTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACGCTGCTGC  >  minE/86196‑86330

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: