Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,069,5391.G59.1% 42.1 / 29.0 22coding (442/1974 nt)glgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base . (4/5);  new base G (13/0);  total (17/5)
Fisher's exact test for biased strand distribution p-value = 4.78e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGCAGCACCATCATTTATGA  >  W3110S.gb/4069477‑4069561
                                                               |                      
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:2984321/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:2247470/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:2331182/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:2453009/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:2498991/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:2518429/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:2555231/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:2589352/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:2814807/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:2933434/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:1076895/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:2993217/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:3072713/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:3097360/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:533223/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:584861/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:678245/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:72549/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:772649/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:937597/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:1632356/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:1047348/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:120622/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:1227259/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:1620769/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:2222459/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:1753837/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:181692/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:1850700/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:1995141/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTgg                      >  1:2104752/1‑65 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTg                       >  1:573534/1‑64 (MQ=255)
cGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTg                       >  1:1239606/1‑64 (MQ=255)
                  ttGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGCAGCACCATCATTTa     >  1:2374385/1‑64 (MQ=255)
                  ttGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGCAGCACCATCATTTa     >  1:1480963/1‑64 (MQ=255)
                  ttGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGCAGCACCATCATTTa     >  1:167819/1‑64 (MQ=255)
                  ttGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGCAGCACCATCATTTa     >  1:262864/1‑64 (MQ=255)
                    gATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGCAGCACCATCATTTATGa  <  1:829918/65‑1 (MQ=255)
                    gATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGCAGCACCATCATTTATGa  <  1:1281134/65‑1 (MQ=255)
                    gATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGCAGCACCATCATTTATGa  <  1:2297401/65‑1 (MQ=255)
                    gATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGGGGGGCAGCACCATCATTTATGa  <  1:3089678/65‑1 (MQ=255)
                     aTCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGGCAGCACCATCATTTATGa  >  1:2329147/1‑65 (MQ=255)
                     aTCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGGCAGCACCATCATTTATGa  >  1:406513/1‑65 (MQ=255)
                     aTCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGGCAGCACCATCATTTATGa  >  1:667166/1‑65 (MQ=255)
                     aTCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGGCAGCACCATCATTTATGa  >  1:862816/1‑65 (MQ=255)
                     aTCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGGCAGCACCATCATTTATGa  >  1:577733/1‑65 (MQ=255)
                     aTCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGGCAGCACCATCATTTATGa  >  1:568006/1‑65 (MQ=255)
                     aTCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGGCAGCACCATCATTTATGa  >  1:1439980/1‑65 (MQ=255)
                     aTCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGGCAGCACCATCATTTATGa  >  1:1585033/1‑65 (MQ=255)
                     aTCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGGCAGCACCATCATTTATGa  >  1:2871039/1‑65 (MQ=255)
                     aTCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGGCAGCACCATCATTTATGa  >  1:2000866/1‑65 (MQ=255)
                     aTCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGGCAGCACCATCATTTATGa  >  1:2711922/1‑65 (MQ=255)
                     aTCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGGCAGCACCATCATTTATGa  >  1:1797708/1‑65 (MQ=255)
                     aTCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGGCAGCACCATCATTTATGa  >  1:2477337/1‑65 (MQ=255)
                            tGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGCAGCACCATCATTTATGa  <  1:2352681/57‑1 (MQ=255)
                                                               |                      
CGAAATGCGTAGTGGTGGTTGATCACTATGACTGGGAAGATGATGCCCCGCCGCGCACGCCGTGGGGCAGCACCATCATTTATGA  >  W3110S.gb/4069477‑4069561

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: