Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,974,3760CT66.7% 29.9 / 25.1 30L10L (CTG→CTAtarmethyl‑accepting chemotaxis protein II
Reads supporting (aligned to +/- strand):  ref base C (2/8);  new base T (0/20);  total (2/28)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-02
Rejected as consensus: Frequency below/above cutoff threshold.

TGCGAATACCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCT  >  W3110S.gb/1974349‑1974434
                           |                                                          
tGCGAATACCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCAcc                       >  1:473879/1‑65 (MQ=255)
tGCGAATACCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCAcc                       >  1:1665770/1‑65 (MQ=255)
             aGCTCCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACaa          <  1:691633/65‑1 (MQ=255)
             aGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACaa          <  1:1097280/65‑1 (MQ=255)
             aGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACaa          <  1:1882024/65‑1 (MQ=255)
             aGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACaa          <  1:677992/65‑1 (MQ=255)
             aGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACaa          <  1:703213/65‑1 (MQ=255)
             aGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACaa          <  1:1476131/65‑1 (MQ=255)
             aGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACaa          <  1:715795/65‑1 (MQ=255)
             aGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACaa          <  1:1187017/65‑1 (MQ=255)
                     tACCAATAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGAtct  <  1:789959/65‑1 (MQ=255)
                     tACCAATAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGAtct  <  1:100158/65‑1 (MQ=255)
                     tACCAATAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGAtct  <  1:590022/65‑1 (MQ=255)
                     tACCAATAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGAtct  <  1:922826/65‑1 (MQ=255)
                     tACCAATAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGAtct  <  1:46246/65‑1 (MQ=255)
                     tACCAATAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGAtct  <  1:388730/65‑1 (MQ=255)
                     tACCAATAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGAtct  <  1:2864919/65‑1 (MQ=255)
                     tACCAATAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGAtct  <  1:285744/65‑1 (MQ=255)
                     tACCAATAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGAtct  <  1:2765305/65‑1 (MQ=255)
                     tACCAATAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGAtct  <  1:2647858/65‑1 (MQ=255)
                     tACCAATAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGAtct  <  1:2468738/65‑1 (MQ=255)
                     tACCAATAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGAtct  <  1:221959/65‑1 (MQ=255)
                     tACCAATAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGAtct  <  1:1938907/65‑1 (MQ=255)
                     tACCAATAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGAtct  <  1:1789820/65‑1 (MQ=255)
                     tACCAATAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGAtct  <  1:1605696/65‑1 (MQ=255)
                     tACCAATAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGAtct  <  1:138582/65‑1 (MQ=255)
                     tACCAATAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGAtct  <  1:1376435/65‑1 (MQ=255)
                     tACCAATAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGAtct  <  1:1368903/65‑1 (MQ=255)
                      aCCAATAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGAtct  <  1:2011414/64‑1 (MQ=255)
                      aCCAATAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGAtct  <  1:1813047/64‑1 (MQ=255)
                           |                                                          
TGCGAATACCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCT  >  W3110S.gb/1974349‑1974434

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: