Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,679,4310GT72.6% 62.7 / 25.2 62R657R (CGC→CGAyphGconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/17);  new base T (45/0);  total (45/17)
Fisher's exact test for biased strand distribution p-value = 1.35e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTT  >  W3110S.gb/2679379‑2679475
                                                    |                                            
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:66404/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:50515/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:51813/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:53606/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:54718/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:56136/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:57261/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:61619/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:61812/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:62669/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:63298/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:64455/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:64776/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:66042/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:15685/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:67836/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:67868/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:68918/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:69563/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:73150/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:76434/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:76447/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:81134/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:81479/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:82566/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:85670/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:87113/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:87896/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:1535/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:47852/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:18893/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:20975/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:25630/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:457/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:26835/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:28110/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:44833/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTgc                                     >  1:31044/1‑62 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTg                                      >  1:34570/1‑61 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTg                                      >  1:9618/1‑61 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTg                                      >  1:1870/1‑61 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTg                                      >  1:34413/1‑61 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTg                                      >  1:35086/1‑61 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTg                                      >  1:47740/1‑61 (MQ=255)
tCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACAAAGACCTCGCCACTg                                      >  1:55353/1‑61 (MQ=255)
                                                gACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGtt  <  1:4501/49‑1 (MQ=255)
                                                gACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGtt  <  1:50049/49‑1 (MQ=255)
                                                gACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGtt  <  1:50646/49‑1 (MQ=255)
                                                gACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGtt  <  1:47788/49‑1 (MQ=255)
                                                gACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGtt  <  1:20157/49‑1 (MQ=255)
                                                gACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGtt  <  1:46054/49‑1 (MQ=255)
                                                gACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGtt  <  1:79596/49‑1 (MQ=255)
                                                gACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGtt  <  1:77803/49‑1 (MQ=255)
                                                gACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGtt  <  1:35510/49‑1 (MQ=255)
                                                gACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGtt  <  1:61773/49‑1 (MQ=255)
                                                gACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGtt  <  1:29123/49‑1 (MQ=255)
                                                gACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGtt  <  1:40307/49‑1 (MQ=255)
                                                gACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGtt  <  1:33040/49‑1 (MQ=255)
                                                gACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGtt  <  1:33387/49‑1 (MQ=255)
                                                gACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGtt  <  1:37889/49‑1 (MQ=255)
                                                gACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGtt  <  1:50326/49‑1 (MQ=255)
                                                gACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGtt  <  1:34893/49‑1 (MQ=255)
                                                    |                                            
TCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTT  >  W3110S.gb/2679379‑2679475

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: