Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 791,603 G→A 59.5% G171S (GGC→AGC)  dadX → alanine racemase 2, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE791,6030GA59.5% 25.2 / 48.3 69G171S (GGC→AGC) dadXalanine racemase 2, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base G (10/18);  new base A (14/27);  total (24/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

GCGAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTC  >  minE/791537‑791667
                                                                  |                                                                
gCGAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGgc                                                                >  1:1741835/1‑69 (MQ=255)
gCGAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGgc                                                                >  1:1049851/1‑69 (MQ=255)
gCGAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGgc                                                                >  1:304899/1‑69 (MQ=255)
gCGAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGgc                                                                >  1:2201137/1‑69 (MQ=255)
gCGAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGgc                                                                >  1:3393809/1‑69 (MQ=255)
gCGAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGgc                                                                >  1:2828597/1‑69 (MQ=255)
gCGAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGgc                                                                >  1:1687584/1‑69 (MQ=255)
gCGAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGgc                                                                >  1:1379540/1‑69 (MQ=255)
  gAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCTAAGCGGAACATCCTAATGGAATTTCCGgcgc                                                              <  1:2623161/69‑1 (MQ=255)
  gAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGgcgc                                                              <  1:1932496/69‑1 (MQ=255)
  gAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGgcgc                                                              <  1:463621/69‑1 (MQ=255)
  gAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGgcgc                                                              <  1:3088939/69‑1 (MQ=255)
  gAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGgcgc                                                              <  1:2853622/69‑1 (MQ=255)
  gAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGgcgc                                                              <  1:996127/69‑1 (MQ=255)
  gAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGgcgc                                                              <  1:3047584/69‑1 (MQ=255)
            gAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGcgcg                                                     <  1:1654602/68‑1 (MQ=255)
            gAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGt                                                    <  1:3370823/69‑1 (MQ=255)
            gAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGt                                                    <  1:2984395/69‑1 (MQ=255)
            gAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGt                                                    <  1:1930245/69‑1 (MQ=255)
            gAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGt                                                    <  1:1288629/69‑1 (MQ=255)
            gAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGt                                                    <  1:867354/69‑1 (MQ=255)
            gAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGt                                                    <  1:785741/69‑1 (MQ=255)
            gAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGt                                                    <  1:1917042/69‑1 (MQ=255)
            gAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGt                                                    <  1:1038038/69‑1 (MQ=255)
            gAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGt                                                    <  1:987056/69‑1 (MQ=255)
            gAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGt                                                    <  1:1620215/69‑1 (MQ=255)
            gAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGt                                                    <  1:2243017/69‑1 (MQ=255)
            gAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGt                                                    <  1:1734894/69‑1 (MQ=255)
            gAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGt                                                    <  1:47255/69‑1 (MQ=255)
            gAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGt                                                    <  1:486227/69‑1 (MQ=255)
            gAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGt                                                    <  1:2044247/69‑1 (MQ=255)
            gAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGt                                                    <  1:1986982/69‑1 (MQ=255)
            gAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGt                                                    <  1:530790/69‑1 (MQ=255)
            gAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGt                                                    <  1:702862/69‑1 (MQ=255)
                           tCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAggcgg                                       >  1:2802052/1‑67 (MQ=255)
                                       gAAGCGTAACATCCTGATGGATTTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGAAGGGCTGGa                          <  1:1209027/68‑1 (MQ=255)
                                       gAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGAAGGGCTGGa                          <  1:757673/68‑1 (MQ=255)
                                       gAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGAAGGGCTGGa                          <  1:360336/68‑1 (MQ=255)
                                       gAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGAAGGGCTGGa                          <  1:551920/68‑1 (MQ=255)
                                       gAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGAAGGGCTGGa                          <  1:994911/68‑1 (MQ=255)
                                       gAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGAAGGGCTGGa                          <  1:998266/68‑1 (MQ=255)
                                       gAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGAAGGGCTGGa                          <  1:1001629/68‑1 (MQ=255)
                                       gAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGAAGGGCTGGa                          <  1:555895/68‑1 (MQ=255)
                                       gAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGAAGGGCTGGa                          <  1:1113885/68‑1 (MQ=255)
                                       gAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGAAGGGCTGGa                          <  1:1366444/68‑1 (MQ=255)
                                       gAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGAAGGGCTGGa                          <  1:2499255/68‑1 (MQ=255)
                                       gAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGAAGGGCTGGa                          <  1:1727833/68‑1 (MQ=255)
                                       gAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGAAGGGCTGGa                          <  1:2120013/68‑1 (MQ=255)
                                       gAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGAAGGGCTGGa                          <  1:193352/68‑1 (MQ=255)
                                       gAAGCGGAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGAAGGGCTGGa                          <  1:3007311/68‑1 (MQ=255)
                                            ggAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGaag                                 >  1:1984765/1‑54 (MQ=39)
                                            ggAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGaag                                 >  1:116046/1‑54 (MQ=39)
                                            ggAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGaag                                 >  1:1299520/1‑54 (MQ=39)
                                            ggAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGaag                                 >  1:1359274/1‑54 (MQ=39)
                                            ggAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGaag                                 >  1:1723489/1‑54 (MQ=39)
                                            ggAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGaag                                 >  1:1831860/1‑54 (MQ=39)
                                            ggAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGaag                                 >  1:638889/1‑54 (MQ=39)
                                            ggAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGaag                                 >  1:514389/1‑54 (MQ=39)
                                            ggAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGaag                                 >  1:2194098/1‑54 (MQ=39)
                                            ggAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGaag                                 >  1:2479974/1‑54 (MQ=39)
                                            ggAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGaag                                 >  1:2606594/1‑54 (MQ=39)
                                            ggAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGaag                                 >  1:3354382/1‑54 (MQ=39)
                                            ggAACATCCTGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGaag                                 >  1:3109484/1‑54 (MQ=39)
                                            ggAACATCCAGATGGAATTTCCAGCGCGATGGCGCGTATTGAGCAGGCGGCGGaag                                 >  1:1916194/1‑54 (MQ=38)
                                                  tCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGcgttc              >  1:2321493/1‑69 (MQ=255)
                                                              ttCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTc  <  1:3362768/69‑1 (MQ=255)
                                                              ttCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTc  <  1:878266/69‑1 (MQ=255)
                                                              ttCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTc  <  1:3125717/69‑1 (MQ=255)
                                                              ttCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTc  <  1:2634175/69‑1 (MQ=255)
                                                                  |                                                                
GCGAATGTTGGCGAAATGACCCTGATGTCGCATTTTGCCGAAGCGGAACATCCTGATGGAATTTCCGGCGCGATGGCGCGTATTGAGCAGGCGGCGGAGGGGCTGGAGTGTCGGCGTTCGTTGTCCAATTC  >  minE/791537‑791667

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: