Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,819,1310CG56.9% 19.5 / 65.3 51intergenic (‑142/‑116)gmk/ligBguanylate kinase/DNA ligase, NAD(+)‑dependent
Reads supporting (aligned to +/- strand):  ref base C (4/18);  new base G (29/0);  total (33/18)
Fisher's exact test for biased strand distribution p-value = 2.62e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGCGCAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCG  >  W3110S.gb/3819084‑3819192
                                               |                                                             
tGCGCAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGGAAt                                                   <  1:2278180/60‑1 (MQ=255)
tGCGCAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAAt                                                   <  1:1144271/60‑1 (MQ=255)
tGCGCAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAAt                                                   <  1:956637/60‑1 (MQ=255)
tGCGCAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAAt                                                   <  1:3492984/60‑1 (MQ=255)
tGCGCAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAAt                                                   <  1:288313/60‑1 (MQ=255)
tGCGCAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAAt                                                   <  1:2098365/60‑1 (MQ=255)
tGCGCAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAAt                                                   <  1:1918751/60‑1 (MQ=255)
tGCGCAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAAt                                                   <  1:1792576/60‑1 (MQ=255)
 gcgcAATGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTAACTGCAATaa                                                 <  1:988399/61‑1 (MQ=255)
 gcgcAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATaa                                                 <  1:3097497/61‑1 (MQ=255)
 gcgcAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATaa                                                 <  1:4197838/61‑1 (MQ=255)
 gcgcAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATaa                                                 <  1:2687438/61‑1 (MQ=255)
 gcgcAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATaa                                                 <  1:474770/61‑1 (MQ=255)
  cgcAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATAAAg                                               >  1:1456174/1‑62 (MQ=255)
    cAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATAAAGAc                                             >  1:3013843/1‑62 (MQ=255)
    cAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATAAAGAc                                             >  1:2926115/1‑62 (MQ=255)
    cAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATAAAGAc                                             >  1:4330535/1‑62 (MQ=255)
         ttgttgCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATaa                                                 <  1:335654/53‑1 (MQ=255)
                     gCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATACAGACGATTTTCAGGAAGCGt                             <  1:1089014/61‑1 (MQ=255)
                     gCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATAAAGACGATTTTCAGGAAGCGt                             <  1:1235007/61‑1 (MQ=255)
                     gCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATAAAGACGATTTTCAGGAAGCGt                             <  1:4366890/61‑1 (MQ=255)
                      cGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATAAAGACGATTTTCAGGAAGCGt                             <  1:2646656/60‑1 (MQ=255)
                       gTAGCTGCTTCCTCATTTGCCTGTGATTCACTGCAATAAAGACg                                            >  1:3508874/1‑44 (MQ=38)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACg                                            >  1:2622653/1‑44 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:553055/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:3105082/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:3537547/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:1189299/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:947686/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:3941448/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:4044024/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:84386/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:656492/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:4647466/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:4379188/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:487358/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:4388242/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:4482382/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:4523667/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:1685809/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:3209070/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:1665868/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:1752751/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:1951153/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:2999852/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:2962097/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:1194838/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:2821064/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:2071781/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:2241661/1‑45 (MQ=255)
                       gTAGCTGCTTCCGCATTTGCCTGTGATTCACTGCAATAAAGACGa                                           >  1:3237381/1‑45 (MQ=255)
                                               cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTc                         >  1:2305222/1‑39 (MQ=255)
                                               cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg  >  1:2718133/1‑62 (MQ=255)
                                               cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg  >  1:2225101/1‑62 (MQ=255)
                                               cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg  >  1:213247/1‑62 (MQ=255)
                                               cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg  >  1:2080887/1‑62 (MQ=255)
                                               cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg  >  1:2718149/1‑62 (MQ=255)
                                               cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg  >  1:1988692/1‑62 (MQ=255)
                                               cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg  >  1:3711458/1‑62 (MQ=255)
                                               cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg  >  1:152128/1‑62 (MQ=255)
                                               cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg  >  1:1594118/1‑62 (MQ=255)
                                               cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg  >  1:3639992/1‑62 (MQ=255)
                                               cATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCg  >  1:119440/1‑62 (MQ=255)
                                               |                                                             
TGCGCAAAGTTGTTGCGTTTTGCGTAGCTGCTTCCGCATTTGCCTGTCATTCACTGCAATAAAGACGATTTTCAGGAAGCGTCCTCGCAGACTCCGTTTTTGCCACTCG  >  W3110S.gb/3819084‑3819192

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: