Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,068,6000GA33.3% 89.0 / 70.6 75A397V (GCC→GTC) aerfused signal transducer for aerotaxis sensory component and methyl accepting chemotaxis component
Reads supporting (aligned to +/- strand):  ref base G (50/0);  new base A (25/0);  total (75/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 = 8.87e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTC  >  minE/2068552‑2068649
                                                |                                                 
gCATCAATCAGCTTGCGTATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:3300219/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTaa                                               >  1:2774918/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTaa                                               >  1:833413/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTaa                                               >  1:77500/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTaa                                               >  1:676103/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTaa                                               >  1:447434/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTaa                                               >  1:3630915/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTaa                                               >  1:3167316/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTaa                                               >  1:3121612/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTaa                                               >  1:3117355/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTaa                                               >  1:2763090/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTaa                                               >  1:2333533/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTaa                                               >  1:1817597/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTaa                                               >  1:1756842/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTaa                                               >  1:1029324/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTaa                                               >  1:1562149/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTaa                                               >  1:1559298/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTaa                                               >  1:1242719/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTaa                                               >  1:1305392/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTa                                                >  1:308433/1‑52 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTa                                                >  1:543561/1‑52 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:2826944/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:1432529/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:2924048/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:1576901/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:136137/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:320530/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:335141/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:3431879/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:3502720/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:3571616/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:380489/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:411391/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:418489/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:573530/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:589008/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:1072236/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:2331010/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:2148547/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:1740849/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:154912/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:2403825/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:1659867/1‑53 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGACTa                                                >  1:1981319/1‑52 (MQ=255)
gCATCAATCAGCTTGCGAATATCGTAGGCAGCATTAGCGCTGCGGCTGACTaa                                               >  1:482471/1‑53 (MQ=39)
                                    gcgcTGCGGCTGGTTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:1341899/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGccacca                        >  1:2241204/1‑40 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCGAAACCTTTGCCCTGTTc  >  1:3143369/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGtt   >  1:1591642/1‑61 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:816099/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:1942023/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:1934711/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:1791014/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:548777/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:3449628/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:1668811/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:1181347/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:776130/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:703089/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:2121200/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:2164658/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:1012771/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:2171812/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:3383288/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:2180816/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:3322382/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:3290433/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:3133882/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:2379769/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:2887887/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:2842046/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:1518741/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCCTTGCCCTGTTc  >  1:1364651/1‑62 (MQ=255)
                                    gcgcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAAACTTTGCCCTGTTc  >  1:188794/1‑62 (MQ=255)
                                    gagcTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTc  >  1:2839588/3‑62 (MQ=255)
                                                |                                                 
GCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTC  >  minE/2068552‑2068649

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: