Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,263,6370GA55.6% 3.9 / 31.9 27A273V (GCT→GTT) ispE4‑diphosphocytidyl‑2‑C‑methylerythritol kinase
Reads supporting (aligned to +/- strand):  ref base G (1/11);  new base A (0/15);  total (1/26)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCATTGAGCCATTCCGGGGCTTGCTCT  >  W3110S.gb/1263584‑1263681
                                                     |                                            
gAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAACGCCTTTc                                      <  1:1241733/62‑1 (MQ=255)
gAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAACGCCTTTc                                      <  1:897990/62‑1 (MQ=255)
gAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAACGCCTTTc                                      <  1:529234/62‑1 (MQ=255)
gAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAACGCCTTTc                                      <  1:4640286/62‑1 (MQ=255)
gAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAACGCCTTTc                                      <  1:3672142/62‑1 (MQ=255)
gAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAACGCCTTTc                                      <  1:3619632/62‑1 (MQ=255)
gAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAACGCCTTTc                                      <  1:3279335/62‑1 (MQ=255)
gAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAACGCCTTTc                                      <  1:2484940/62‑1 (MQ=255)
gAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAACGCCTTTc                                      <  1:217610/62‑1 (MQ=255)
gAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAACGCCTTTc                                      <  1:131284/62‑1 (MQ=255)
gAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAACGCCTTTc                                      <  1:1507792/62‑1 (MQ=255)
gAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAACGCCTTTc                                      <  1:1916286/62‑1 (MQ=255)
gAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAACGCCTTTc                                      <  1:152797/62‑1 (MQ=255)
gAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAACGCCTTTc                                      <  1:1520437/62‑1 (MQ=255)
    cTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCca                                  <  1:2436763/62‑1 (MQ=255)
    cTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCca                                  <  1:3112131/62‑1 (MQ=255)
    cTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCca                                  <  1:4610079/62‑1 (MQ=255)
    cTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCca                                  <  1:2308467/62‑1 (MQ=255)
    cTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCca                                  <  1:369100/62‑1 (MQ=255)
     tCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAACGCCTTTc                                      <  1:3648070/57‑1 (MQ=255)
         cTTGCCCGGCTTATAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAg                             <  1:3437850/62‑1 (MQ=255)
         cTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAg                             <  1:594044/62‑1 (MQ=255)
           tGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGcc                           <  1:4190597/62‑1 (MQ=255)
           tGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGcc                           <  1:1101386/62‑1 (MQ=255)
                             gCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCATTGAGCCATTCCGGGGc         <  1:302478/62‑1 (MQ=255)
                             gCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCATTGAGCCATTCCGGGGc         <  1:726576/62‑1 (MQ=255)
                                    gCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCATTGAGCCATTCCGGGGCTTGttct  >  1:3264518/1‑62 (MQ=255)
                                                     |                                            
GAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCATTGAGCCATTCCGGGGCTTGCTCT  >  W3110S.gb/1263584‑1263681

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: