Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb486,7970TA72.5% 105.2 / 51.4 69S346S (TCT→TCAkefAfused mechanosensitive channel proteins
Reads supporting (aligned to +/- strand):  ref base T (12/7);  new base A (50/0);  total (62/7)
Fisher's exact test for biased strand distribution p-value = 4.67e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGA  >  W3110S.gb/486761‑486855
                                    |                                                          
gCAGATTGCCTTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                   >  1:32144/1‑62 (MQ=255)
gCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                   >  1:1005107/1‑62 (MQ=255)
gCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                   >  1:2926466/1‑62 (MQ=255)
gCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                   >  1:547985/1‑62 (MQ=255)
gCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                   >  1:585375/1‑62 (MQ=255)
gCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                   >  1:1663614/1‑62 (MQ=255)
gCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                   >  1:1587723/1‑62 (MQ=255)
gCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                   >  1:1354128/1‑62 (MQ=255)
gCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                   >  1:1289050/1‑62 (MQ=255)
gCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                   >  1:809879/1‑62 (MQ=255)
gCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                   >  1:827800/1‑62 (MQ=255)
gCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                   >  1:2290987/1‑62 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACaaa                                   >  1:600309/1‑37 (MQ=37)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGa                >  1:1855020/1‑56 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2425898/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2385515/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2601543/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2278927/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2724623/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2759978/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2884081/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2188077/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2205327/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:364425/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:386113/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:392385/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:420020/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:503862/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:508112/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:542518/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:587652/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:665059/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:845169/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:904239/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:911840/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:948094/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1997249/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1871057/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1882546/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1798240/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1894827/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1910216/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1927762/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1937820/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1946007/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:195024/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1963184/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1965338/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1227499/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1217898/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1213648/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2089688/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2095060/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1188783/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2161700/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2175270/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1184895/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2190367/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTg            >  1:976194/1‑60 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTg            >  1:1867712/1‑60 (MQ=255)
                         ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTAGGCAGATGAACTgg           >  1:540106/1‑61 (MQ=255)
                         ctgctgTTGTCACGTATCCTCTACCAGCAACAACAAACACTCCCCTCGGCAGATGAAc              >  1:20110/1‑58 (MQ=255)
                                 gtCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGa  <  1:1564645/62‑1 (MQ=255)
                                 gtCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGa  <  1:1556384/62‑1 (MQ=255)
                                 gtCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGa  <  1:69363/62‑1 (MQ=255)
                                 gtCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGa  <  1:2787516/62‑1 (MQ=255)
                                 gtCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGa  <  1:21047/62‑1 (MQ=255)
                                 gtCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGa  <  1:2187547/62‑1 (MQ=255)
                                 gtCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGa  <  1:100084/62‑1 (MQ=255)
                                  tctcGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGa  <  1:1261116/61‑1 (MQ=255)
                                  tctcGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGa  <  1:2742204/61‑1 (MQ=255)
                                    |                                                          
GCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGA  >  W3110S.gb/486761‑486855

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