Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,389,530 C→T S52S (AGC→AGT aaeR → predicted DNA‑binding transcriptional regulator, efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,389,5300CT88.4% 112.2 / 6.2 43S52S (AGC→AGTaaeRpredicted DNA‑binding transcriptional regulator, efflux system
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base T (38/0);  total (41/2)
Fisher's exact test for biased strand distribution p-value = 1.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

TCAGTCAGACGGTATCAAAACTGGAAGATGAGTTGCAGGTAAAGCTGTTAAACCGTAGCACACGCAGCATTGGCCTGACCGAAGCCGG  >  W3110S.gb/3389472‑3389559
                                                          |                             
tcggtcagACGGTATCAAAACTGGAAGATGAGTTGCAGGTAAAGCTGTTTAACCGTAGcac                             <  1:980637/58‑1 (MQ=255)
tcagtcagACGGTATCAAAACTGGAAGATGAGTTGCAGGTAAAGCTGTTAAACCGTAGcac                             <  1:2431595/61‑1 (MQ=255)
tcagtcagACGGTATCAAAACTGGAAGATGAGTTGCAGGTAAAGCTGTTAAACCGTAGcac                             <  1:3126235/61‑1 (MQ=255)
tcagtcagACGGTATCAAAACTGGAAGATGAGTTGCAGGTAAAGCTGTTAAACCGTAGcac                             <  1:2504719/61‑1 (MQ=255)
tcagtcagACGGTATCAAAACTGGAAGATGAGTTGCAGGTAAAGCTGTTAAACCGTAGcac                             <  1:1703125/61‑1 (MQ=255)
tcagtcagACGGTATCAAAACTGGAAGATGAGTTGCAGGTAAAGCTGTTAAACCGTAGcac                             <  1:2825669/61‑1 (MQ=255)
       gACGGTATCAAAACTGGAAGATGAGTTGCAGGTAAAGCTGTTAAACCGTAGCACACgcagca                     >  1:1683297/1‑62 (MQ=255)
       gACGGTATCAAAACTGGAAGATGAGTTGCAGGTAAAGCTGTTAAACCGTAGCACACgcagca                     >  1:1821296/1‑62 (MQ=255)
       gACGGTATCAAAACTGGAAGATGAGTTGCAGGTAAAGCTGTTAAACCGTAGCACACgcagca                     >  1:67752/1‑62 (MQ=255)
                      ggAAGATGAGTTGCAGGTAAAGCTGTTAAACCGTAGCACACGCAGCATTGGCCTGACCGAAg      <  1:949998/62‑1 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACGCGCAGCATTg                  >  1:1488433/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCAtt                   >  1:2150370/1‑48 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCAtag                  >  1:2127271/1‑47 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:2943336/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:3032542/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:3045633/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:3111358/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:115392/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:2919326/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:37553/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:434518/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:508476/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:559609/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:650524/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:695286/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:761491/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:894010/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:966638/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:978119/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:1974825/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:1052943/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:1158160/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:1614835/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:1631313/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:1684449/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:1746957/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:1761911/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:1882073/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:285164/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:1999287/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:2139493/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:2152367/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:218770/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:2478731/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:255773/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:265561/1‑49 (MQ=255)
                       gAAGATGAGTTACAGGTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:2812731/1‑49 (MQ=255)
                       gAAGATCAGTTACAGCTAAAGCTGTTAAACCGTAGTACACGCAGCATTg                  >  1:97114/1‑49 (MQ=255)
                          gATGAGTTGCAGGTAAAGCTGTTAAACCGTAGCACACGCAGCATTGGCCTGACCGAAGCCgg  <  1:2623298/62‑1 (MQ=255)
                                                          |                             
TCAGTCAGACGGTATCAAAACTGGAAGATGAGTTGCAGGTAAAGCTGTTAAACCGTAGCACACGCAGCATTGGCCTGACCGAAGCCGG  >  W3110S.gb/3389472‑3389559

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