Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,301,360 A→G 37.5% S237S (AGT→AGC 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,3600AG37.5% 31.6 / 40.0 64S237S (AGT→AGCmetLfused aspartokinase II and homoserine dehydrogenase II
Reads supporting (aligned to +/- strand):  ref base A (18/22);  new base G (14/10);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

TCGCTGGCCTCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGT  >  minE/2301295‑2301426
                                                                 |                                                                  
tCGCTGGCCTCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                 <  1:3245799/69‑1 (MQ=255)
tCGCTGGCCTCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                 <  1:2030712/69‑1 (MQ=255)
tCGCTGGCCTCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                 <  1:1575875/69‑1 (MQ=255)
tCGCTGGCCTCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                 <  1:131666/69‑1 (MQ=255)
   cTGGCCTCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGta                                                               <  1:618804/68‑1 (MQ=255)
        ttCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCGCTATATAcccc                                                         >  1:1790432/2‑69 (MQ=255)
        ttCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCGCTATATAcccc                                                         >  1:1120486/2‑69 (MQ=255)
        ttCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCGCTATATAcccc                                                         >  1:171681/2‑69 (MQ=255)
        ttCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCGCTATATAcccc                                                         >  1:3055827/2‑69 (MQ=255)
        ttCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCGCTATATAcccc                                                         >  1:28218/2‑69 (MQ=255)
             ccAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                 <  1:2662726/56‑1 (MQ=255)
             ccAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                 <  1:2828466/56‑1 (MQ=255)
             ccAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                 <  1:3371694/56‑1 (MQ=255)
             ccAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                 <  1:276454/56‑1 (MQ=255)
             ccAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                 <  1:1883804/56‑1 (MQ=255)
             ccAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTg                                                                 <  1:678141/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:1409557/1‑69 (MQ=255)
                        aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                         >  1:1946340/1‑69 (MQ=255)
                        aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                         >  1:2897768/1‑69 (MQ=255)
                        aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                         >  1:1544740/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:1284631/1‑69 (MQ=255)
                        aaCGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAg                                         >  1:905144/1‑69 (MQ=255)
                                         aTCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCgaga                           <  1:3316026/66‑1 (MQ=255)
                                         aTCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCgaga                           <  1:2560977/66‑1 (MQ=255)
                                         aTCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCgaga                           <  1:2433075/66‑1 (MQ=255)
                                         aTCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCgaga                           <  1:2056318/66‑1 (MQ=255)
                                         aTCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCgaga                           <  1:368714/66‑1 (MQ=255)
                                         aTCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCgaga                           <  1:114602/66‑1 (MQ=255)
                                         aTCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCgaga                           <  1:1474981/66‑1 (MQ=255)
                                            tttCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACAcc                     >  1:2671538/1‑69 (MQ=255)
                                            tttCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACAcc                     >  1:1249806/1‑69 (MQ=255)
                                                        cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc         >  1:2403478/1‑69 (MQ=255)
                                                        cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc         >  1:1037018/1‑69 (MQ=255)
                                                        cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc         >  1:1081267/1‑69 (MQ=255)
                                                        cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc         >  1:316225/1‑69 (MQ=255)
                                                        cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc         >  1:2862053/1‑69 (MQ=255)
                                                        cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc         >  1:1806371/1‑69 (MQ=255)
                                                        cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc         >  1:1876609/1‑69 (MQ=255)
                                                        cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc         >  1:2537776/1‑69 (MQ=255)
                                                        cgGGTCGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCAcc         >  1:2293764/1‑69 (MQ=255)
                                                             cGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt    <  1:555310/69‑1 (MQ=255)
                                                             cGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt    <  1:471174/69‑1 (MQ=255)
                                                             cGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt    <  1:363356/69‑1 (MQ=255)
                                                             cGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt    <  1:960774/69‑1 (MQ=255)
                                                             cGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt    <  1:1395086/69‑1 (MQ=255)
                                                             cGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt    <  1:1357382/69‑1 (MQ=255)
                                                             cGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt    <  1:1139438/69‑1 (MQ=255)
                                                             cGGCGCTATATACCCCGGCGACGTCACTCCAAATGGTTACGCGAGAAACACCCGCCAGCGCACCGAttt    <  1:1132467/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:2311460/1‑69 (MQ=255)
                                                              ggCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTtg   >  1:1886709/1‑69 (MQ=255)
                                                              ggCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTtg   >  1:1790112/1‑69 (MQ=255)
                                                              ggCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTtg   >  1:613260/1‑69 (MQ=255)
                                                              ggCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTtg   >  1:1062728/1‑69 (MQ=255)
                                                              ggCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCCCAGCGCACCGATTtgt  >  1:773366/1‑69 (MQ=255)
                                                                 |                                                                  
TCGCTGGCCTCATCCAGACGCAGCAACGGCAGCAGGCAGGCATCTTTCACTTTACGCGGGTCGGCACTGTATACCCCGGCGACGTCGCTCCAGATGGTTACGCGAGAAACACCCGCCAGCGCACCGATTTGT  >  minE/2301295‑2301426

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: