Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb26,9970CA71.0% 69.5 / 65.5 62R241S (CGC→AGC) ispH1‑hydroxy‑2‑methyl‑2‑(E)‑butenyl 4‑diphosphate reductase, 4Fe‑4S protein
Reads supporting (aligned to +/- strand):  ref base C (16/2);  new base A (0/44);  total (16/46)
Fisher's exact test for biased strand distribution p-value = 5.60e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAACTCCAACCGTCTGGCGGAGCTGGCCCAGCGTATGGGCAAACGCGCGTTTTTGATTGACGATGCGAAAGACATCCAGGAAG  >  W3110S.gb/26954‑27036
                                           |                                       
caactccaacCGTCTGGCGGAGCTGGCCCAGCGTATGGGCAAAc                                         >  1:188593/1‑44 (MQ=255)
caactccaacCGTCTGGCGGAGCTGGCCCAGCGTATGGGCAAACGCGCGTTTTTGATTGACGa                      >  1:77487/1‑63 (MQ=255)
caactccaacCGTCTGGCGGAGCTGGCCCAGCGTATGGGCAAACGCGCGTTTTTGATTGACGAt                     >  1:2013814/1‑64 (MQ=255)
caactccaacCGTCTGGCGGAGCTGGCCCAGCGTATGGGCAAACGCGCGTTTTTGATTGACGAt                     >  1:475852/1‑64 (MQ=255)
caactccaacCGTCTGGCGGAGCTGGCCCAGCGTATGGGCAAACGCGCGTTTTTGATTGACGAt                     >  1:174172/1‑64 (MQ=255)
caactccaacCGTCTGGCGGAGCTGGCCCAGCGTATGGGCAAACGCGCGTTTTTGATTGACGAt                     >  1:1683995/1‑64 (MQ=255)
caactccaacCGTCTGGCGGAGCTGGCCCAGCGTATGGGCAAACGCGCGTTTTTGATTGACGAt                     >  1:695908/1‑64 (MQ=255)
caactccaacCGTCTGGCGGAGCTGGCCCAGCGTATGGGCAAACGCGCGTTTTTGATTGACGAt                     >  1:1991048/1‑64 (MQ=255)
caactccaacCGTCTGGCGGAGCTGGCCCAGCGTATGGGCAAACGCGCGTTTTTGATTGACGAt                     >  1:204356/1‑64 (MQ=255)
caactccaacCGTCTGGCGGAGCTGGCCCAGCGTATGGGCAAACGCGCGTTTTTGATTGACGAt                     >  1:2072176/1‑64 (MQ=255)
caactccaacCGTCTGGCGGAGCTGGCCCAGCGTATGGGCAAACGCGCGTTTTTGATTGACGAt                     >  1:146795/1‑64 (MQ=255)
caactccaacCGTCTGGCGGAGCTGGCCCAGCGTATGGGCAAACGCGCGTTTTTGATTGACGAt                     >  1:1044820/1‑64 (MQ=255)
caactccaacCGTCTGGCGGAGCTGGCCCAGCGTATGGGCAAACGCGCGTTTTTGATTGACGAt                     >  1:959358/1‑64 (MQ=255)
caactccaacCGTCTGGCGGAGCTGGCCCAGCGTATGGGCAAACGCGCGTTTTTGAATGACGAt                     >  1:829407/1‑64 (MQ=255)
   ctccaacCGTCTGGCGGAGCTGGCCCAGCGTATGGGCAAACTCGCGTTTTTGATTGACGATGc                   >  1:2525147/1‑63 (MQ=255)
   ctccaacCGTCTGGCGGAGCTGGCCCAGCGTATGGGCAAACGCGCGTTTTTGATTGACGATGc                   >  1:834648/1‑63 (MQ=255)
   ctccaacCGTCTGGCGGAGCTGGCCCAGCGTATGGGCAAACGCGCGTTTTTGATTGACGATGc                   >  1:1843118/1‑63 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGTCATCCAg      <  1:2447684/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:2403057/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:369956/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:403633/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:2477776/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:438932/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:2083504/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:468523/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:563377/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:594597/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:718956/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:781414/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:78269/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:801313/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:825212/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:843361/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:868826/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:906513/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:977370/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:1614618/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:1120994/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:116875/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:1203349/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:1311078/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:1311768/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:1324872/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:1340673/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:1375375/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:1497843/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:15376/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:1561088/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:2332998/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:1646872/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:1657191/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:180090/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:1827578/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:1891839/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:1915649/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:2233888/65‑1 (MQ=255)
              tGGCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:2266259/65‑1 (MQ=255)
                gCGGAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:65419/63‑1 (MQ=255)
                  ggAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:1561089/61‑1 (MQ=255)
                  ggAGCTGGCCCAGCGTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:2138992/61‑1 (MQ=255)
                   gAGCTGGCCCAGCGTATGGGCAAACGCGCGTTTTTGATTGACGATGCGAAAGACATCCAGGAag  <  1:1562091/64‑1 (MQ=255)
                   gAGCTGGCCCAGCGTATGGGCAAACGCGCGTTTTTGATTGACGATGCGAAAGACATCCAGGAag  <  1:1955945/64‑1 (MQ=255)
                                gTATGGGCAAAAGCGCGTTTTTGATTGACGATGCGAAAGACATCCAg      <  1:428522/47‑1 (MQ=255)
                                           |                                       
CAACTCCAACCGTCTGGCGGAGCTGGCCCAGCGTATGGGCAAACGCGCGTTTTTGATTGACGATGCGAAAGACATCCAGGAAG  >  W3110S.gb/26954‑27036

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: ‑