Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,680,1000TC73.7% 58.6 / 21.4 38L417L (CTA→CTGilvAthreonine deaminase
Reads supporting (aligned to +/- strand):  ref base T (6/4);  new base C (0/28);  total (6/32)
Fisher's exact test for biased strand distribution p-value = 7.61e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATGCGATGGACGTCCGCCGACCATATAGCGCACGTGTAGCTTCGCCATTTCGTCGTCGGAGAGATCAACCACGCTGTAGCCGCCGTCGTTGAGC  >  W3110S.gb/3680064‑3680157
                                    |                                                         
aTGCGATGGACGTCCGCCGACCATATAGCGCACGTGTAGCTTCGCCATTTCGTCGTCGgaga                                  <  1:270738/62‑1 (MQ=255)
  gCGATGGACGTCCGCCGACCATATAGCGCACGTGTAGCTTCGCCATTTCGTCGTCGgagaga                                <  1:2564957/62‑1 (MQ=255)
  gCGATGGACGTCCGCCGACCATATAGCGCACGTGTAGCTTCGCCATTTCGTCGTCGgagaga                                <  1:2588050/62‑1 (MQ=255)
  gCGATGGACGTCCGCCGACCATATAGCGCACGTGTAGCTTCGCCATTTCGTCGTCGgagaga                                <  1:2234844/62‑1 (MQ=255)
               gccgACCATATAGCGCACGTGTAGCTTCGCCATTTCGTCGTCGGAGAGATCAACCACGCTGt                   >  1:1161914/1‑62 (MQ=255)
               gccgACCATATAGCGCACGTGTAGCTTCGCCATTTCGTCGTCGGAGAGATCAACCACGCTGt                   >  1:351397/1‑62 (MQ=255)
               gccgACCATATAGCGCACGTGTAGCTTCGCCATTTCGTCGTCGGAGAGATCAACCACGCTGt                   >  1:298267/1‑62 (MQ=255)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:252932/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:923263/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:2620943/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:2629773/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:2757118/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:2839618/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:2884248/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:331570/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:469961/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:713338/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:859617/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:1086693/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:251089/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:1258613/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:129774/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:1410935/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:1438553/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:1539144/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:1691464/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:1695017/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:1740407/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:17487/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:1759219/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:1964011/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:2217352/43‑1 (MQ=38)
                  gACCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:2518783/43‑1 (MQ=38)
                   aCCATATAGCGTACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:554074/42‑1 (MQ=37)
                   aCCATATAGCGCACATGCAGCTTCGCCATTTCGTCGTCGgag                                   <  1:696329/42‑1 (MQ=38)
                                cGTGTAGCTTCGCCATTTCGTCGTCGGAGAGATCAACCACGCTGTAGCCGCCGTCGTTGAGc  >  1:2663082/1‑62 (MQ=255)
                                cGTGTAGCTTCGCCATTTCGTCGTCGGAGAGATCAACCACGCTGTAGCCGCCGTCGTTGAGc  >  1:2453878/1‑62 (MQ=255)
                                cGTGTAGCTTCGCCATTTCGTCGTCGGAGAGATCAACCACGATGTAGCCGCCGTCGTTGAGc  >  1:1174695/1‑62 (MQ=255)
                                    |                                                         
ATGCGATGGACGTCCGCCGACCATATAGCGCACGTGTAGCTTCGCCATTTCGTCGTCGGAGAGATCAACCACGCTGTAGCCGCCGTCGTTGAGC  >  W3110S.gb/3680064‑3680157

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: