Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 78,700 T→C 28.0% I265I (ATT→ATC murE → UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate:meso‑ diaminopimelate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE78,7000TC28.0% 90.5 / 37.4 89I265I (ATT→ATCmurEUDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate:meso‑ diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base T (39/25);  new base C (10/15);  total (49/40)
Fisher's exact test for biased strand distribution p-value = 9.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-01

TGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGA  >  minE/78639‑78767
                                                             |                                                                   
tGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGa                                                              <  1:2447306/69‑1 (MQ=255)
tGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGa                                                              <  1:21288/69‑1 (MQ=255)
tGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGa                                                              <  1:1656341/69‑1 (MQ=255)
  ggCCGCCGCTGGCTGTCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAAc                                                            <  1:2549785/69‑1 (MQ=255)
  ggCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAAc                                                            <  1:2138665/69‑1 (MQ=255)
  ggCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAAc                                                            <  1:1571844/69‑1 (MQ=255)
     cgccgcTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTTATCCGAACTGt                                                         <  1:46390/69‑1 (MQ=255)
     cgccgcTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGt                                                         <  1:1220632/69‑1 (MQ=255)
     cgccgcTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGt                                                         <  1:233159/69‑1 (MQ=255)
     cgccgcTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGt                                                         <  1:2049469/69‑1 (MQ=255)
     cgccgcTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGt                                                         <  1:1763523/69‑1 (MQ=255)
     cgccgcTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGt                                                         <  1:1297413/69‑1 (MQ=255)
                       cTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGaaa                                       >  1:1223592/1‑69 (MQ=255)
                       cTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGaaa                                       >  1:128070/1‑69 (MQ=255)
                       cTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGaaa                                       >  1:2614575/1‑69 (MQ=255)
                       cTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGaa                                        >  1:353790/1‑68 (MQ=255)
                       cTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGaa                                        >  1:3034559/1‑68 (MQ=255)
                         gCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGc                                     <  1:2328883/69‑1 (MQ=255)
                         gCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGc                                     <  1:726957/69‑1 (MQ=255)
                              aCGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCg                                <  1:494715/69‑1 (MQ=255)
                              aCGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCg                                <  1:3140414/69‑1 (MQ=255)
                              aCGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCg                                <  1:2652709/69‑1 (MQ=255)
                                    ttGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGa                          <  1:3314140/69‑1 (MQ=255)
                                    ttGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGa                          <  1:2134560/69‑1 (MQ=255)
                                    ttGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGa                          <  1:2527065/69‑1 (MQ=255)
                                    ttGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGa                          <  1:648582/69‑1 (MQ=255)
                                    ttGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGa                          <  1:485643/69‑1 (MQ=255)
                                    ttGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGa                          <  1:2955001/69‑1 (MQ=255)
                                     tGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGa                          <  1:2453648/68‑1 (MQ=255)
                                     tGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGa                          <  1:3182933/68‑1 (MQ=255)
                                          tttCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                       <  1:2989467/64‑1 (MQ=255)
                                          tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTTGTTGAAAGCGACCGACGTGAACt                       <  1:3220212/66‑1 (MQ=255)
                                          tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                       <  1:1128825/66‑1 (MQ=255)
                                          tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                       <  1:3127287/66‑1 (MQ=255)
                                          tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                       <  1:2656285/66‑1 (MQ=255)
                                          tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                       <  1:2335882/66‑1 (MQ=255)
                                          tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                       <  1:2532731/66‑1 (MQ=255)
                                          tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                       <  1:1652857/66‑1 (MQ=255)
                                          tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                       <  1:202632/66‑1 (MQ=255)
                                          tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                       <  1:1415905/66‑1 (MQ=255)
                                          tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                       <  1:1449135/66‑1 (MQ=255)
                                          tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                       <  1:1512196/66‑1 (MQ=255)
                                          tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                       <  1:2030250/66‑1 (MQ=255)
                                          tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                       <  1:1551625/66‑1 (MQ=255)
                                           aTCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                       <  1:1602527/65‑1 (MQ=255)
                                            tCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCAc                  >  1:3429244/1‑69 (MQ=255)
                                              aaTGGATGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACTATCAc                  >  1:643/1‑67 (MQ=255)
                                              aaTGGAAGATCATATTAATCCGAACTGTCACGGACGc                                                >  1:1495137/1‑37 (MQ=255)
                                              aaTGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCAc                  >  1:1832636/1‑67 (MQ=255)
                                              aaTGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCAc                  >  1:982987/1‑67 (MQ=255)
                                              aaTGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCAc                  >  1:1536722/1‑67 (MQ=255)
                                              aaTGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCAc                  >  1:140527/1‑67 (MQ=255)
                                              aaTGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCAc                  >  1:648595/1‑67 (MQ=255)
                                              aaTGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCAc                  >  1:870183/1‑67 (MQ=255)
                                              aaTGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCAc                  >  1:2307748/1‑67 (MQ=255)
                                              aaTGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCAc                  >  1:2660915/1‑67 (MQ=255)
                                              aaTGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCAc                  >  1:2078403/1‑67 (MQ=255)
                                              aaTGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCAc                  >  1:2269069/1‑67 (MQ=255)
                                              aaTGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCAc                  >  1:2035680/1‑67 (MQ=255)
                                              aaTGGAAGATCATATTAATCCGAACTGTCAGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGAc               >  1:1972871/1‑69 (MQ=255)
                                              aaTGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACTATCAc                  >  1:2673934/1‑67 (MQ=255)
                                              aaTGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACTATCAc                  >  1:1353118/1‑67 (MQ=255)
                                              aaTGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACTATCAc                  >  1:800661/1‑67 (MQ=255)
                                              aaTGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACTATCAc                  >  1:360024/1‑67 (MQ=255)
                                              aaTGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACTATCAc                  >  1:832200/1‑67 (MQ=255)
                                              aaTGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACTATCAc                  >  1:854083/1‑67 (MQ=255)
                                              aaTGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACTATCAc                  >  1:1978736/1‑67 (MQ=255)
                                              aaTGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACTATCAc                  >  1:17513/1‑67 (MQ=255)
                                              aaTGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACAGACGTGAACTATCAc                  >  1:259788/1‑67 (MQ=255)
                                                  gAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTAt                     >  1:1485784/1‑60 (MQ=255)
                                                  gAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTAt                     >  1:2337353/1‑60 (MQ=255)
                                                  gAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTAt                     >  1:687971/1‑60 (MQ=255)
                                                  gAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTAt                     >  1:524526/1‑60 (MQ=255)
                                                  gAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTAt                     >  1:3274918/1‑60 (MQ=255)
                                                  gAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTAt                     >  1:2916029/1‑60 (MQ=255)
                                                  gAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTAt                     >  1:3429936/1‑60 (MQ=255)
                                                  gAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTAt                     >  1:2566439/1‑60 (MQ=255)
                                                       tCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGc       >  1:2485892/1‑69 (MQ=255)
                                                       tCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGc       >  1:583417/1‑69 (MQ=255)
                                                        cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCTAAGTGAACTATCACGACAGCGGTGcg      >  1:442973/1‑69 (MQ=255)
                                                        cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcg      >  1:919314/1‑69 (MQ=255)
                                                        cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcg      >  1:2705210/1‑69 (MQ=255)
                                                        cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcg      >  1:710699/1‑69 (MQ=255)
                                                        cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcg      >  1:2174981/1‑69 (MQ=255)
                                                        cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcg      >  1:2121828/1‑69 (MQ=255)
                                                        cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcg      >  1:3325574/1‑69 (MQ=255)
                                                        cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcg      >  1:1621233/1‑69 (MQ=255)
                                                        cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcg      >  1:2714882/1‑69 (MQ=255)
                                                        cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGc       >  1:518309/1‑68 (MQ=255)
                                                            ttAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcgacga  >  1:2497210/1‑69 (MQ=255)
                                                            ttAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcgacga  >  1:2542290/1‑69 (MQ=255)
                                                            ttAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcgacg   >  1:2333493/1‑68 (MQ=255)
                                                             |                                                                   
TGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGA  >  minE/78639‑78767

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: