Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 370,232 G→A 38.3% E393E (GAG→GAA entF → enterobactin synthase multienzyme complex component, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE370,2320GA38.3% 40.0 / 45.0 73E393E (GAG→GAAentFenterobactin synthase multienzyme complex component, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (38/7);  new base A (23/5);  total (61/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TTAATGACCTTGAACTGGCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGAAAATG  >  minE/370177‑370298
                                                       |                                                                  
ttAATGACCTTGAACTGGCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATa                                                       <  1:285940/69‑1 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATaaa                                                     >  1:1204462/1‑59 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGt                                               >  1:1019910/1‑65 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:3332373/1‑68 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:2717090/1‑68 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:2909629/1‑68 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:2283278/1‑68 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:2197721/1‑68 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:2193751/1‑68 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:2189874/1‑68 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:2934442/1‑68 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:2995588/1‑68 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:1953514/1‑68 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:416198/1‑68 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:1579225/1‑68 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:487565/1‑68 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:543680/1‑68 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:862325/1‑68 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:1060568/1‑68 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:595170/1‑68 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:1522210/1‑68 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:1524177/1‑68 (MQ=255)
            aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:531243/1‑68 (MQ=255)
            aaCTGGACCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                                            >  1:4530/1‑68 (MQ=255)
             aCTGGCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTAc                                            <  1:2660115/67‑1 (MQ=255)
             aCTGGCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTAc                                            <  1:3065079/67‑1 (MQ=255)
             aCTGGCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTAc                                            <  1:2255740/67‑1 (MQ=255)
             aCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTACg                                           <  1:1007917/68‑1 (MQ=255)
             aCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTACg                                           <  1:1817471/68‑1 (MQ=255)
             aCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTACg                                           <  1:3198189/68‑1 (MQ=255)
             aCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTACg                                           <  1:2417327/68‑1 (MQ=255)
             aCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTACg                                           <  1:1365916/68‑1 (MQ=255)
                ggCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATg                                        >  1:1109412/1‑68 (MQ=255)
                 gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATg                                        >  1:1804682/1‑67 (MQ=255)
                 gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                                       >  1:73429/1‑68 (MQ=255)
                 gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                                       >  1:637007/1‑68 (MQ=255)
                 gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                                       >  1:602488/1‑68 (MQ=255)
                 gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                                       >  1:427149/1‑68 (MQ=255)
                 gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                                       >  1:565487/1‑68 (MQ=255)
                 gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                                       >  1:3244935/1‑68 (MQ=255)
                 gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                                       >  1:262478/1‑68 (MQ=255)
                 gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                                       >  1:2616922/1‑68 (MQ=255)
                 gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                                       >  1:1416495/1‑68 (MQ=255)
                 gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                                       >  1:1668498/1‑68 (MQ=255)
                 gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                                       >  1:2172252/1‑68 (MQ=255)
                 gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                                       >  1:1156995/1‑68 (MQ=255)
                 gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                                       >  1:2492057/1‑68 (MQ=255)
                 gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                                       >  1:2026858/1‑68 (MQ=255)
                 gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                                       >  1:1791278/1‑68 (MQ=255)
                 gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                                       >  1:2437281/1‑68 (MQ=255)
                      gTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATa                                                       >  1:156599/1‑47 (MQ=255)
                      gTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATa                                                       >  1:3043109/1‑47 (MQ=255)
                      gTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATa                                                       >  1:2880753/1‑47 (MQ=255)
                      gTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATa                                                       >  1:2788787/1‑47 (MQ=255)
                      gTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATa                                                       >  1:1726588/1‑47 (MQ=255)
                               tGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATcc                         >  1:1797272/1‑68 (MQ=255)
                               tGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATcc                         >  1:78516/1‑68 (MQ=255)
                               tGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATcc                         >  1:418560/1‑68 (MQ=255)
                               tGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATcc                         >  1:2492228/1‑68 (MQ=255)
                               tGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATcc                         >  1:2940374/1‑68 (MQ=255)
                               tGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATcc                         >  1:2582137/1‑68 (MQ=255)
                               tGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATcc                         >  1:2145911/1‑68 (MQ=255)
                                       gTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTg                <  1:2505991/69‑1 (MQ=255)
                                       gTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTg                <  1:1045032/69‑1 (MQ=255)
                                        tGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTg                <  1:3389511/68‑1 (MQ=255)
                                                   ttGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCACCGTTAATCCAg                       >  1:508837/1‑50 (MQ=255)
                                                   ttGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGaa      >  1:2499147/1‑67 (MQ=255)
                                                   ttGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGaa      >  1:2975915/1‑67 (MQ=255)
                                                   ttGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGaa      >  1:2629321/1‑67 (MQ=255)
                                                   ttGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGaa      >  1:1135837/1‑67 (MQ=255)
                                                   ttGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGaa      >  1:3259978/1‑67 (MQ=255)
                                                   ttGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGaa      >  1:1533986/1‑67 (MQ=255)
                                                   ttGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGaa      >  1:3426116/1‑67 (MQ=255)
                                                     gagaTCCTCTCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGc           <  1:1573009/60‑1 (MQ=255)
                                                     gagaTCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGc           <  1:2706657/60‑1 (MQ=255)
                                                     gagaTCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGc           <  1:3242526/60‑1 (MQ=255)
                                                     gagaTCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGc           <  1:436726/60‑1 (MQ=255)
                                                     gagaTCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGc           <  1:3014112/60‑1 (MQ=255)
                                                     gagaTCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGAAAAt   >  1:380967/1‑68 (MQ=255)
                                                     gagaTCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGAAAATg  >  1:863394/1‑69 (MQ=255)
                                                     gagaTCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGAAAATg  >  1:1479182/1‑69 (MQ=255)
                                                     gagaTCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGAAAATg  >  1:80316/1‑69 (MQ=255)
                                                     gagaTCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGAAAATg  >  1:1416090/1‑69 (MQ=255)
                                                     gagaTCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGAAAATg  >  1:489977/1‑69 (MQ=255)
                                                     gagaTCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGAAAATg  >  1:3346143/1‑69 (MQ=255)
                                                     gagaTCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGAAAATg  >  1:2061851/1‑69 (MQ=255)
                                                     gagaTCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGAAAATg  >  1:2159764/1‑69 (MQ=255)
                                                     gagaTCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGAAAATg  >  1:2316455/1‑69 (MQ=255)
                                                     gagaTCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGAAAATg  >  1:2908322/1‑69 (MQ=255)
                                                     gagaTCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGAAAATg  >  1:2862685/1‑69 (MQ=255)
                                                     gagaTCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGAAAATg  >  1:2852040/1‑69 (MQ=255)
                                                     gagaTCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGAAAATg  >  1:2460236/1‑69 (MQ=255)
                                                     gagaTCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGAAAATg  >  1:2548277/1‑69 (MQ=255)
                                                       |                                                                  
TTAATGACCTTGAACTGGCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCCAGCATGCTGAACGCCTGAAAATG  >  minE/370177‑370298

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: