Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,622,2930TA51.7% ‑5.7 / 42.0 29I187F (ATT→TTT) eamAcysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base T (7/7);  new base A (15/0);  total (22/7)
Fisher's exact test for biased strand distribution p-value = 2.20e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGA  >  W3110S.gb/1622239‑1622353
                                                      |                                                            
tgtgAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAAAGag                                                           >  1:18948/1‑57 (MQ=255)
tgtgAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAAAGag                                                           >  1:3782042/1‑57 (MQ=255)
tgtgAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAAAGag                                                           >  1:1103348/1‑57 (MQ=255)
tgtgAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAAAGag                                                           >  1:4510999/1‑57 (MQ=255)
tgtgAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAAAGag                                                           >  1:2633532/1‑57 (MQ=255)
tgtgAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAAAGag                                                           >  1:585457/1‑57 (MQ=255)
tgtgAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAAAGag                                                           >  1:588851/1‑57 (MQ=255)
tgtgAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAAAGag                                                           >  1:642584/1‑57 (MQ=255)
tgtgAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAAAGag                                                           >  1:1894653/1‑57 (MQ=255)
tgtgAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAAAGag                                                           >  1:649939/1‑57 (MQ=255)
tgtgAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAAAGag                                                           >  1:1607551/1‑57 (MQ=255)
tgtgAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAAAGag                                                           >  1:1499838/1‑57 (MQ=255)
tgtgAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAAAGag                                                           >  1:705374/1‑57 (MQ=255)
tgtgAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAAAGag                                                           >  1:1222422/1‑57 (MQ=255)
tgtgAATCATGGTTGCGGAACCATAGAGAATCAGCGAGGCAACAAAGAAGGGAAAGag                                                           >  1:1704998/1‑57 (MQ=255)
         tGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAgcg                                              >  1:3452074/1‑62 (MQ=255)
         tGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAgcg                                              >  1:3905533/1‑62 (MQ=255)
         tGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAgcg                                              >  1:211913/1‑62 (MQ=255)
         tGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAgcg                                              >  1:1640128/1‑62 (MQ=255)
               cGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGAt                                                     >  1:2692796/1‑49 (MQ=255)
               cGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAg                                        >  1:4265186/1‑62 (MQ=255)
                   aCCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTa                                    <  1:526886/62‑1 (MQ=255)
                   aCCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTa                                    <  1:245523/62‑1 (MQ=255)
                   aCCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTa                                    <  1:92827/62‑1 (MQ=255)
                      aTCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCa                                 <  1:3569824/62‑1 (MQ=255)
                      aTCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCa                                 <  1:4320890/62‑1 (MQ=255)
                      aTCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCa                                 <  1:1316831/62‑1 (MQ=255)
                                aGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCAc                       >  1:948874/1‑62 (MQ=255)
                                                     aTGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGa  <  1:2416915/62‑1 (MQ=255)
                                                      |                                                            
TGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGA  >  W3110S.gb/1622239‑1622353

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