Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,312,3310CA57.9% ‑1.3 / 92.1 57A259E (GCG→GAG) menA1,4‑dihydroxy‑2‑naphthoate octaprenyltransferase
Reads supporting (aligned to +/- strand):  ref base C (22/2);  new base A (0/33);  total (22/35)
Fisher's exact test for biased strand distribution p-value = 7.91e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCTGGTGTGTCTGGCGCTGTTTAATCTCTTTTCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAATGGACCCGGTGGCGATGCGACCAA  >  minE/2312263‑2312399
                                                                    |                                                                    
gCTGGTGTGTCTGGCGCTGTTTAATCTCTTTTCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTggc                                                                      <  1:3268976/69‑1 (MQ=255)
gCTGGTGTGTCTGGCGCTGTTTAATCTCTTTTCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTggc                                                                      <  1:1828871/69‑1 (MQ=255)
gCTGGTGTGTCTGGCGCTGTTTAATCTCTTTTCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTggc                                                                      <  1:2616172/69‑1 (MQ=255)
gCTGGTGTGTCTGGCGCTGTTTAATCTCTTTTCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTggc                                                                      <  1:813255/69‑1 (MQ=255)
gCTGGTGTGTCTGGCGCTGTTTAATCTCTTTTCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTggc                                                                      <  1:1917227/69‑1 (MQ=255)
gCTGGTGTGTCTGGCGCTGTTTAATCTCTTTTCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTggc                                                                      <  1:1957313/69‑1 (MQ=255)
gCTGGTGTGTCTGGCGCTGTTTAATCTCTTTTCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTggc                                                                      <  1:2448802/69‑1 (MQ=255)
gCTGGTGTGTCTGGCGCTGTTTAATCTCTTTTCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTggc                                                                      <  1:381661/69‑1 (MQ=255)
gCTGGTGTGTCTGGCGCTGTTTAATCTCTTTTCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTggc                                                                      <  1:1258950/69‑1 (MQ=255)
gCTGGTGTGTCTGGCGCTGTTTAATCTCTTTTCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTggc                                                                      <  1:61033/69‑1 (MQ=255)
gCTGGTGTGTCTGGCGCTGTTTAATCTCTTTTCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTggc                                                                      <  1:3107690/69‑1 (MQ=255)
gCTGGTGTGTCTGGCGCTGTTTAATCTCTTTTCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTggc                                                                      <  1:1177951/69‑1 (MQ=255)
gCTGGTGTGTCTGGCGCTGTTTAATCTCTTTTCGCTGCATAGCCTGTGGGGCTGGCTGTTCCGGCTggc                                                                      <  1:420771/69‑1 (MQ=255)
gCTGGTGTCTCTGGCGCTGTTTAATCTCTTTTCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTggc                                                                      <  1:336033/69‑1 (MQ=255)
 cTGGTGTGTCTGGCGCTGTTTAATCTCTTTTCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTTCTggc                                                                      <  1:3053663/68‑1 (MQ=255)
              cgcTGTTTAATCTCTTTTCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGCGGCACCATTACTGGTg                                                      <  1:316260/71‑1 (MQ=255)
              cgcTGTTTAATCTCTTTTCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGCGGCACCATTACTGGTg                                                      <  1:1655574/71‑1 (MQ=255)
                     tAATCTCTTTTCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTggc                                                                      <  1:1268825/48‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:558146/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:2472116/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:2221920/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:2383760/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:2476429/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:2537116/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:2638125/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:2730838/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:2802827/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:2956055/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:2959780/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:3292288/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:3273430/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:839545/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:1252910/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:1255475/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:1330045/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:1386700/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:1418694/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:849580/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:1480940/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:1525636/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:795055/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:1209581/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:2039064/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:63482/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:2370880/69‑1 (MQ=255)
                              ttCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCACGTTa                                        <  1:1493729/69‑1 (MQ=255)
                              ttCGCTGCATAGACTGTGGGGCTGGCTGTTCATGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:1579206/69‑1 (MQ=255)
                              ttCGCTGCAGAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:1460122/69‑1 (MQ=255)
                               tCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:3155333/68‑1 (MQ=255)
                               tCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:2508015/68‑1 (MQ=255)
                                            tgtgGGGCTGGCTGTTCCTGCTGGAGGCACCATTACTGGTGAAGCAAGCCCGTTa                                        <  1:1008734/55‑1 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:894773/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:887729/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:534655/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:667629/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:2491592/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:605024/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:2367424/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:1131254/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:1154888/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:1541300/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:1561683/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:1648066/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:1804144/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:2187554/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:3227814/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:240341/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:2449648/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:257059/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:2955710/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:2986208/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:302809/1‑57 (MQ=255)
                                                         gTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAAt                         >  1:3076835/1‑57 (MQ=255)
                                                                  ggcggcACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAATGGACCCGGTGGCGATGCGACCaa  >  1:14504/1‑71 (MQ=255)
                                                                    |                                                                    
GCTGGTGTGTCTGGCGCTGTTTAATCTCTTTTCGCTGCATAGCCTGTGGGGCTGGCTGTTCCTGCTGGCGGCACCATTACTGGTGAAGCAAGCCCGTTATGTGATGCGGGAAATGGACCCGGTGGCGATGCGACCAA  >  minE/2312263‑2312399

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: