Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,703,017 T→C 70.0% D136D (GAT→GAC aroB → 3‑dehydroquinate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,703,0170TC70.0% 71.5 / 31.6 73D136D (GAT→GACaroB3‑dehydroquinate synthase
Reads supporting (aligned to +/- strand):  ref base T (7/15);  new base C (44/7);  total (51/22)
Fisher's exact test for biased strand distribution p-value = 8.25e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

CGGCGGCGAGTTATCAGCGCGGTGTCCGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGATTCCTCCGTTGGCGGCAAAACTGCGGTCAACCATCCCCTCGGTAAAAACATGATTGGCGCGTTCT  >  minE/2702950‑2703081
                                                                   |                                                                 
cggcggcgAGTTATCAGCGCGGTGTCCGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGAt                                                                   <  1:837856/68‑1 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGTTCGActcctc                                                              >  1:1398791/1‑62 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACtcctcc                                                             >  1:1963269/1‑63 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACtcctcc                                                             >  1:2657887/1‑63 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACtcctcc                                                             >  1:2614778/1‑63 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACtcctcc                                                             >  1:2592158/1‑63 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACtcctcc                                                             >  1:2830621/1‑63 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACtcctcc                                                             >  1:2869760/1‑63 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACtcctcc                                                             >  1:2424265/1‑63 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACtcctcc                                                             >  1:3182739/1‑63 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACtcctcc                                                             >  1:3248296/1‑63 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACtcctcc                                                             >  1:424673/1‑63 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACtcctcc                                                             >  1:1971842/1‑63 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACtcctcc                                                             >  1:1908209/1‑63 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACtcctcc                                                             >  1:991748/1‑63 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACtcctcc                                                             >  1:564133/1‑63 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACtcctcc                                                             >  1:904084/1‑63 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACtcctcc                                                             >  1:127697/1‑63 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACtcctcc                                                             >  1:891043/1‑63 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACtcctcc                                                             >  1:1359903/1‑63 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACtcctcc                                                             >  1:1368784/1‑63 (MQ=255)
           tATCAGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACtcctcc                                                             >  1:1806597/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:3276695/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:1083368/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:2896988/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:856908/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:800176/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:3379288/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:402361/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:75599/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:427111/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:506923/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:649999/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:2737730/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:561572/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:1443449/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:2518364/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:2481261/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:1723283/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:1780599/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:1858587/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:200453/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:2136352/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:2188028/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:2401122/1‑63 (MQ=255)
               aGCGCGGTGTACGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTg                                                         >  1:2191849/1‑63 (MQ=255)
                       gtCCGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGATTCCTCCGTTGGCGGCAAAACTGCg                                          <  1:1299302/69‑1 (MQ=255)
                       gtCCGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGATTCCTCCGTTGGCGGCAAAACTGCg                                          >  1:2753968/1‑69 (MQ=255)
                       gtCCGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGATTCCTCCGTTGGCGGCAAAACTGCg                                          <  1:809807/69‑1 (MQ=255)
                        tCCGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGATTCCTCCGTTGGCGGCAAAACTGCg                                          >  1:110386/1‑68 (MQ=255)
                        tCCGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGATTCCTCCGTTGGCGGCAAAACTGCg                                          <  1:1646648/68‑1 (MQ=255)
                            tttCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGATTCCTCCGTTGGCGGCAAAACTGCGGTCa                                      <  1:2851223/68‑1 (MQ=255)
                            tttCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGATTCCTCCGTTGGCGGCAAAACTGCGGTCa                                      <  1:466942/68‑1 (MQ=255)
                            tttCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGATTCCTCCGTTGGCGGCAAAACTGCGGTCa                                      <  1:512085/68‑1 (MQ=255)
                             ttcattcaAGTCCCGACGACGTTACTGTCGCAGGTCGATTCCTCCGTTGGCGGCAAAACTGCGGTCAAc                                    <  1:1189717/69‑1 (MQ=255)
                             ttcattcaAGTCCCGACGACGTTACTGTCGCAGGTCGATTCCTCCGTTGGCGGCAAAACTGCGGTCAAc                                    <  1:1501460/69‑1 (MQ=255)
                             ttcattcaAGTCCCGACGACGTTACTGTCGCAGGTCGATTCCTCCGTTGGCGGCAAAACTGCGGTCAAc                                    <  1:3230579/69‑1 (MQ=255)
                             ttcattcaAGTCCCGACGACGTTACTGTCGCAGGTCGATTCCTCCGTTGGCGGCAAAACTGCGGTCAAc                                    <  1:2113349/69‑1 (MQ=255)
                              tcattcaAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTGGCGGCAAAACTGCGGTCAAcc                                   <  1:2405737/69‑1 (MQ=255)
                              tcattcaAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTGGCGGCAAAACTGCGGTCAAcc                                   <  1:2625493/69‑1 (MQ=255)
                              tcattcaAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTGGCGGCAAAACTGCGGTCAAcc                                   <  1:2346006/69‑1 (MQ=255)
                              tcattcaAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTGGCGGCAAAACTGCGGTCAAcc                                   <  1:3209987/69‑1 (MQ=255)
                              tcattcaAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTGGCGGCAAAACTGCGGTCAAcc                                   <  1:485455/69‑1 (MQ=255)
                              tcattcaAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTGGCGGCAAAACTGCGGTCAAcc                                   <  1:3229149/69‑1 (MQ=255)
                              tcattcaAGTCCCGACGACGTTACTGTCGCAGGTCGACTCCTCCGTTGGCGGCAAAACTGCGGTCAAcc                                   <  1:345789/69‑1 (MQ=255)
                               cattcaAGTCCCGACGACGTTACTGTCGCAGGTCGATTCCTCCGTTGGCGGCAAAAACTGCGGTCAAc                                    <  1:2704106/68‑1 (MQ=255)
                                                cgTTACTGTCGCAGGTCGATTCCTCCGTTGGCGGCAAAACTGCGGTCAAc                                    <  1:947373/50‑1 (MQ=255)
                                                  ttACTGTCGCAGGTCGATTCCTCCGTTGGCGGCAAAACTGCGGTCAACCATCCCCTCGGTAAAAACATg               <  1:2548445/69‑1 (MQ=255)
                                                   tACTGTCGCAGGTCGATTCCTCCGTTGGCGGCAAAACTGCGGTCAACCATCCCCTCGGTAAAAACATg               <  1:577078/68‑1 (MQ=255)
                                                            aGGTCGATTCCTCCGTTGGCGGCAAAACTGCGGTCa                                      <  1:2975427/36‑1 (MQ=255)
                                                              gTCGATTCCTCCGTTGGCGGCAAAACTGCGGTCAACCATCCCCTCGGTAAAAACATGATTGGCGCGTTc   >  1:2657442/1‑69 (MQ=255)
                                                               tCGATTCCTCCGTTGGCGGCAAAACTGCGGTCAACCATCCCCTCGGTAAAAAc                  >  1:1538983/1‑53 (MQ=255)
                                                                cGATTCCTCCGTTGGCGGCAAAACTGCGGTCAACCATCCCCTCGGTAAAAACATGATTGGCGCGTTCt  >  1:3094384/1‑68 (MQ=255)
                                                                cGATTCCTCCGTTGGCGGCAAAACTGCGGTCAACCATCCCCTCGGTAAAAACATGATTGGCGCGTTCt  >  1:3254970/1‑68 (MQ=255)
                                                                cGATTCCTCCGTTGGCGGCAAAACTGCGGTCAACCATCCCCTCGGTAAAAACATGATTGGCGCGTTCt  >  1:3345252/1‑68 (MQ=255)
                                                                   |                                                                 
CGGCGGCGAGTTATCAGCGCGGTGTCCGTTTCATTCAAGTCCCGACGACGTTACTGTCGCAGGTCGATTCCTCCGTTGGCGGCAAAACTGCGGTCAACCATCCCCTCGGTAAAAACATGATTGGCGCGTTCT  >  minE/2702950‑2703081

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: