Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 593,454 T→G I180I (ATA→ATC cusS ← sensory histidine kinase in two‑component regulatory system with CusR, senses copper ions

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb593,4540TG84.1% 197.2 / 19.2 69I180I (ATA→ATCcusSsensory histidine kinase in two‑component regulatory system with CusR, senses copper ions
Reads supporting (aligned to +/- strand):  ref base T (0/11);  new base G (0/58);  total (0/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTAGA  >  W3110S.gb/593393‑593494
                                                             |                                        
tAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTTTTCATCAAATCAttt                                          >  1:1712588/1‑62 (MQ=255)
tAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCAttt                                          >  1:229596/1‑62 (MQ=255)
tAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCAttt                                          >  1:65170/1‑62 (MQ=255)
tAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCAttt                                          >  1:2230158/1‑62 (MQ=255)
tAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCAttt                                          >  1:537214/1‑62 (MQ=255)
tAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCAttt                                          >  1:2226817/1‑62 (MQ=255)
tAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCAttt                                          >  1:1671526/1‑62 (MQ=255)
tAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCAttt                                          >  1:495330/1‑62 (MQ=255)
tAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCAttt                                          >  1:2530516/1‑62 (MQ=255)
tAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCAttt                                          >  1:2393666/1‑62 (MQ=255)
tAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCAttt                                          >  1:1925112/1‑62 (MQ=255)
tAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCAttt                                          >  1:360746/1‑62 (MQ=255)
tAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCAttt                                          >  1:2380762/1‑62 (MQ=255)
tAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCActt                                          >  1:360302/1‑62 (MQ=255)
                   gATTACCGATGCTGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1943379/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGCTGAAAat                       <  1:2725491/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGCTGAAAat                       <  1:2722244/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:2368961/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:2290767/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:2368384/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:2388323/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:2361816/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:23349/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:2325845/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:2418126/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:2266647/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:2495209/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:2498566/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:2752529/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:2795190/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:2830043/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:2868367/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:2895110/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:314551/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:372697/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:467226/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:483942/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:520677/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:53158/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:597847/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:853858/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:900978/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1674198/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1385407/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1381247/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1485407/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1496220/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1273054/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:152168/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1552643/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1229819/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1599994/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1184372/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1183317/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1647269/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1173403/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1150246/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1833396/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1888210/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1950323/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:200291/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:2007545/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1113518/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:2118153/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:2149125/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:2261948/62‑1 (MQ=255)
                   gATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1051799/62‑1 (MQ=255)
                    aTTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:22915/61‑1 (MQ=255)
                    aTTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1390731/61‑1 (MQ=255)
                            tGCGGTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:342059/53‑1 (MQ=255)
                               ggTCATAATAAGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:1068745/50‑1 (MQ=255)
                                    taataaGTTTATTCATCAAATCATTGATGTAATGAAGATGAAAat                       <  1:2627022/45‑1 (MQ=255)
                                        aaGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTaga  <  1:1431736/62‑1 (MQ=255)
                                        aaGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTaga  <  1:1296660/62‑1 (MQ=255)
                                        aaGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTaga  <  1:1514743/62‑1 (MQ=255)
                                        aaGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTaga  <  1:484814/62‑1 (MQ=255)
                                        aaGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTaga  <  1:1554102/62‑1 (MQ=255)
                                        aaGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTaga  <  1:519554/62‑1 (MQ=255)
                                        aaGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTaga  <  1:2711332/62‑1 (MQ=255)
                                        aaGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTaga  <  1:1772964/62‑1 (MQ=255)
                                        aaGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTaga  <  1:1156684/62‑1 (MQ=255)
                                        aaGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTaga  <  1:1821045/62‑1 (MQ=255)
                                        aaGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTaga  <  1:2026870/62‑1 (MQ=255)
                                                             |                                        
TAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTAGA  >  W3110S.gb/593393‑593494

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: