Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,301,363 G→A 47.0% Y236Y (TAC→TAT 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,3630GA47.0% 2.6 / 41.1 51Y236Y (TAC→TATmetLfused aspartokinase II and homoserine dehydrogenase II
Reads supporting (aligned to +/- strand):  ref base G (18/9);  new base A (14/10);  total (32/19)
Fisher's exact test for biased strand distribution p-value = 5.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCTGGCCTCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTCG  >  minE/2301295‑2301430
                                                                    |                                                                   
tCGCTGGCCTCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                     <  1:2030712/69‑1 (MQ=255)
tCGCTGGCCTCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                     <  1:1575875/69‑1 (MQ=255)
tCGCTGGCCTCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                     <  1:131666/69‑1 (MQ=255)
tCGCTGGCCTCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                     <  1:3245799/69‑1 (MQ=255)
   cTGGCCTCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGta                                                                   <  1:618804/68‑1 (MQ=255)
        ttCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCGCTATATAcccc                                                             >  1:1120486/2‑69 (MQ=255)
        ttCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCGCTATATAcccc                                                             >  1:28218/2‑69 (MQ=255)
        ttCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCGCTATATAcccc                                                             >  1:3055827/2‑69 (MQ=255)
        ttCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCGCTATATAcccc                                                             >  1:1790432/2‑69 (MQ=255)
        ttCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCGCTATATAcccc                                                             >  1:171681/2‑69 (MQ=255)
             ccAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                     <  1:276454/56‑1 (MQ=255)
             ccAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                     <  1:678141/56‑1 (MQ=255)
             ccAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                     <  1:3371694/56‑1 (MQ=255)
             ccAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                     <  1:1883804/56‑1 (MQ=255)
             ccAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                     <  1:2828466/56‑1 (MQ=255)
             ccAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                     <  1:2662726/56‑1 (MQ=255)
             ccAGACGCAGCAACGGCAGAAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                     <  1:3291266/56‑1 (MQ=255)
               aGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                     <  1:3066877/54‑1 (MQ=255)
                       caACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                     <  1:880994/46‑1 (MQ=255)
                        aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                             >  1:1544740/1‑69 (MQ=255)
                        aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                             >  1:905144/1‑69 (MQ=255)
                        aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                             >  1:1409557/1‑69 (MQ=255)
                        aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                             >  1:2897768/1‑69 (MQ=255)
                        aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                             >  1:560808/1‑69 (MQ=255)
                        aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                             >  1:1322447/1‑69 (MQ=255)
                        aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                             >  1:1310821/1‑69 (MQ=255)
                        aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                             >  1:1946340/1‑69 (MQ=255)
                        aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                             >  1:1284631/1‑69 (MQ=255)
                                         aTCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCgaga                               <  1:1474981/66‑1 (MQ=255)
                                         aTCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCgaga                               <  1:2560977/66‑1 (MQ=255)
                                         aTCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCgaga                               <  1:2433075/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:368714/66‑1 (MQ=255)
                                            tttCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACAcc                         >  1:1249806/1‑69 (MQ=255)
                                            tttCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACAcc                         >  1:2671538/1‑69 (MQ=255)
                                                        cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc             >  1:2862053/1‑69 (MQ=255)
                                                        cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc             >  1:316225/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:2293764/1‑69 (MQ=255)
                                                        cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc             >  1:1806371/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:960774/69‑1 (MQ=255)
                                                             cGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt        <  1:471174/69‑1 (MQ=255)
                                                             cGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt        <  1:555310/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:1790112/1‑69 (MQ=255)
                                                              ggCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTtg       >  1:1886709/1‑69 (MQ=255)
                                                              ggCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTtg       >  1:2311460/1‑69 (MQ=255)
                                                              ggCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCCCAGCGCACCGATTtgt      >  1:773366/1‑69 (MQ=255)
                                                                    gTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTcg  >  1:187046/1‑68 (MQ=255)
                                                                    gTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTcg  >  1:349446/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)
                                                                    |                                                                   
TCGCTGGCCTCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGTGTCG  >  minE/2301295‑2301430

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: