Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,057,322 G→A 29.3% G49G (GGG→GGA ttdA → L‑tartrate dehydratase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,057,3220GA29.3% 68.7 / 28.0 68G49G (GGG→GGAttdAL‑tartrate dehydratase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (8/40);  new base A (12/8);  total (20/48)
Fisher's exact test for biased strand distribution p-value = 8.58e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

GTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTG  >  minE/2057275‑2057387
                                               |                                                                 
gtggATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGAt                                               >  1:2176399/1‑68 (MQ=255)
gtggATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATg                                              >  1:3339071/1‑69 (MQ=255)
gtggATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATg                                              >  1:1112123/1‑69 (MQ=255)
gtggATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATg                                              >  1:2133974/1‑69 (MQ=255)
gtggATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATg                                              <  1:2002868/69‑1 (MQ=255)
gtggATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATg                                              >  1:1971809/1‑69 (MQ=255)
gtggATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATg                                              <  1:3313533/69‑1 (MQ=255)
gtggATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATg                                              <  1:1078419/69‑1 (MQ=255)
gtggATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATg                                              >  1:2681250/1‑69 (MQ=255)
gtggATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATg                                              >  1:2678908/1‑69 (MQ=255)
      aaaCTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTc                                           <  1:2646539/66‑1 (MQ=255)
      aaaCTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTc                                           <  1:2186976/66‑1 (MQ=255)
      aaaCTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTc                                           <  1:1729130/66‑1 (MQ=255)
          tAAAACAGCTAATGGATGCCGAAACATCGTCGATGGGAAAAATCATCTACCACACGATGTTc                                           >  1:145860/1‑62 (MQ=255)
          tAAAACAGCTAAAGGATGCCGAAACATCGTCGATGGGAAAAATCATCTACCACACGATGtt                                            >  1:2314110/1‑61 (MQ=255)
          tAAAACAGCTAAAGGATGCCGAAACATCGTCGATGGGAAAAATCATCTACCACACGATGTTc                                           >  1:2715345/1‑62 (MQ=255)
          tAAAACAGCTAAAGGATGCCGAAACATCGTCGATGGGAAAAATCATCTACCACACGATGTTc                                           >  1:2932323/1‑62 (MQ=255)
          tAAAACAGCTAAAGGATGCCGAAACATCGTCGATGGGAAAAATCATCTACCACACGATGTTc                                           >  1:3385945/1‑62 (MQ=255)
          tAAAACAGCTAAAGGATGCCGAAACATCGTCGATGGGAAAAATCATCTACCACACGATGTTc                                           >  1:1192590/1‑62 (MQ=255)
          tAAAACAGCTAAAGGATGCCGAAACATCGTCGATGGGAAAAATCATCTACCACACGATGTTc                                           >  1:619726/1‑62 (MQ=255)
          tAAAACAGCTAAAGGATGCCGAAACATCGTCGATGGGAAAAATCATCTACCACACGATGTTc                                           >  1:2704720/1‑62 (MQ=255)
          tAAAACAGCTAAAGGATGCCGAAACATCGTCGATGGGAAAAATCATCTACCACACGATGTTc                                           >  1:2350392/1‑62 (MQ=255)
          tAAAACAGCTAAAGGATGCCGAAACATCGTCGATGGGAAAAATCATCTACCACACGATGTTc                                           >  1:701881/1‑62 (MQ=255)
          tAAAACAGCTAAAGGATGCCGAAACATCGTCGATGGGAAAAATCATCTACCACACGATGTTc                                           >  1:448849/1‑62 (MQ=255)
          tAAAACAGCTAAAGGATGCCGAAACATCGTCGATGGGAAAAATCATCTACCACACGATGTTc                                           >  1:3243969/1‑62 (MQ=255)
                       ggATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACGACATGCAAAAAGCGa                        <  1:2295570/68‑1 (MQ=255)
                       ggATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGa                        <  1:2739394/68‑1 (MQ=255)
                       ggATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGa                        <  1:3171202/68‑1 (MQ=255)
                       ggATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGa                        <  1:3169532/68‑1 (MQ=255)
                       ggATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGa                        <  1:2263454/68‑1 (MQ=255)
                       ggATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGa                        <  1:2210459/68‑1 (MQ=255)
                       ggATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGa                        <  1:2137091/68‑1 (MQ=255)
                       ggATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGa                        <  1:1945227/68‑1 (MQ=255)
                       ggATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGa                        <  1:923352/68‑1 (MQ=255)
                       ggATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGa                        <  1:753770/68‑1 (MQ=255)
                       ggATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGa                        <  1:1750540/68‑1 (MQ=255)
                       ggATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGa                        <  1:676031/68‑1 (MQ=255)
                       ggATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGa                        <  1:163602/68‑1 (MQ=255)
                        gATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGa                        <  1:3374817/67‑1 (MQ=255)
                         aTGCCGAAACATCGTCGATGGGAAAAATCATCTACCACACGATGTTCGATAACATGCAAAAAGCGAtcg                     <  1:2292540/69‑3 (MQ=255)
                         aTGCCGAAACATCGTCGATGGGAAAAATCATCTACCACACGATGTTCGATAACATGCAAAAAGCGAtcg                     <  1:986892/69‑3 (MQ=255)
                         aTGCCGAAACATCGTCGATGGGAAAAATCATCTACCACACGATGTTCGATAACATGCAAAAAGCGAtcg                     <  1:1235983/69‑3 (MQ=255)
                         aTGCCGAAACATCGTCGATGGGAAAAATCATCTACCACACGATGTTCGATAACATGCAAAAAGCGAtcg                     <  1:1381529/69‑3 (MQ=255)
                         aTGCCGAAACATCGTCGATGGGAAAAATCATCTACCACACGATGTTCGATAACATGCAAAAAGCGAtcg                     <  1:457099/69‑3 (MQ=255)
                         aTGCCGAAACATCGTCGATGGGAAAAATCATCTACCACACGATGTTCGATAACATGCAAAAAGCGAtcg                     <  1:1892023/69‑3 (MQ=255)
                         aTGCCGAAACATCGTCGATGGGAAAAATCATCTACCACACGATGTTCGATAACATGCAAAAAGCGAtcg                     <  1:2531897/69‑3 (MQ=255)
                         aTGCCGAAACATCGTCGATGGGAAAAATCATCTACCACACGATGTTCGATAACATGCAAAAAGCGAtcg                     <  1:193624/69‑3 (MQ=255)
                             cGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGa                        <  1:2564279/62‑1 (MQ=255)
                                   gtcgtcGATGGGGATAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGa                        <  1:1649993/56‑1 (MQ=255)
                                    tcgtcgATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGt          <  1:13653/69‑1 (MQ=255)
                                    tcgtcgATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGt          <  1:42215/69‑1 (MQ=255)
                                    tcgtcgATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGt          <  1:3377644/69‑1 (MQ=255)
                                      gtcgATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcc        >  1:1523254/1‑69 (MQ=255)
                                       tcgATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGa                        <  1:1676378/52‑1 (MQ=255)
                                        cgATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcctg      <  1:1089296/69‑1 (MQ=255)
                                        cgATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcctg      <  1:575733/69‑1 (MQ=255)
                                        cgATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcctg      <  1:339634/69‑1 (MQ=255)
                                        cgATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcctg      <  1:2223423/69‑1 (MQ=255)
                                        cgATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcctg      <  1:329106/69‑1 (MQ=255)
                                        cgATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcctg      <  1:2323538/69‑1 (MQ=255)
                                        cgATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcct       <  1:3156365/68‑1 (MQ=255)
                                        cgATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcct       <  1:794749/68‑1 (MQ=255)
                                        cgATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcct       <  1:3028901/68‑1 (MQ=255)
                                        cgATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcct       <  1:2649686/68‑1 (MQ=255)
                                        cgATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcct       <  1:1823255/68‑1 (MQ=255)
                                        cgATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcct       <  1:38067/68‑1 (MQ=255)
                                         gATTGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcctg      <  1:464905/68‑1 (MQ=255)
                                         gATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcctg      <  1:1218020/68‑1 (MQ=255)
                                            ggggAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcctgcctg  <  1:186368/69‑1 (MQ=255)
                                            ggggAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcctgcctg  <  1:2901710/69‑1 (MQ=255)
                                            ggggAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcctgcctg  <  1:2205985/69‑1 (MQ=255)
                                            ggggAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcctgcctg  <  1:2195571/69‑1 (MQ=255)
                                            ggggAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcctgcctg  <  1:820131/69‑1 (MQ=255)
                                            ggggAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcctgcctg  <  1:920520/69‑1 (MQ=255)
                                            ggggAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTcctgcctg  <  1:3374923/69‑1 (MQ=255)
                                               |                                                                 
GTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGATTGACCTGAATCGTCCTGCCTG  >  minE/2057275‑2057387

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: