Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,365,0570AT65.1% 69.4 / 69.3 66pseudogene (18/480 nt)
pseudogene (18/1745 nt)
yhcE
yhcE
ECK3207:JW3184:b3217; hypothetical protein, N‑ter fragment
ECK3207:JW3184+JW3187:b4569; hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (23/0);  new base T (0/43);  total (23/43)
Fisher's exact test for biased strand distribution p-value = 2.87e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTGCCAATGAAACGGATAATAACAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATATCTCCAGTAGCAT  >  W3110S.gb/3365034‑3365119
                       |                                                              
gTGCCAATGACACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:371183/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCTTGGCCGCTCAGGCaga                       <  1:1228946/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCATGCaga                       <  1:1580099/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCtga                       <  1:326051/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:2217631/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:600986/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:588899/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:566323/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:2069604/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:2146331/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:2165213/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:50880/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:481028/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:2217266/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:646407/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:401262/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:2275849/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:2294453/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:2327520/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:2366485/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:2458342/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:14600/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:1066164/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:1116766/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:999886/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:983557/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:1265667/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:1297576/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:1305421/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:1319973/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:915994/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:1864889/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:1475240/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:1495315/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:1546771/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:850302/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:655035/65‑1 (MQ=255)
gTGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:1737330/65‑1 (MQ=255)
 tGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:2042064/64‑1 (MQ=255)
 tGCCAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:1764272/64‑1 (MQ=255)
   ccAATGAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:1082701/62‑1 (MQ=255)
        gAAACGGATAATAACTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:558768/57‑1 (MQ=255)
            cGGATAATAACAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATAtc             >  1:1993120/1‑63 (MQ=255)
            cGGATAATAACAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATAtc             >  1:977650/1‑63 (MQ=255)
            cGGATAATAACAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATAtc             >  1:1348464/1‑63 (MQ=255)
            cGGATAATAACAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATAtc             >  1:852195/1‑63 (MQ=255)
            cGGATAATAACAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATAtc             >  1:1611809/1‑63 (MQ=255)
            cGGATAATAACAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATAtc             >  1:83310/1‑63 (MQ=255)
            cGGATAATAACAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATAtc             >  1:780360/1‑63 (MQ=255)
            cGGATAATAACAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATAtc             >  1:1652955/1‑63 (MQ=255)
            cGGATAATAACAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATAtc             >  1:1785796/1‑63 (MQ=255)
            cGGATAATAACAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATAtc             >  1:631771/1‑63 (MQ=255)
            cGGATAATAACAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATAtc             >  1:61409/1‑63 (MQ=255)
            cGGATAATAACAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATAtc             >  1:1174054/1‑63 (MQ=255)
            cGGATAATAACAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATAtc             >  1:2168533/1‑63 (MQ=255)
            cGGATAATAACAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATAtc             >  1:2193331/1‑63 (MQ=255)
            cGGATAATAACAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATAtc             >  1:430704/1‑63 (MQ=255)
            cGGATAATAACAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATAtc             >  1:2251771/1‑63 (MQ=255)
            cGGATAATAACAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATAtc             >  1:2367779/1‑63 (MQ=255)
            cGGATAATAACAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATAtc             >  1:2390172/1‑63 (MQ=255)
                   taaCTGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCaga                       <  1:368368/46‑1 (MQ=255)
                      cAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAgtgt                   >  1:503084/1‑47 (MQ=255)
                      cAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATATCTCCAGTAGCAt  >  1:117066/1‑64 (MQ=255)
                      cAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATATCTCCAGTAGCAt  >  1:308688/1‑64 (MQ=255)
                      cAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATATCTCCAGTAGCAt  >  1:1939886/1‑64 (MQ=255)
                      cAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGACTGTaa                 >  1:200517/1‑49 (MQ=255)
                       |                                                              
GTGCCAATGAAACGGATAATAACAGGATGCCTGCTGCTGAACTTCGCCATGGCCGCTCAGGCAGAGTGTAATATCTCCAGTAGCAT  >  W3110S.gb/3365034‑3365119

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: