Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,301,381 G→A 37.8% S230S (AGC→AGT metL ← fused aspartokinase II and homoserine dehydrogenase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,301,3810GA37.8% 36.5 / 39.9 66S230S (AGC→AGTmetLfused aspartokinase II and homoserine dehydrogenase II
Reads supporting (aligned to +/- strand):  ref base G (30/11);  new base A (14/11);  total (44/22)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTCGCGGAATAGTCGGAACCGT  >  minE/2301319‑2301448
                                                              |                                                                   
aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                                               >  1:1409557/1‑69 (MQ=255)
aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                                               >  1:905144/1‑69 (MQ=255)
aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                                               >  1:560808/1‑69 (MQ=255)
aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                                               >  1:2897768/1‑69 (MQ=255)
aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                                               >  1:1946340/1‑69 (MQ=255)
aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                                               >  1:1322447/1‑69 (MQ=255)
aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                                               >  1:1544740/1‑69 (MQ=255)
aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                                               >  1:1284631/1‑69 (MQ=255)
aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                                               >  1:1310821/1‑69 (MQ=255)
                 aTCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCgaga                                                 <  1:368714/66‑1 (MQ=255)
                 aTCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCgaga                                                 <  1:1474981/66‑1 (MQ=255)
                 aTCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCgaga                                                 <  1:3316026/66‑1 (MQ=255)
                 aTCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCgaga                                                 <  1:114602/66‑1 (MQ=255)
                 aTCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCgaga                                                 <  1:2056318/66‑1 (MQ=255)
                 aTCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCgaga                                                 <  1:2433075/66‑1 (MQ=255)
                 aTCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCgaga                                                 <  1:2560977/66‑1 (MQ=255)
                    tttCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACAcc                                           >  1:2671538/1‑69 (MQ=255)
                    tttCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACAcc                                           >  1:1249806/1‑69 (MQ=255)
                                cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc                               >  1:2293764/1‑69 (MQ=255)
                                cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc                               >  1:1876609/1‑69 (MQ=255)
                                cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc                               >  1:2537776/1‑69 (MQ=255)
                                cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc                               >  1:2403478/1‑69 (MQ=255)
                                cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc                               >  1:2862053/1‑69 (MQ=255)
                                cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc                               >  1:1806371/1‑69 (MQ=255)
                                cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc                               >  1:316225/1‑69 (MQ=255)
                                cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc                               >  1:1081267/1‑69 (MQ=255)
                                cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc                               >  1:1037018/1‑69 (MQ=255)
                                     cGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt                          <  1:363356/69‑1 (MQ=255)
                                     cGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt                          <  1:471174/69‑1 (MQ=255)
                                     cGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt                          <  1:555310/69‑1 (MQ=255)
                                     cGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt                          <  1:960774/69‑1 (MQ=255)
                                     cGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt                          <  1:1132467/69‑1 (MQ=255)
                                     cGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt                          <  1:1139438/69‑1 (MQ=255)
                                     cGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt                          <  1:1357382/69‑1 (MQ=255)
                                     cGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt                          <  1:1395086/69‑1 (MQ=255)
                                     cGGCGCTATATACCCCGGCGACGCCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt                          <  1:1820038/69‑1 (MQ=255)
                                      ggCGCTATATACCCCGGCGACGTAACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt                          <  1:2964449/68‑1 (MQ=255)
                                      ggCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCgaga                                                 <  1:3430247/45‑1 (MQ=255)
                                      ggCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt                          >  1:1247251/1‑68 (MQ=255)
                                      ggCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTtg                         >  1:1062728/1‑69 (MQ=255)
                                      ggCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTtg                         >  1:613260/1‑69 (MQ=255)
                                      ggCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTtg                         >  1:1886709/1‑69 (MQ=255)
                                      ggCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTtg                         >  1:2311460/1‑69 (MQ=255)
                                      ggCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTtg                         >  1:1790112/1‑69 (MQ=255)
                                      ggCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCCCAGCGCACCGATTtgt                        >  1:773366/1‑69 (MQ=255)
                                            gTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTcg                    >  1:349446/1‑68 (MQ=255)
                                            gTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTcg                    >  1:187046/1‑68 (MQ=255)
                                            gTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTcg                    >  1:2127582/1‑68 (MQ=255)
                                            gTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTcg                    >  1:2654841/1‑68 (MQ=255)
                                            gTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTcg                    >  1:2751045/1‑68 (MQ=255)
                                              ataCCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAaca                                             >  1:601731/1‑41 (MQ=255)
                                              ataCCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTCGCgg                 >  1:405599/1‑69 (MQ=255)
                                              ataCCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTCGCgg                 >  1:811752/1‑69 (MQ=255)
                                              ataCCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTCGCgg                 >  1:1102316/1‑69 (MQ=255)
                                              ataCCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTCGCgg                 >  1:93398/1‑69 (MQ=255)
                                                        gACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTCGCGGAATAGTCGGa       >  1:2086835/1‑69 (MQ=255)
                                                        gACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTCGCGGAATAGTCGGa       >  1:2469621/1‑69 (MQ=255)
                                                        gACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTCGCGGAATAGTCGGa       >  1:2452164/1‑69 (MQ=255)
                                                        gACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTCGCGGAATAGTCGGa       >  1:88065/1‑69 (MQ=255)
                                                        gACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTCGCGGAATAGTCGGa       >  1:3217258/1‑69 (MQ=255)
                                                         aCGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt                          <  1:1133108/49‑1 (MQ=39)
                                                            tCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTCGCGGAATAGTCGGAAcc    >  1:212154/1‑68 (MQ=255)
                                                            tCGCTCCAGATGGTTACGCGAGAAACACCCGACAGCGCACCGATTTGTGTCGCGGAATAGTCGGAAcc    >  1:792087/1‑68 (MQ=255)
                                                             cGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTCGCGGAATAGTCGGAACCGt  <  1:1920577/69‑1 (MQ=255)
                                                             cGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTCGCGGAATAGTCGGAACCGt  <  1:374549/69‑1 (MQ=255)
                                                             cGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTCGCGGAATAGTCGGAACCGt  <  1:41085/69‑1 (MQ=255)
                                                              gCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTCGCGGAATAGTCGGAACCGt  <  1:271457/68‑1 (MQ=255)
                                                              |                                                                   
AACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTCGCGGAATAGTCGGAACCGT  >  minE/2301319‑2301448

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: