Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 370,210 T→A 39.9% V386E (GTA→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,2100TA39.9% 40.9 / 63.4 70V386E (GTA→GAA) entFenterobactin synthase multienzyme complex component, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (38/4);  new base A (23/5);  total (61/9)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCCATACCCTGGCAACCGGTCCGGTTAATGACCTTGAACTGGCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCC  >  minE/370151‑370275
                                                           |                                                                 
aaCCCATACCCTGGCAACCGGTCCGGTTAATGACCTTGAACTGGCCCTGTTCCCGGATGTACACGGTGa                                                          >  1:2721858/1‑69 (MQ=255)
aaCCCATACCCTGGCAACCGGTCCGGTTAATGACCTTGAACTGGCCCTGTTCCCGGATGTACACGGTGa                                                          >  1:1620842/1‑69 (MQ=255)
aaCCCATACCCTGGCAACCGGTCCGGTTAATGACCTTGAACTGGCCCTGTTCCCGGATGTACACGGTGa                                                          >  1:1559413/1‑69 (MQ=255)
  cccATACCCTGGCAACCGGTCCGGTTAATGACCTTGAACTGGCCCTGTTCCCGGATGTACACGGTGa                                                          >  1:971245/1‑67 (MQ=255)
  cccATACCCTGGCAACCGGTCCGGTTAATGACCTTGAACTGGCCCTGTTCCCGGATGTACACGGTGa                                                          >  1:1074194/1‑67 (MQ=255)
  cccATACCCTGGCAACCGGTCCGGTTAATGACCTTGAACTGGCCCTGTTCCCGGATGTACACGGTGa                                                          >  1:3341343/1‑67 (MQ=255)
  cccATACCCTGGCAACCGGTCCGGTTAATGACCTTGAACTGGCCCTGTTCCCGGATGTACACGGTGa                                                          >  1:2633958/1‑67 (MQ=255)
        cccTGGCAACCGGTCCGGTTAATGACCTTGAACTGGCCCTGTTCCCGGATGTACACGGTGATTTGAGTa                                                  >  1:3390953/1‑69 (MQ=255)
                          ttAATGACCTTGAACTGGCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATa                                <  1:285940/69‑1 (MQ=255)
                                      aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATaaa                              >  1:1204462/1‑59 (MQ=255)
                                      aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGt                        >  1:1019910/1‑65 (MQ=255)
                                      aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                     >  1:416198/1‑68 (MQ=255)
                                      aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                     >  1:2995588/1‑68 (MQ=255)
                                      aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                     >  1:2934442/1‑68 (MQ=255)
                                      aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                     >  1:2909629/1‑68 (MQ=255)
                                      aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                     >  1:2283278/1‑68 (MQ=255)
                                      aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                     >  1:1524177/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:3332373/1‑68 (MQ=255)
                                      aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                     >  1:1060568/1‑68 (MQ=255)
                                      aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                     >  1:1953514/1‑68 (MQ=255)
                                      aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                     >  1:487565/1‑68 (MQ=255)
                                      aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                     >  1:531243/1‑68 (MQ=255)
                                      aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                     >  1:543680/1‑68 (MQ=255)
                                      aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                     >  1:595170/1‑68 (MQ=255)
                                      aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                     >  1:862325/1‑68 (MQ=255)
                                      aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                     >  1:1522210/1‑68 (MQ=255)
                                      aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                     >  1:1579225/1‑68 (MQ=255)
                                      aaCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTAc                     >  1:2717090/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:3198189/68‑1 (MQ=255)
                                       aCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTACg                    <  1:1007917/68‑1 (MQ=255)
                                       aCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTACg                    <  1:1365916/68‑1 (MQ=255)
                                       aCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTACg                    <  1:1817471/68‑1 (MQ=255)
                                       aCTGGCCCTGTTCCCGGATGAACACGGTGATTTGAGTATTGAAATCCTCGCCAATAAACAGCGTTACg                    <  1:2417327/68‑1 (MQ=255)
                                          ggCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATg                 >  1:1109412/1‑68 (MQ=255)
                                           gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATg                 >  1:1804682/1‑67 (MQ=255)
                                           gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                >  1:2437281/1‑68 (MQ=255)
                                           gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                >  1:427149/1‑68 (MQ=255)
                                           gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                >  1:1791278/1‑68 (MQ=255)
                                           gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                >  1:2026858/1‑68 (MQ=255)
                                           gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                >  1:565487/1‑68 (MQ=255)
                                           gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                >  1:2172252/1‑68 (MQ=255)
                                           gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                >  1:3244935/1‑68 (MQ=255)
                                           gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                >  1:2492057/1‑68 (MQ=255)
                                           gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                >  1:602488/1‑68 (MQ=255)
                                           gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                >  1:637007/1‑68 (MQ=255)
                                           gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                >  1:73429/1‑68 (MQ=255)
                                           gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                >  1:1668498/1‑68 (MQ=255)
                                           gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                >  1:1156995/1‑68 (MQ=255)
                                           gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                >  1:1416495/1‑68 (MQ=255)
                                           gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                >  1:2616922/1‑68 (MQ=255)
                                           gCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGa                >  1:262478/1‑68 (MQ=255)
                                                gTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATa                                >  1:156599/1‑47 (MQ=255)
                                                gTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATa                                >  1:1726588/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)
                                                         tGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATcc  >  1:1797272/1‑68 (MQ=255)
                                                         tGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATcc  >  1:418560/1‑68 (MQ=255)
                                                         tGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATcc  >  1:2145911/1‑68 (MQ=255)
                                                         tGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATcc  >  1:2940374/1‑68 (MQ=255)
                                                         tGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATcc  >  1:2492228/1‑68 (MQ=255)
                                                         tGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATcc  >  1:2582137/1‑68 (MQ=255)
                                                         tGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATcc  >  1:78516/1‑68 (MQ=255)
                                                           |                                                                 
AACCCATACCCTGGCAACCGGTCCGGTTAATGACCTTGAACTGGCCCTGTTCCCGGATGTACACGGTGATTTGAGTATTGAGATCCTCGCCAATAAACAGCGTTACGATGAGCCAACGTTAATCC  >  minE/370151‑370275

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: