Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb716,4311.T61.3% 58.9 / 43.0 31coding (349/411 nt)ybfGhypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (12/0);  new base T (0/19);  total (12/19)
Fisher's exact test for biased strand distribution p-value = 7.09e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGG  >  W3110S.gb/716367‑716473
                                                                 |                                          
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAAc                                             >  1:2253790/1‑65 (MQ=255)
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAAc                                             >  1:909877/1‑65 (MQ=255)
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAAc                                             >  1:89059/1‑65 (MQ=255)
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAAc                                             >  1:705621/1‑65 (MQ=255)
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAAc                                             >  1:672981/1‑65 (MQ=255)
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAAc                                             >  1:495765/1‑65 (MQ=255)
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAAc                                             >  1:459028/1‑65 (MQ=255)
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAAc                                             >  1:2270631/1‑65 (MQ=255)
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAAc                                             >  1:1008311/1‑65 (MQ=255)
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAAc                                             >  1:2067878/1‑65 (MQ=255)
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAAc                                             >  1:2029416/1‑65 (MQ=255)
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAAc                                             >  1:188675/1‑65 (MQ=255)
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAAc                                             >  1:1689292/1‑65 (MQ=255)
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAAc                                             >  1:1650122/1‑65 (MQ=255)
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAAc                                             >  1:1513151/1‑65 (MQ=255)
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAAc                                             >  1:1069638/1‑65 (MQ=255)
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAAc                                             >  1:1099780/1‑65 (MQ=255)
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAAc                                             >  1:1256009/1‑65 (MQ=255)
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAAc                                             >  1:1349117/1‑65 (MQ=255)
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAAc                                             >  1:1427477/1‑65 (MQ=255)
tCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGATCAGTAAGCGAAATAAc                                             >  1:1786707/1‑65 (MQ=255)
        cGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTttcttg                                    >  1:569763/1‑65 (MQ=255)
        cGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTttcttg                                    >  1:1414510/1‑65 (MQ=255)
        cGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTttcttg                                    >  1:736915/1‑65 (MQ=255)
        cGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTttcttg                                    >  1:1830248/1‑65 (MQ=255)
        cGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTttctt                                     >  1:172350/1‑64 (MQ=255)
        cGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTttctt                                     >  1:1298087/1‑64 (MQ=255)
        cGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTttctt                                     >  1:1323593/1‑64 (MQ=255)
        cGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTttctt                                     >  1:264322/1‑64 (MQ=255)
        cGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTttctt                                     >  1:72640/1‑64 (MQ=255)
                  tctttcATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGcttc                          >  1:1379810/1‑65 (MQ=255)
                  tctttcATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGcttc                          >  1:1306599/1‑65 (MQ=255)
                  tctttcATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGcttc                          >  1:1136320/1‑65 (MQ=255)
                                           cgagTTCAGTAAGCGAAATAACTTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAgg  <  1:2325212/64‑1 (MQ=255)
                                           agagTTCAGTAAGCGAAATAACTTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAgg  <  1:70393/65‑1 (MQ=255)
                                           agagTTCAGTAAGCGAAATAACTTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAgg  <  1:731827/65‑1 (MQ=255)
                                           agagTTCAGTAAGCGAAATAACTTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAgg  <  1:286936/65‑1 (MQ=255)
                                           agagTTCAGTAAGCGAAATAACTTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAgg  <  1:861668/65‑1 (MQ=255)
                                           agagTTCAGTAAGCGAAATAACTTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAgg  <  1:863317/65‑1 (MQ=255)
                                           agagTTCAGTAAGCGAAATAACTTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAgg  <  1:907839/65‑1 (MQ=255)
                                           agagTTCAGTAAGCGAAATAACTTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAgg  <  1:1569333/65‑1 (MQ=255)
                                           agagTTCAGTAAGCGAAATAACTTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAgg  <  1:630078/65‑1 (MQ=255)
                                           agagTTCAGTAAGCGAAATAACTTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAgg  <  1:2512469/65‑1 (MQ=255)
                                           agagTTCAGTAAGCGAAATAACTTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAgg  <  1:2415895/65‑1 (MQ=255)
                                           agagTTCAGTAAGCGAAATAACTTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAgg  <  1:2415133/65‑1 (MQ=255)
                                           agagTTCAGTAAGCGAAATAACTTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAgg  <  1:2413671/65‑1 (MQ=255)
                                           agagTTCAGTAAGCGAAATAACTTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAgg  <  1:2315765/65‑1 (MQ=255)
                                           agagTTCAGTAAGCGAAATAACTTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAgg  <  1:2300982/65‑1 (MQ=255)
                                           agagTTCAGTAAGCGAAATAACTTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAgg  <  1:1813524/65‑1 (MQ=255)
                                           agagTTCAGTAAGCGAAATAACTTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAgg  <  1:2137632/65‑1 (MQ=255)
                                           agagTTCAGTAAGCGAAATAACTTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAgg  <  1:2019496/65‑1 (MQ=255)
                                                   gTAAGGGAAATAACTTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAgg  <  1:1931550/57‑1 (MQ=255)
                                                                 |                                          
TCTTATTGCGTGAATGTTTCTTTCATTTGATACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGG  >  W3110S.gb/716367‑716473

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: