Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 78,739 A→C 37.9% E278D (GAA→GAC 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,7390AC37.9% 61.0 / 107.4 116E278D (GAA→GACmurEUDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate:meso‑ diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (57/15);  new base C (10/34);  total (67/49)
Fisher's exact test for biased strand distribution p-value = 2.27e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

TTGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATT  >  minE/78675‑78802
                                                                |                                                               
ttGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGa                                                             <  1:2527065/69‑1 (MQ=255)
ttGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGa                                                             <  1:2134560/69‑1 (MQ=255)
ttGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGa                                                             <  1:2955001/69‑1 (MQ=255)
ttGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGa                                                             <  1:648582/69‑1 (MQ=255)
ttGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGa                                                             <  1:3314140/69‑1 (MQ=255)
ttGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGa                                                             <  1:485643/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:2335882/66‑1 (MQ=255)
      tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                                                          <  1:2532731/66‑1 (MQ=255)
      tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                                                          <  1:2656285/66‑1 (MQ=255)
      tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                                                          <  1:202632/66‑1 (MQ=255)
      tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                                                          <  1:3127287/66‑1 (MQ=255)
      tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                                                          <  1:1652857/66‑1 (MQ=255)
      tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                                                          <  1:1551625/66‑1 (MQ=255)
      tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                                                          <  1:1415905/66‑1 (MQ=255)
      tATCAATGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACt                                                          <  1:1512196/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)
          aaTGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCAc                                                     >  1:2307748/1‑67 (MQ=255)
          aaTGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCAc                                                     >  1:140527/1‑67 (MQ=255)
          aaTGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCAc                                                     >  1:982987/1‑67 (MQ=255)
          aaTGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCAc                                                     >  1:2078403/1‑67 (MQ=255)
          aaTGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCAc                                                     >  1:2035680/1‑67 (MQ=255)
          aaTGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCAc                                                     >  1:648595/1‑67 (MQ=255)
          aaTGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCAc                                                     >  1:1536722/1‑67 (MQ=255)
          aaTGGAAGATCATATTAATCCGAACTGTCAGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGAc                                                  >  1:1972871/1‑69 (MQ=255)
          aaTGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACTATCAc                                                     >  1:854083/1‑67 (MQ=255)
          aaTGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACTATCAc                                                     >  1:2673934/1‑67 (MQ=255)
          aaTGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACTATCAc                                                     >  1:832200/1‑67 (MQ=255)
          aaTGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACTATCAc                                                     >  1:800661/1‑67 (MQ=255)
          aaTGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACTATCAc                                                     >  1:360024/1‑67 (MQ=255)
          aaTGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACTATCAc                                                     >  1:1353118/1‑67 (MQ=255)
          aaTGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACTATCAc                                                     >  1:17513/1‑67 (MQ=255)
          aaTGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGACGTGAACTATCAc                                                     >  1:1978736/1‑67 (MQ=255)
          aaTGGAAGATCATATCAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACAGACGTGAACTATCAc                                                     >  1:259788/1‑67 (MQ=255)
              gAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTAt                                                        >  1:2337353/1‑60 (MQ=255)
              gAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTAt                                                        >  1:1485784/1‑60 (MQ=255)
              gAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTAt                                                        >  1:687971/1‑60 (MQ=255)
              gAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTAt                                                        >  1:524526/1‑60 (MQ=255)
              gAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTAt                                                        >  1:3429936/1‑60 (MQ=255)
              gAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTAt                                                        >  1:3274918/1‑60 (MQ=255)
                   tCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGc                                          >  1:2485892/1‑69 (MQ=255)
                    cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCTAAGTGAACTATCACGACAGCGGTGcg                                         >  1:442973/1‑69 (MQ=255)
                    cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcg                                         >  1:3325574/1‑69 (MQ=255)
                    cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcg                                         >  1:710699/1‑69 (MQ=255)
                    cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcg                                         >  1:1621233/1‑69 (MQ=255)
                    cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcg                                         >  1:2121828/1‑69 (MQ=255)
                    cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcg                                         >  1:919314/1‑69 (MQ=255)
                    cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcg                                         >  1:2714882/1‑69 (MQ=255)
                    cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcg                                         >  1:2705210/1‑69 (MQ=255)
                    cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcg                                         >  1:2174981/1‑69 (MQ=255)
                    cATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGc                                          >  1:518309/1‑68 (MQ=255)
                        ttAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcgacga                                     >  1:2542290/1‑69 (MQ=255)
                        ttAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGcgacg                                      >  1:2333493/1‑68 (MQ=255)
                              cGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGc                                >  1:7663/1‑68 (MQ=255)
                              cGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGc                                >  1:1897334/1‑68 (MQ=255)
                              cGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGc                                >  1:165386/1‑68 (MQ=255)
                              cGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGc                                >  1:3094674/1‑68 (MQ=255)
                                aaCTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCttt                             <  1:992882/69‑1 (MQ=255)
                                    gTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGc                          >  1:565074/1‑68 (MQ=255)
                                    gTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGc                          >  1:568928/1‑68 (MQ=255)
                                       aCGGACGCTGGTTGAAAGCGACCGACGTGAACTATCACTACAGCGGTGCGACCATTCGCTTTAGCTc                        <  1:1840372/67‑1 (MQ=255)
                                       aCGGACGCTGGTTGAAAGCGACCGACGTGAACTATCACGACAGCGGTGCGACCATTCGCTTTAGCTc                        <  1:2609511/67‑1 (MQ=255)
                                       aCGGACGCTGGTTGAAAGCGACCGACGTGAACTATCACGACAGCGGTGCGACCATTCGCTTTAGCTc                        <  1:1687030/67‑1 (MQ=255)
                                       aCGGACGCTGGTTGAAAGCGACCGACGTGAACTATCACGACAGCGGTGCGACCATTCGCTTTAGCTc                        <  1:588071/67‑1 (MQ=255)
                                       aCGGACGCTGGTTGAAAGCGACCGACGTGAACTATCACGACAGCGGTGCGACCATTCGCTTTAGCTc                        <  1:2574305/67‑1 (MQ=255)
                                       aCGGACGCTGGTTGAAAGCGACCGACGTGAACTATCACGACAGCGGTGCGACCATTCGCTTTAGCTc                        <  1:1518066/67‑1 (MQ=255)
                                       aCGGACGCTGGTTGAAAGCGACCGACGTGAACTATCACGACAGCGGTGCGACCATTCGCTTTAGCTc                        <  1:679383/67‑1 (MQ=255)
                                       aCGGACGCTGGTTGAAAGCGACCGACGTGAACTATCACGACAGCGGTGCGACCATTCGCTTTAGCTc                        <  1:3135840/67‑1 (MQ=255)
                                       aCGGACGCTGGTTGAAAGCGACCGACGTGAACTATCACGACAGCGGTGCGACCATTCGCTTTAGCTc                        <  1:3087737/67‑1 (MQ=255)
                                       aCGGACGCTGGTTGAAAGCGACCGACGTGAACTATCACGACAGCGGTGCGACCATTCGCTTTAGCTc                        <  1:715529/67‑1 (MQ=255)
                                       aCGGACGCTGGTTGAAAGCGACCGACGTGAACTATCACGACAGCGGTGCGACCATTCGCTTTAGCTc                        <  1:3041765/67‑1 (MQ=255)
                                       aCGGACGCTGGTTGAAAGCGACCGACGTGAACTATCACGACAGCGGTGCGACCATTCGCTTTAGCTc                        <  1:1321199/67‑1 (MQ=255)
                                       aCGGACGCTGGTTGAAAGCGACCGACGTGAACTATCACGACAGCGGTGCGACCATTCGCTTTAGCTc                        <  1:2854883/67‑1 (MQ=255)
                                       aCGGACGCTGGTTGAAAGCGACCGACGTGAACTATCACGACAGCGGTGCGACCATTCGCTTTAGCTc                        <  1:92580/67‑1 (MQ=255)
                                        cGGACGCTGGTTGAAAGCGACCGACGTGAACTATCACGACAGCGGTGCGACAATTCGCTTTAGCTc                        <  1:1746965/66‑1 (MQ=255)
                                              cTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTggg                >  1:88602/1‑68 (MQ=255)
                                              cTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTggg                >  1:359592/1‑68 (MQ=255)
                                              cTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTggg                >  1:2227373/1‑68 (MQ=255)
                                              cTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTggg                >  1:1967073/1‑68 (MQ=255)
                                              cTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTggg                >  1:1072180/1‑68 (MQ=255)
                                              cTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTggg                >  1:3043270/1‑68 (MQ=255)
                                              cTGGCTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTggg                >  1:2751619/1‑68 (MQ=255)
                                               tGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGc              <  1:2627904/69‑1 (MQ=255)
                                                ggTTGAGAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGt                    >  1:439786/1‑62 (MQ=255)
                                                ggTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGc              <  1:1205874/68‑1 (MQ=255)
                                                ggTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCg             >  1:260058/1‑69 (MQ=255)
                                                ggTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCg             >  1:2075583/1‑69 (MQ=255)
                                                ggTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCg             >  1:910938/1‑69 (MQ=255)
                                                ggTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCg             >  1:3141204/1‑69 (MQ=255)
                                                ggTTGAAAGCGACAGAAGTTAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGc              <  1:2511721/68‑1 (MQ=255)
                                                        gCGACCGACGTGAACTATCACGACAGCGGTGCGACCATTCGCTTTAGCTc                        <  1:2498265/50‑1 (MQ=255)
                                                        gCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGc              <  1:2508043/60‑1 (MQ=255)
                                                        gCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGaa     >  1:2808998/1‑69 (MQ=255)
                                                             cGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAAtt  <  1:29054/67‑1 (MQ=255)
                                                                |                                                               
TTGCGGTATCAATGGAAGATCATATTAATCCGAACTGTCACGGACGCTGGTTGAAAGCGACCGAAGTGAACTATCACGACAGCGGTGCGACGATTCGCTTTAGCTCAAGTTGGGGCGATGGCGAAATT  >  minE/78675‑78802

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: