Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,821,1530TA71.4% 65.1 / 32.6 42E48V (GAG→GTG) recXregulatory protein for RecA
Reads supporting (aligned to +/- strand):  ref base T (9/3);  new base A (30/0);  total (39/3)
Fisher's exact test for biased strand distribution p-value = 1.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCAC  >  W3110S.gb/2821123‑2821213
                              |                                                            
cgcTCGTGATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:965155/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCtcttc                                                          >  1:1124030/1‑35 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:1218752/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:764522/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:409254/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:1656871/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGc                               >  1:2338598/1‑62 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCGATAATCGGTGcc                              >  1:2236382/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:1033102/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:2209831/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:2230643/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:226395/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:2357382/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:2470863/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:2514639/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:339600/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:457905/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:7520/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:828137/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:982048/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:986739/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:2101939/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:1051106/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:105328/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:1065499/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:1194854/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:1319118/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:1569984/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:1582567/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:193546/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:1943825/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGcc                              >  1:1951277/1‑63 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGc                               >  1:1696552/1‑62 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGc                               >  1:2396128/1‑62 (MQ=255)
cgcTCGTAATCTTCTGCCGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATAATCGGTGc                               >  1:2357898/1‑62 (MQ=255)
cgcTCGTAATCTTCTGCAGTAGCATCAATCACTTCTGGGCCATTTTTGCCCATaa                                      >  1:375359/1‑55 (MQ=255)
              tGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCa              >  1:298605/1‑65 (MQ=255)
              tGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCa              >  1:1355564/1‑65 (MQ=255)
                   tAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTc          >  1:1746685/1‑64 (MQ=255)
                   tAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCt         >  1:2125627/1‑65 (MQ=255)
                          aaTCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCAc  <  1:1643721/65‑1 (MQ=255)
                          aaTCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCAc  <  1:1590410/65‑1 (MQ=255)
                          aaTCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCAc  <  1:468919/65‑1 (MQ=255)
                              |                                                            
CGCTCGTAATCTTCTGCCGTAGCATCAATCTCTTCTGGGCCATTTTTGCCCATAATCGGTGCCGCGAGTTTACGTCGCAGTTCTTGCTCAC  >  W3110S.gb/2821123‑2821213

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: ‑