Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb63,9480GA73.9% 32.8 / 10.7 23T611T (ACC→ACTpolBDNA polymerase II
Reads supporting (aligned to +/- strand):  ref base G (0/6);  new base A (17/0);  total (17/6)
Fisher's exact test for biased strand distribution p-value = 9.91e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AACACCATCCGCTGCTTGTCGCCCTCCTGAATCAGTCCGGCATAACGCTTTTTACTGCCGGTATCGGCTCCGCGAATGGTTGGCATCAGAAAACGGCAGAAATGGGTTTCATA  >  W3110S.gb/63889‑64001
                                                           |                                                     
aaCACCATCCGCTGCTTGTCGCCCTCCTGAATCAGTCCGGCATAACGCTTTTTACTGCCgg                                                      <  1:672084/61‑1 (MQ=255)
aaCACCATCCGCTGCTTGTCGCCCTCCTGAATCAGTCCGGCATAACGCTTTTTACTGCCgg                                                      <  1:190691/61‑1 (MQ=255)
aaCACCATCCGCTGCTTGTCGCCCTCCTGAATCAGTCCGGCATAACGCTTTTTACTGCCgg                                                      <  1:2605965/61‑1 (MQ=255)
                                 aGTGCGGCATAACGCTTTTTACTGCCAGTATCGGCTCCGCGAATGGTTGGCATCAGAAAAc                     >  1:363943/1‑61 (MQ=255)
                                 aGTCCGGCATAACGCTTTTTACTGGCAGTATCGGCTCCGCGAATGGTTGGCATCAGAAAAc                     >  1:1346438/1‑61 (MQ=255)
                                 aGTCCGGCATAACGCTTTTTACTGCCAGTATCGGCTCCGCGAATGGTTGTCATCAGAAAAc                     >  1:1160028/1‑61 (MQ=255)
                                 aGTCCGGCATAACGCTTTTTACTGCCAGTATCGGCTCCGCGAATGGTTGGCATCAGAAAAc                     >  1:891212/1‑61 (MQ=255)
                                 aGTCCGGCATAACGCTTTTTACTGCCAGTATCGGCTCCGCGAATGGTTGGCATCAGAAAAc                     >  1:1069032/1‑61 (MQ=255)
                                 aGTCCGGCATAACGCTTTTTACTGCCAGTATCGGCTCCGCGAATGGTTGGCATCAGAAAAc                     >  1:784398/1‑61 (MQ=255)
                                 aGTCCGGCATAACGCTTTTTACTGCCAGTATCGGCTCCGCGAATGGTTGGCATCAGAAAAc                     >  1:722324/1‑61 (MQ=255)
                                 aGTCCGGCATAACGCTTTTTACTGCCAGTATCGGCTCCGCGAATGGTTGGCATCAGAAAAc                     >  1:374776/1‑61 (MQ=255)
                                 aGTCCGGCATAACGCTTTTTACTGCCAGTATCGGCTCCGCGAATGGTTGGCATCAGAAAAc                     >  1:368154/1‑61 (MQ=255)
                                 aGTCCGGCATAACGCTTTTTACTGCCAGTATCGGCTCCGCGAATGGTTGGCATCAGAAAAc                     >  1:2933827/1‑61 (MQ=255)
                                 aGTCCGGCATAACGCTTTTTACTGCCAGTATCGGCTCCGCGAATGGTTGGCATCAGAAAAc                     >  1:2678174/1‑61 (MQ=255)
                                 aGTCCGGCATAACGCTTTTTACTGCCAGTATCGGCTCCGCGAATGGTTGGCATCAGAAAAc                     >  1:2625848/1‑61 (MQ=255)
                                 aGTCCGGCATAACGCTTTTTACTGCCAGTATCGGCTCCGCGAATGGTTGGCATCAGAAAAc                     >  1:2159622/1‑61 (MQ=255)
                                 aGTCCGGCATAACGCTTTTTACTGCCAGTATCGGCTCCGCGAATGGTTGGCATCAGAAAAc                     >  1:1655174/1‑61 (MQ=255)
                                 aGTCCGGCATAACGCTTTTTACTGCCAGTATCGGCTCCGCGAATGGTTGGCATCAGAAAAc                     >  1:1448618/1‑61 (MQ=255)
                                 aGTCCGGCATAACGCTTTTTACTGCCAGTATCGGCTCCGCGAATGGTTGGCATCAGAAAAc                     >  1:1440962/1‑61 (MQ=255)
                                 aGTCCGGCATAACGCTTTTTACTGCCAGTATCGGCTCCGCGAATGGTTGGCATCAGAAAAc                     >  1:1180095/1‑61 (MQ=255)
                                                    tACTGCCGGTATCGGCTCCGCGAATGGTTGGCATCAGAAAACGGCAGAAATGGGTTTCATa  <  1:2884113/61‑1 (MQ=255)
                                                    tACTGCCGGTATCGGCTCCGCGAATGGTTGGCATCAGAAAACGGCAGAAATGGGTTTCATa  <  1:1505298/61‑1 (MQ=255)
                                                    tACTGCCGGTATCGGCTCCGCGAATGGTTGGCATCAGAAAACGGCAGAAATGGGTTTCATa  <  1:428311/61‑1 (MQ=255)
                                                    tACTGCCGGTATCGGCTCCGCGAATGGTTGGCATCAGAAAACGGCAGAAATGGGTTTCATa  <  1:660103/61‑1 (MQ=255)
                                                    tACTGCCGGTATCGGCTCCGCGAATGGTTGGCATCAGAAAACGGCAGAAATGGGTTTCATa  <  1:1321207/61‑1 (MQ=255)
                                                    tACTGCCGGTATCGGCTCCGCGAATGGTTGGCATCAGAAAACGGCAGAAATGGGTTTCATa  <  1:840314/61‑1 (MQ=255)
                                                           |                                                     
AACACCATCCGCTGCTTGTCGCCCTCCTGAATCAGTCCGGCATAACGCTTTTTACTGCCGGTATCGGCTCCGCGAATGGTTGGCATCAGAAAACGGCAGAAATGGGTTTCATA  >  W3110S.gb/63889‑64001

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: