Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE717,2300GA71.1% 66.3 / 33.2 45E82E (GAG→GAAyceHconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (10/3);  new base A (16/16);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCAC  >  minE/717163‑717295
                                                                   |                                                                 
aTAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCaa                                                                >  1:731974/1‑71 (MQ=255)
aTAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCaa                                                                >  1:151261/1‑71 (MQ=255)
aTAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGCAATCGGGTCACCAAATATGAGCaa                                                                >  1:1290321/1‑71 (MQ=255)
     cTGGTCAAACGTCATTATCTACTCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:278361/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCTCCAAATATGAg                                                                   >  1:1718413/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:3265765/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:602774/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:550782/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:429493/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:34546/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:1978571/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:327636/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:2116106/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:1024244/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:2187370/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:2211861/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:2329714/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:236302/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:2406945/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:2961619/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:2715863/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:1863756/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:1173035/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:1008027/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:925950/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:915087/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:1375042/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:806403/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:802643/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:1619399/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:780007/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:1667021/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:1783298/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAg                                                                   >  1:635874/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTCGTAATCGGGTCACCAAATATGAg                                                                   >  1:1921635/1‑63 (MQ=255)
     cTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTATGGTAATCGGGTCACCAAATATGAg                                                                   >  1:3142892/1‑63 (MQ=255)
              cgccATTATCTGCGCACGGTGAGTGGCTTTGGTAATCGGGTCACCAAATATGAACAACGTTTTTGTaa                                                     >  1:661247/4‑68 (MQ=255)
              cgccATTATCTGCGCACGGTGAGTGGCTTTGGTAATCGGGTCACCAAATATGAACAACGTTTTTGTAAtt                                                   >  1:327632/4‑70 (MQ=255)
              cgccATTATCTGCGCACGGTGAGTGGCTTTGGTAATCGGGTCACCAAATATGAACAACGTTTTTGTAAtt                                                   >  1:3247634/4‑70 (MQ=255)
              cgccATTATCTGCGCACGGTGAGTGGCTTTGGTAATCGGGTCACCAAATATGAACAACGTTTTTGTAAtt                                                   >  1:929201/4‑70 (MQ=255)
              cgccATTATCTGCGCACGGTGAGTGGCTTTGGTAATCGGGTCACCAAATATGAACAACGTTTTTGTAAtt                                                   >  1:3111631/4‑70 (MQ=255)
              cgccATTATCTGCGCACGGTGAGTGGCTTTGGTAATCGGGTCACCAAATATGAACAACGTTTTTGTAAtt                                                   >  1:2831169/4‑70 (MQ=255)
              cgccATTATCTGCGCACGGTGAGTGGCTTTGGTAATCGGGTCACCAAATATGAACAACGTTTTTGTAAtt                                                   >  1:96016/4‑70 (MQ=255)
              cgccATTATCTGCGCACGGTGAGTGGCTTTGGTAATCGGGTCACCAAATATGAACAACGTTTTTGTAAtt                                                   >  1:1816543/4‑70 (MQ=255)
              cgccATTATCTGCGCACGGTGAGTGGCTTTGGTAATCGGGTCACCAAATATGAACAACGTTTTTGTAAtt                                                   >  1:2657592/4‑70 (MQ=255)
              cgccATTATCTGCGCACGGTGAGTGGCTTTGGTAATCGGGTCACCAAATATGAACAACGTTTTTGTAAtt                                                   >  1:1635595/4‑70 (MQ=255)
              cgccATTATCTGCGCACGGTGAGTGGCTTTGGTAATCGGGTCACCAAATATGAACAACGTTTTTGTAAtt                                                   >  1:1956637/4‑70 (MQ=255)
              cgccATTATCTGCGCACGGTGAGTGGCTTTGGTAATCGGGTCACCAAATATGAACAACGTTTTTGTAAtt                                                   >  1:2562564/4‑70 (MQ=255)
              cgccATTATCTGCGCACGGTGAGTGGCTTTGGTAATCGGGTCACCAAATATGAACAACGTTTTTGTAAtt                                                   >  1:2385938/4‑70 (MQ=255)
              cgccATTATCTGCGCACGGTGAGTGGCTTTGGTAATCGGGTCACCAAATATGAACAACGTTTTTGTAAtt                                                   >  1:2070471/4‑70 (MQ=255)
              cgccATTATCTGCGCACGGTGAGTGGCTTTGGTAATCGGGTCACCAAATATGAACAACGTTTTTGTAAtt                                                   >  1:1379576/4‑70 (MQ=255)
              cgccATTATCTGCGCACGGTGAGTGGCTTTGGTAATCGGGTCACCAAATATGAACAACGTATTTGTAAtt                                                   >  1:1316568/4‑70 (MQ=255)
                                                aTCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGa                 <  1:1331204/70‑1 (MQ=255)
                                                aTCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGa                 <  1:2269687/70‑1 (MQ=255)
                                                aTCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGa                 <  1:2154795/70‑1 (MQ=255)
                                                 tCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGaa                >  1:1359975/1‑70 (MQ=255)
                                                 tCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAg               >  1:2680930/1‑71 (MQ=255)
                                                 tCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAg               >  1:1819956/1‑71 (MQ=255)
                                                 tCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAg               >  1:1479605/1‑71 (MQ=255)
                                                 tCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAg               >  1:1970730/1‑71 (MQ=255)
                                                 tCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAg               >  1:2931232/1‑71 (MQ=255)
                                                 tCGGGTCACCAAATATGAGCAACGTTTTTATAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAg               >  1:649624/1‑71 (MQ=255)
                                                              tatGAACAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCCGCAGAAGTGGCGTTAATcac  <  1:879847/71‑1 (MQ=255)
                                                              tatGAACAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCCGCAGAAGTGGCGTTAATcac  <  1:931625/71‑1 (MQ=255)
                                                              tatGAACAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCCGCAGAAGTGGCGTTAATcac  <  1:1637731/71‑1 (MQ=255)
                                                              tatGAACAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCCGCAGAAGTGGCGTTAATcac  <  1:1240668/71‑1 (MQ=255)
                                                              tatGAACAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCCGCAGAAGTGGCGTTAATcac  <  1:2678255/71‑1 (MQ=255)
                                                              tatGAACAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCCGCAGAAGTGGCGTTAATcac  <  1:722962/71‑1 (MQ=255)
                                                              tatGAACAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCCGCAGAAGTGGCGTTAATcac  <  1:1954470/71‑1 (MQ=255)
                                                              tatGAACAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCCGCAGAAGTGGCGTTAATcac  <  1:195567/71‑1 (MQ=255)
                                                              tatGAACAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCCGCAGAAGTGGCGTTAATcac  <  1:3292399/71‑1 (MQ=255)
                                                              tatGAACAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCCGCAGAAGTGGCGTTAATcac  <  1:3268871/71‑1 (MQ=255)
                                                              tatGAACAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCCGCAGAAGTGGCGTTAATcac  <  1:2185020/71‑1 (MQ=255)
                                                              tatGAACAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCCGCAGAAGTGGCGTTAATcac  <  1:3079802/71‑1 (MQ=255)
                                                              tatGAACAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCCGCAGAAGTGGCGTTAATcac  <  1:3012152/71‑1 (MQ=255)
                                                              tatGAACAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCCGCAGAAGTGGCGTTAATcac  <  1:2566699/71‑1 (MQ=255)
                                                              tatGAACAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCCGCAGAAGTGGCGTTAATcac  <  1:2916895/71‑1 (MQ=255)
                                                               atGAACAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCCGCAGAAGTGGCGTTAATcac  <  1:360193/70‑1 (MQ=255)
                                                                   |                                                                 
ATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCAC  >  minE/717163‑717295

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: