Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,039,489 A→C 60.6% A142A (GCT→GCG ribB ← 3,4 dihydroxy‑2‑butanone‑4‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,039,4890AC60.6% 71.2 / 71.0 91A142A (GCT→GCGribB3,4 dihydroxy‑2‑butanone‑4‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (27/9);  new base C (33/22);  total (60/31)
Fisher's exact test for biased strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCAATAGTTGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAGATCTGACGGTTTTGCGCCATCGGCAATCGC  >  minE/2039437‑2039554
                                                    |                                                                 
gATCGATAGTCGCTTCAGTATGCCCGCCACGCGTCAGTACACCACCTGCCTgcgca                                                                >  1:644066/1‑55 (MQ=255)
gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGa                                                             >  1:3205240/1‑59 (MQ=255)
gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGt                                                          >  1:512270/1‑62 (MQ=255)
gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTg                                                         >  1:2357946/1‑63 (MQ=255)
gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   >  1:3051977/1‑69 (MQ=255)
gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   >  1:470323/1‑69 (MQ=255)
gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   >  1:2114579/1‑69 (MQ=255)
gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   >  1:522846/1‑69 (MQ=255)
gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   >  1:546245/1‑69 (MQ=255)
gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   >  1:3217824/1‑69 (MQ=255)
gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   >  1:2482290/1‑69 (MQ=255)
gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   >  1:1275433/1‑69 (MQ=255)
gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   >  1:697736/1‑69 (MQ=255)
gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   >  1:3166001/1‑69 (MQ=255)
gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   >  1:1259955/1‑69 (MQ=255)
gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   >  1:3123109/1‑69 (MQ=255)
gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAACTGGGaaaa                                                   >  1:3054708/1‑69 (MQ=255)
gATCGAAAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   >  1:435884/1‑69 (MQ=255)
 aTCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   <  1:2628769/68‑1 (MQ=255)
 aTCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   <  1:3327770/68‑1 (MQ=255)
 aTCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   <  1:2176274/68‑1 (MQ=255)
 aTCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   <  1:3015684/68‑1 (MQ=255)
 aTCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   <  1:2037603/68‑1 (MQ=255)
 aTCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   <  1:1084050/68‑1 (MQ=255)
 aTCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   <  1:44283/68‑1 (MQ=255)
 aTCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   <  1:466737/68‑1 (MQ=255)
 aTCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   <  1:1112096/68‑1 (MQ=255)
 aTCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa                                                   <  1:501995/68‑1 (MQ=255)
   cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACg                                                 >  1:284666/3‑68 (MQ=255)
   cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt                                                >  1:657071/3‑69 (MQ=255)
   cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt                                                >  1:1073749/3‑69 (MQ=255)
   cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt                                                >  1:1469614/3‑69 (MQ=255)
   cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt                                                >  1:3230613/3‑69 (MQ=255)
   cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt                                                >  1:531306/3‑69 (MQ=255)
   cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt                                                >  1:1882218/3‑69 (MQ=255)
   cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt                                                >  1:1326178/3‑69 (MQ=255)
   cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt                                                >  1:3361660/3‑69 (MQ=255)
   cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt                                                >  1:2336842/3‑69 (MQ=255)
   cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt                                                >  1:2277825/3‑69 (MQ=255)
   cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt                                                >  1:224411/3‑69 (MQ=255)
   cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt                                                >  1:2214028/3‑69 (MQ=255)
   cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt                                                >  1:1646717/3‑69 (MQ=255)
   cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt                                                >  1:1686953/3‑69 (MQ=255)
   cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt                                                >  1:1823740/3‑69 (MQ=255)
        gTTGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGcc                                            >  1:2836692/1‑68 (MQ=255)
               cAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGATAACGTGGCCAGGACGa                                     >  1:73569/1‑68 (MQ=255)
               cAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGa                                     >  1:2796655/1‑68 (MQ=255)
                   aTGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGTCGATTCAg                                <  1:3332930/69‑1 (MQ=255)
                   aTGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAg                                <  1:1633184/69‑1 (MQ=255)
                   aTGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAg                                <  1:2296216/69‑1 (MQ=255)
                   aTGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAg                                <  1:2212691/69‑1 (MQ=255)
                   aTGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAg                                <  1:2128639/69‑1 (MQ=255)
                   aTGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAg                                <  1:2838746/69‑1 (MQ=255)
                   aTGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAg                                <  1:2029147/69‑1 (MQ=255)
                   aTGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAg                                <  1:1247840/69‑1 (MQ=255)
                   aTGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAg                                <  1:1946475/69‑1 (MQ=255)
                    tGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACg                                                 >  1:695006/1‑51 (MQ=255)
                    tGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACg                                                 >  1:757953/1‑51 (MQ=255)
                    tGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACg                                                 >  1:2427689/1‑51 (MQ=255)
                    tGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACg                                                 >  1:1620142/1‑51 (MQ=255)
                    tGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACg                                                 >  1:2423645/1‑51 (MQ=255)
                    tGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACg                                                 >  1:2398858/1‑51 (MQ=255)
                    tGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAGa                               >  1:2993775/1‑69 (MQ=255)
                    tGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAGa                               >  1:1630962/1‑69 (MQ=255)
                    tGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAGa                               >  1:1731172/1‑69 (MQ=255)
                    tGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAGa                               >  1:3354820/1‑69 (MQ=255)
                    tGACCGCCACGCGTCAGTACACCACCTGCCTGAGAGCGAAGTGGGAAAACGTGGCCAGGACGATTCAGa                               >  1:1927757/1‑69 (MQ=255)
                      aCCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAGAt                              >  1:1867957/1‑68 (MQ=255)
                                         ccaccTGCCTGAGCGCGAAGTGGGAAAACGTGGCTAGGACGATTCAGATCTGACGGTTTTGCGCCATc           >  1:2948814/1‑68 (MQ=255)
                                         ccaccTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAGATCTGACGGTTTTGCGCCATc           >  1:1172480/1‑68 (MQ=255)
                                         ccaccTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAGATCTGACGGTTTTGCGCCATc           >  1:445393/1‑68 (MQ=255)
                                         ccaccTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAGATCTGACGGTTTTGCGCCATc           >  1:3258596/1‑68 (MQ=255)
                                         ccaccTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAGATCTGACGGTTTTGCGCCATc           >  1:1364981/1‑68 (MQ=255)
                                         ccaccTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAGATCTGACGGTTTTGCGCCATc           >  1:1368858/1‑68 (MQ=255)
                                         ccaccTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAGATCTGACGGTTTTGCGCCATc           >  1:2337069/1‑68 (MQ=255)
                                         ccaccTGCCTGAGCGCGAAGTGGGAAAACGTGGCAAGGACGATTCAGATCTGACGGTTTTGCGCCATc           >  1:1692217/1‑68 (MQ=255)
                                                 ctgCGCACGAAGTGGGAAAACGTGGCCCGTACGATTCAg                                <  1:575920/39‑1 (MQ=21)
                                                 ctgCGCACGAAGTGGGAAAACGTGGCCCGGACGATTCat                                <  1:85148/39‑2 (MQ=25)
                                                 ctgCGCACGAAGTGGGAAAACGTGGCCCGGACGATTCAg                                <  1:528875/39‑1 (MQ=37)
                                                 ctgCGCACGAAGTGGGAAAACGTGGCCCGGACGATTCAg                                <  1:1270260/39‑1 (MQ=37)
                                                 ctgCGCACGAAGTGGGAAAACGTGGCCCGGACGATTCAg                                <  1:1272388/39‑1 (MQ=37)
                                                 ctgCGCACGAAGTGGGAAAACGTGGCCCGGACGATTCAg                                <  1:1282311/39‑1 (MQ=37)
                                                 ctgCGCACGAAGTGGGAAAACGTGGCCCGGACGATTCAg                                <  1:1481927/39‑1 (MQ=37)
                                                 ctgCGCACGAAGTGGGAAAACGTGGCCCGGACGATTCAg                                <  1:2646450/39‑1 (MQ=37)
                                                 ctgCGCACGAAGTGGGAAAACGTGGCCCGGACGATTCAg                                <  1:2427070/39‑1 (MQ=37)
                                                 ctgCGCACGAAGTGGGAAAACGTGGCCCGGACGATTCAg                                <  1:1853939/39‑1 (MQ=37)
                                                 ctgCGCACGAAGTGGGAAAACGTGGCCCGGACGATTCAg                                <  1:1911636/39‑1 (MQ=37)
                                                 ctgCGCACGAAGTGGGAAAACGTGGCCCGGACGATTCAg                                <  1:3383763/39‑1 (MQ=37)
                                                 ctgAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAGATCTGACGGTTTTGCGCCATCGGCAATCGc  >  1:2694111/1‑69 (MQ=255)
                                                 ctgAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAGATCTGACGGTTTTGCGCCATCGGCAATCGc  >  1:73455/1‑69 (MQ=255)
                                                 ctgAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAGATCTGACGGTTTTGCGCCATCGGCAATCGc  >  1:2404874/1‑69 (MQ=255)
                                                 ctgAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAGATCTGACGGTTTTGCGCCATCGGCAATCGc  >  1:38024/1‑69 (MQ=255)
                                                    |                                                                 
GATCAATAGTTGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGCCAGGACGATTCAGATCTGACGGTTTTGCGCCATCGGCAATCGC  >  minE/2039437‑2039554

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: