Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,403,448 G→A 52.6% A135V (GCG→GTG)  dapF ← diaminopimelate epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,403,4480GA52.6% 3.0 / 86.1 97A135V (GCG→GTG) dapFdiaminopimelate epimerase
Reads supporting (aligned to +/- strand):  ref base G (27/19);  new base A (33/18);  total (60/37)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

Only 100 of 101 total aligned reads displayed.

TTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCGGACCAG  >  minE/2403383‑2403515
                                                                 |                                                                   
ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                                  >  1:2159124/1‑69 (MQ=255)
ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                                  >  1:1000210/1‑69 (MQ=255)
ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                                  >  1:2440756/1‑69 (MQ=255)
ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                                  >  1:1954590/1‑69 (MQ=255)
ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                                  >  1:2666315/1‑69 (MQ=255)
ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                                  >  1:2823932/1‑69 (MQ=255)
ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                                  >  1:1517806/1‑69 (MQ=255)
ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                                  >  1:2108569/1‑69 (MQ=255)
ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                                  >  1:1345605/1‑69 (MQ=255)
ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                                  >  1:122122/1‑69 (MQ=255)
ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                                  >  1:2175336/1‑69 (MQ=255)
ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                                  >  1:3347166/1‑69 (MQ=255)
ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                                  >  1:355092/1‑69 (MQ=255)
ttCCCATCGACACCACGCCACATAAGATTCTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                                  >  1:626887/1‑69 (MQ=255)
     aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGtt                                                              >  1:1331406/1‑68 (MQ=255)
     aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGtt                                                              >  1:1337554/1‑68 (MQ=255)
     aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             >  1:3127892/1‑69 (MQ=255)
     aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             >  1:2596683/1‑69 (MQ=255)
     aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             >  1:2599015/1‑69 (MQ=255)
     aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             >  1:2747737/1‑69 (MQ=255)
     aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             >  1:2830041/1‑69 (MQ=255)
     aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             >  1:2995312/1‑69 (MQ=255)
     aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             >  1:2546491/1‑69 (MQ=255)
     aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             >  1:2516793/1‑69 (MQ=255)
     aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             >  1:3159866/1‑69 (MQ=255)
     aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             >  1:791706/1‑69 (MQ=255)
     aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             >  1:2413749/1‑69 (MQ=255)
     aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             >  1:3218376/1‑69 (MQ=255)
     aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             >  1:1200624/1‑69 (MQ=255)
     aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTACACTTTGTTa                                                             >  1:791870/1‑69 (MQ=255)
     aTCGACACCACGCCACATAAGATTGTCTGCTAGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             >  1:790327/1‑69 (MQ=255)
     aTCGACACCACGACACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTg                                                                >  1:3169480/1‑66 (MQ=255)
        gACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAgc                                                           >  1:3020558/1‑68 (MQ=255)
        gACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAgc                                                           >  1:345720/1‑68 (MQ=255)
             cacGCCGCATAAGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGaaa                                                     >  1:1130452/1‑69 (MQ=255)
                gccgcATAAGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAacgg                                                  <  1:3054866/69‑1 (MQ=255)
                gccgcATAAGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAacgg                                                  <  1:1451380/69‑1 (MQ=255)
                        aGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAacggcacggc                                            >  1:2590816/1‑67 (MQ=255)
                        aGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAacggcacggc                                            >  1:257105/1‑67 (MQ=255)
                        aGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAacggcacggc                                            >  1:2234258/1‑67 (MQ=255)
                           ttGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAg                                       <  1:3072584/69‑1 (MQ=255)
                           ttGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAg                                       <  1:2323888/69‑1 (MQ=255)
                           ttGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAg                                       <  1:3191070/69‑1 (MQ=255)
                           ttGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAg                                       <  1:336513/69‑1 (MQ=255)
                           ttGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAg                                       <  1:915114/69‑1 (MQ=255)
                           ttGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAg                                       <  1:1060580/69‑1 (MQ=255)
                           ttGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAg                                       <  1:1459513/69‑1 (MQ=255)
                           ttGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAg                                       <  1:1105262/69‑1 (MQ=255)
                                tGCTCGGCGGCTCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             <  1:729958/42‑1 (MQ=255)
                                tGCTCGGCGGCTCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             <  1:202777/42‑1 (MQ=255)
                                tGCTCGGCGGCGCGCTTGATATAGGTCTTTTCCACTTTGTTa                                                             <  1:2860530/42‑1 (MQ=255)
                                tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             <  1:3283970/42‑1 (MQ=38)
                                tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             <  1:1640199/42‑1 (MQ=38)
                                tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             <  1:1189465/42‑1 (MQ=38)
                                tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             <  1:1152570/42‑1 (MQ=38)
                                tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             <  1:1562143/42‑1 (MQ=38)
                                tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             <  1:1275393/42‑1 (MQ=38)
                                tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             <  1:1812577/42‑1 (MQ=38)
                                tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             <  1:1903481/42‑1 (MQ=38)
                                tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             <  1:1953314/42‑1 (MQ=38)
                                tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             <  1:2025570/42‑1 (MQ=38)
                                tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             <  1:585684/42‑1 (MQ=38)
                                tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             <  1:2075113/42‑1 (MQ=38)
                                tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             <  1:1075029/42‑1 (MQ=38)
                                 gCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             <  1:1901532/41‑1 (MQ=38)
                                  cTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                                             <  1:3274241/40‑1 (MQ=38)
                                  cTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGaa                                >  1:1791726/1‑69 (MQ=255)
                                  cTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGaa                                >  1:3421459/1‑69 (MQ=255)
                                  cTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGaa                                >  1:396122/1‑69 (MQ=255)
                                              cATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGcc                    <  1:1752164/69‑1 (MQ=255)
                                              cATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGcc                    <  1:3268897/69‑1 (MQ=255)
                                              cATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGcc                    <  1:3181556/69‑1 (MQ=255)
                                                   tataGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCata                >  1:2168000/1‑68 (MQ=255)
                                                   tataGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCata                >  1:3173202/1‑68 (MQ=255)
                                                   tataGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCata                >  1:2164129/1‑68 (MQ=255)
                                                   tataGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCata                >  1:2516388/1‑68 (MQ=255)
                                                         tCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTAcg          >  1:1481972/1‑68 (MQ=255)
                                                         tCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTAcg          >  1:2939712/1‑68 (MQ=255)
                                                           ttttCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCgg       >  1:1007983/1‑69 (MQ=255)
                                                           ttttCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCgg       >  1:1445259/1‑69 (MQ=255)
                                                           ttttCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCgg       >  1:1078330/1‑69 (MQ=255)
                                                           ttttCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCgg       >  1:1680935/1‑69 (MQ=255)
                                                             ttCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCgg       >  1:1484303/1‑67 (MQ=255)
                                                             ttCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCgg       >  1:2476589/1‑67 (MQ=255)
                                                             ttCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCgg       >  1:405998/1‑67 (MQ=255)
                                                             ttCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCgg       >  1:2323351/1‑67 (MQ=255)
                                                             ttCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCgg       >  1:2047085/1‑67 (MQ=255)
                                                             ttCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCgg       >  1:617605/1‑67 (MQ=255)
                                                             ttCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCgg       >  1:1104739/1‑67 (MQ=255)
                                                             ttCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCgg       >  1:6570/1‑67 (MQ=255)
                                                               ccGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCGGAc     <  1:1756511/67‑1 (MQ=255)
                                                               ccGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCGGAc     <  1:803092/67‑1 (MQ=255)
                                                               ccGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCGGAc     <  1:2154288/67‑1 (MQ=255)
                                                               ccGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCGGAc     <  1:1589659/67‑1 (MQ=255)
                                                               ccGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCGGAc     <  1:1600244/67‑1 (MQ=255)
                                                               ccGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCGGAc     <  1:1781455/67‑1 (MQ=255)
                                                                 gCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCGGACCAg  <  1:504263/68‑1 (MQ=255)
                                                                 gCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCGGACCAg  <  1:1111092/68‑1 (MQ=255)
                                                                 gCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCGGACCAg  <  1:3310183/68‑1 (MQ=255)
                                                                 gCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCGGACCAg  <  1:1420688/68‑1 (MQ=255)
                                                                 |                                                                   
TTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCCCATATTTACGCGGACCAG  >  minE/2403383‑2403515

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: