Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 364,740 G→C 52.5% V82V (GTC→GTG entD ← phosphopantetheinyltransferase component of enterobactin synthase multienzyme complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE364,7400GC52.5% 0.4 / 45.9 40V82V (GTC→GTGentDphosphopantetheinyltransferase component of enterobactin synthase multienzyme complex
Reads supporting (aligned to +/- strand):  ref base G (5/14);  new base C (19/2);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 5.69e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

GTAGTCCCACAGTGGCTAATACTGCCGTATACCTCCGCAGGCCAGACAGGTTGGCGTAGCTCGCCGATTGCGGGCACACATTTATAGCCATATTCCCGCAAAGCATA  >  minE/364696‑364802
                                            |                                                              
gTAGTCCCACAGTGGCTAATACTGCCGTATACCTCCGCAGGCCAGACAGGTTGGCGTAGCTCGCCGAtt                                        >  1:868267/1‑69 (MQ=255)
   gTCCCACAGTGACTAATGCTGCCGTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGCCGATTGca                                     <  1:842280/69‑2 (MQ=255)
   gTCCCACAGTGACTAATGCTGCCGTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGCCGATTGca                                     <  1:2747059/69‑2 (MQ=255)
             ggCTAATACTGCCGTATACCTCCGCAGGCCAGACAGGTTGGCGTAGCTCGCCGATTGCGGGCACACAtt                           <  1:3050386/69‑1 (MQ=255)
             ggCTAATACTGCCGTATACCTCCGCAGGCCAGACAGGTTGGCGTAGCTCGCCGATTGCGGGCACACAtt                           >  1:3345122/1‑69 (MQ=255)
              gCTAATACTGCCGTATACCTCCGCAGGCCAGACAGGTTGGCGTAGCTCGCCGATTGCGGGCACACAtt                           <  1:3231182/68‑1 (MQ=255)
                     cTGCCGTATACCTCCGCAGGCCAGACAGGTTGGCGTAGCTCGCCGATTGCGGGCACACATTTATAGCCa                   <  1:924579/69‑1 (MQ=255)
                     cTGCCGTATACCTCCGCAGGCCAGACAGGTTGGCGTAGCTCGCCGATTGCGGGCACACATTTATAGCCa                   <  1:1718848/69‑1 (MQ=255)
                          gTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGcc                                            >  1:907978/1‑39 (MQ=38)
                          gTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGCCg                                           >  1:2723589/1‑40 (MQ=255)
                          gTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGCCg                                           >  1:1324089/1‑40 (MQ=255)
                          gTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGCCg                                           >  1:2747654/1‑40 (MQ=255)
                          gTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGCCg                                           >  1:501845/1‑40 (MQ=255)
                          gTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGCCg                                           >  1:535365/1‑40 (MQ=255)
                          gTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGCCg                                           >  1:547107/1‑40 (MQ=255)
                          gTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGCCg                                           >  1:611268/1‑40 (MQ=255)
                          gTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGCCg                                           >  1:627802/1‑40 (MQ=255)
                          gTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGCCg                                           >  1:1307723/1‑40 (MQ=255)
                          gTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGCCg                                           >  1:2412729/1‑40 (MQ=255)
                          gTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGCCg                                           >  1:1398994/1‑40 (MQ=255)
                          gTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGCCg                                           >  1:1538609/1‑40 (MQ=255)
                          gTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGCCg                                           >  1:1611816/1‑40 (MQ=255)
                          gTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGCCg                                           >  1:171937/1‑40 (MQ=255)
                          gTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGCCg                                           >  1:2029426/1‑40 (MQ=255)
                          gTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGCCg                                           >  1:2608165/1‑40 (MQ=255)
                          gTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGCCg                                           >  1:1873060/1‑40 (MQ=255)
                          gTATACCCCCGCAGGCCACACAGGTTGGCGTAGCTCGCCg                                           >  1:1756182/1‑40 (MQ=255)
                               ccTCCGCAGGCCAGACAGGTTGGCGTAGCTCGCCGATTGCGGGCACACATTTATAGCCATATTCCCg           >  1:2796758/1‑67 (MQ=255)
                               ccTCCGCAGGCCAGACAGGTTGGCGTAGCTCGCCGATTGCGGGCACACATTTATAGCCATATTCCCg           >  1:2828972/1‑67 (MQ=255)
                                 tCCGCAGGCCAGACAGGTTGGCGTAGCTCGCCGATTGCGGGCACACATTTATAGCCATATTCCCGCaa        >  1:975566/1‑68 (MQ=255)
                                  ccGCAGGCCAGACAGGTTGGCGTAGCTCGCCGATTGCGGGCACACATTTATAGCCATATTCCCGCAAAg      <  1:1937173/69‑1 (MQ=255)
                                      aGGCCAGACAGGTTGGCGTAGCTCGCCGATTGCGGGCACACATTTATAGCCATATTCCCGCAAAGCATa  <  1:365064/69‑1 (MQ=255)
                                      aGGCCAGACAGGTTGGCGTAGCTCGCCGATTGCGGGCACACATTTATAGCCATATTCCCGCAAAGCATa  <  1:583609/69‑1 (MQ=255)
                                      aGGCCAGACAGGTTGGCGTAGCTCGCCGATTGCGGGCACACATTTATAGCCATATTCCCGCAAAGCATa  <  1:1836675/69‑1 (MQ=255)
                                      aGGCCAGACAGGTTGGCGTAGCTCGCCGATTGCGGGCACACATTTATAGCCATATTCCCGCAAAGCATa  <  1:218543/69‑1 (MQ=255)
                                      aGGCCAGACAGGTTGGCGTAGCTCGCCGATTGCGGGCACACATTTATAGCCATATTCCCGCAAAGCATa  <  1:3224092/69‑1 (MQ=255)
                                      aGGCCAGACAGGTTGGCGTAGCTCGCCGATTGCGGGCACACATTTATAGCCATATTCCCGCAAAGCATa  <  1:2324045/69‑1 (MQ=255)
                                      aGGCCAGACAGGTTGGCGTAGCTCGCCGATTGCGGGCACACATTTATAGCCATATTCCCGCAAAGCATa  <  1:2995032/69‑1 (MQ=255)
                                      aGGCCAGACAGGTTGGCGTAGCTCGCCGATTGCGGGCACACATTTATAGCCATATTCCCGCAAAGCATa  <  1:2434770/69‑1 (MQ=255)
                                      aGGCCAGACAGGTTGGCGTAGCTCGCCGATTGCGGGCACACATTTATAGCCATATTCCCGCAAAGCATa  <  1:2605144/69‑1 (MQ=255)
                                            |                                                              
GTAGTCCCACAGTGGCTAATACTGCCGTATACCTCCGCAGGCCAGACAGGTTGGCGTAGCTCGCCGATTGCGGGCACACATTTATAGCCATATTCCCGCAAAGCATA  >  minE/364696‑364802

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: