Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE958,1740AG60.6% 44.6 / 77.9 66V711A (GTC→GCC) pqqLpredicted peptidase
Reads supporting (aligned to +/- strand):  ref base A (26/0);  new base G (40/0);  total (66/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATC  >  minE/958111‑958229
                                                               |                                                       
gAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGa                                                         >  1:2038009/1‑64 (MQ=255)
gAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACTTTAc                                                   >  1:1524557/1‑70 (MQ=255)
gAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACAtt                                                     >  1:1443605/1‑68 (MQ=255)
gAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAcc                                                  >  1:2366488/1‑71 (MQ=255)
gAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAcc                                                  >  1:2501605/1‑71 (MQ=255)
gAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAcc                                                  >  1:2504428/1‑71 (MQ=255)
gAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAcc                                                  >  1:2245976/1‑71 (MQ=255)
gAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAcc                                                  >  1:3140202/1‑71 (MQ=255)
gAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAcc                                                  >  1:1619031/1‑71 (MQ=255)
gAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAcc                                                  >  1:2132047/1‑71 (MQ=255)
gAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAcc                                                  >  1:1556754/1‑71 (MQ=255)
gAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAcc                                                  >  1:1380135/1‑71 (MQ=255)
gAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAcc                                                  >  1:344986/1‑71 (MQ=255)
gAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAcc                                                  >  1:404621/1‑71 (MQ=255)
gAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAcc                                                  >  1:1232618/1‑71 (MQ=255)
gAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAcc                                                  >  1:416646/1‑71 (MQ=255)
gAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAcc                                                  >  1:509353/1‑71 (MQ=255)
gAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAcc                                                  >  1:596668/1‑71 (MQ=255)
     agagTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACaa                                             >  1:2105656/1‑71 (MQ=255)
     agagTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACaa                                             >  1:2097328/1‑71 (MQ=255)
     agagTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACaa                                             >  1:2044376/1‑71 (MQ=255)
     agagTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACaa                                             >  1:2512489/1‑71 (MQ=255)
     agagTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACaa                                             >  1:1332860/1‑71 (MQ=255)
     agagTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACaa                                             >  1:3187791/1‑71 (MQ=255)
     agagTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACaa                                             >  1:704484/1‑71 (MQ=255)
     agagTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCATCTGCGACATTACCGACaa                                             >  1:3179248/1‑71 (MQ=255)
     agagTGTTAGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACaa                                             >  1:971329/1‑71 (MQ=255)
                                                aGTTTGTCTTGTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:525945/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACaaa                                      >  1:1331645/1‑35 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGatat                              >  1:2893914/1‑43 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGaat               >  1:2019286/1‑57 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAAc             >  1:580544/1‑60 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTgg       >  1:892804/1‑66 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGAt   >  1:2829595/1‑70 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:1034644/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:2024170/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:3244890/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:3245356/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:3271124/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:1261898/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:1245356/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:1215761/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:1211726/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:1173775/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:697386/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:1124614/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:745628/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:747207/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:874188/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:1851436/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:1974647/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:2208805/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:1971371/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:2362633/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:1940591/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:2457588/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:2477190/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:2486619/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:3229179/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:1795662/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:2700070/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:2700790/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:2746185/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:2827976/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:1752838/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATc  >  1:1572825/1‑71 (MQ=255)
                                                aGTTTGTCTTCTGCGGCATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTCGCGATc  >  1:1790573/1‑71 (MQ=255)
                                                               |                                                       
GAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATC  >  minE/958111‑958229

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: