Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,403,432 A→G 63.6% I140I (ATT→ATC dapF ← diaminopimelate epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,403,4320AG63.6% 54.1 / 51.0 85I140I (ATT→ATCdapFdiaminopimelate epimerase
Reads supporting (aligned to +/- strand):  ref base A (18/13);  new base G (36/18);  total (54/31)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

TGAATCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCC  >  minE/2403364‑2403497
                                                                    |                                                                 
tGAATCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCata                                                                   <  1:863145/69‑1 (MQ=255)
tGAATCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCata                                                                   <  1:2715457/69‑1 (MQ=255)
tGAATCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCata                                                                   <  1:2800446/69‑1 (MQ=255)
tGAATCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCata                                                                   <  1:3145253/69‑1 (MQ=255)
tGAATCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCata                                                                   <  1:334088/69‑1 (MQ=255)
 gAATCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCataa                                                                  <  1:3305853/69‑1 (MQ=255)
 gAATCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCataa                                                                  <  1:1958624/69‑1 (MQ=255)
 gAATCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCataa                                                                  <  1:2029809/69‑1 (MQ=255)
 gAATCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCataa                                                                  <  1:2279861/69‑1 (MQ=255)
 gAATCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCataa                                                                  <  1:1317135/69‑1 (MQ=255)
 gAATCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCataa                                                                  <  1:843990/69‑1 (MQ=255)
 gAATCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCataa                                                                  <  1:1007425/69‑1 (MQ=255)
 gAATCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCataa                                                                  <  1:846005/69‑1 (MQ=255)
 gAATCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCataa                                                                  <  1:1197008/69‑1 (MQ=255)
 gAATCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCCCataa                                                                  <  1:72622/69‑1 (MQ=255)
    tCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCATAatat                                                               >  1:3005942/1‑69 (MQ=255)
    tCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCATAatat                                                               >  1:1530981/1‑69 (MQ=255)
    tCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCATAatat                                                               >  1:664685/1‑69 (MQ=255)
    tCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCATAatat                                                               >  1:178303/1‑69 (MQ=255)
              cGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCtttt                                                      >  1:373380/1‑68 (MQ=255)
              cGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCtttt                                                      >  1:805893/1‑68 (MQ=255)
              cGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCtttt                                                      >  1:189934/1‑68 (MQ=255)
              cGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTc                                                     >  1:1282259/1‑69 (MQ=255)
              cGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTc                                                     >  1:874314/1‑69 (MQ=255)
                   ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                >  1:2666315/1‑69 (MQ=255)
                   ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                >  1:1000210/1‑69 (MQ=255)
                   ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                >  1:2108569/1‑69 (MQ=255)
                   ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                >  1:2159124/1‑69 (MQ=255)
                   ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                >  1:2175336/1‑69 (MQ=255)
                   ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                >  1:2440756/1‑69 (MQ=255)
                   ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                >  1:355092/1‑69 (MQ=255)
                   ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                >  1:2823932/1‑69 (MQ=255)
                   ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                >  1:982255/1‑69 (MQ=255)
                   ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                >  1:3347166/1‑69 (MQ=255)
                   ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                >  1:1345605/1‑69 (MQ=255)
                   ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                >  1:1954590/1‑69 (MQ=255)
                   ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                >  1:1517806/1‑69 (MQ=255)
                   ttCCCATCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                >  1:122122/1‑69 (MQ=255)
                   ttCCCATCGACACCACGCCACATAAGATTCTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACtt                                                >  1:626887/1‑69 (MQ=255)
                        aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTccac                                                  >  1:1697805/1‑60 (MQ=255)
                        aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTcca                                                   >  1:2449386/1‑60 (MQ=255)
                        aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGtt                                            >  1:1337554/1‑68 (MQ=255)
                        aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGtt                                            >  1:1331406/1‑68 (MQ=255)
                        aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           >  1:2596683/1‑69 (MQ=255)
                        aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           >  1:3218376/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:3127892/1‑69 (MQ=255)
                        aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           >  1:3159866/1‑69 (MQ=255)
                        aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           >  1:2546491/1‑69 (MQ=255)
                        aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           >  1:2516793/1‑69 (MQ=255)
                        aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           >  1:2413749/1‑69 (MQ=255)
                        aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           >  1:1200624/1‑69 (MQ=255)
                        aTCGACACCACGCCACATAAGATTGTCTGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           >  1:791706/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)
                        aTCGACACAACGCCACATAAGAATGTCTGCTCGGCGGCGCGCATGAATTAGGTCt                                                         >  1:396916/1‑55 (MQ=255)
                           gACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAgc                                         >  1:345720/1‑68 (MQ=255)
                           gACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAgc                                         >  1:3020558/1‑68 (MQ=255)
                                cacGCCGCATAAGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGaaa                                   >  1:1130452/1‑69 (MQ=255)
                                   gccgcATAAGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAacgg                                <  1:1451380/69‑1 (MQ=255)
                                   gccgcATAAGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAacgg                                <  1:3054866/69‑1 (MQ=255)
                                           aGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAacggcacggc                          >  1:257105/1‑67 (MQ=255)
                                           aGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAacggcacggc                          >  1:2234258/1‑67 (MQ=255)
                                           aGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAacggcacggc                          >  1:2590816/1‑67 (MQ=255)
                                              ttGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAg                     <  1:336513/69‑1 (MQ=255)
                                              ttGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAg                     <  1:1105262/69‑1 (MQ=255)
                                              ttGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAg                     <  1:915114/69‑1 (MQ=255)
                                              ttGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAg                     <  1:1060580/69‑1 (MQ=255)
                                              ttGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAg                     <  1:3191070/69‑1 (MQ=255)
                                              ttGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAg                     <  1:1459513/69‑1 (MQ=255)
                                              ttGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAg                     <  1:3072584/69‑1 (MQ=255)
                                              ttGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAg                     <  1:2323888/69‑1 (MQ=255)
                                                   tGCTCGGCGGCTCGCATGATATAGGTCTTTTCCACTTTGTTa                                           <  1:202777/42‑1 (MQ=255)
                                                   tGCTCGGCGGCTCGCATGATATAGGTCTTTTCCACTTTGTTa                                           <  1:729958/42‑1 (MQ=255)
                                                   tGCTCGGCGGCGCGCTTGATATAGGTCTTTTCCACTTTGTTa                                           <  1:2860530/42‑1 (MQ=255)
                                                   tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           <  1:2075113/42‑1 (MQ=38)
                                                   tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           <  1:1189465/42‑1 (MQ=38)
                                                   tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           <  1:1903481/42‑1 (MQ=38)
                                                   tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           <  1:1812577/42‑1 (MQ=38)
                                                   tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           <  1:1152570/42‑1 (MQ=38)
                                                   tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           <  1:1953314/42‑1 (MQ=38)
                                                   tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           <  1:585684/42‑1 (MQ=38)
                                                   tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           <  1:1640199/42‑1 (MQ=38)
                                                   tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           <  1:1275393/42‑1 (MQ=38)
                                                   tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           <  1:1562143/42‑1 (MQ=38)
                                                   tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           <  1:1075029/42‑1 (MQ=38)
                                                   tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           <  1:3283970/42‑1 (MQ=38)
                                                   tGCTCGGCGGCGCGCATGATATAGGTCTTTTCCACTTTGTTa                                           <  1:2025570/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:396122/1‑69 (MQ=255)
                                                     cTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGaa              >  1:3421459/1‑69 (MQ=255)
                                                                 cATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGcc  <  1:1752164/69‑1 (MQ=255)
                                                                 cATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGcc  <  1:3181556/69‑1 (MQ=255)
                                                                 cATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGcc  <  1:3268897/69‑1 (MQ=255)
                                                                    |                                                                 
TGAATCACGCAATGCGGATTTCCCATCGACACCACGCCGCATAAGATTGTCTGCTCGGCGGCGCGCATAATATAGGTCTTTTCCGCTTTGTTAGCGCGAAACGGCACGGCGGAAGGTTCGAAGTTGGGTTCGCC  >  minE/2403364‑2403497

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