Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 1,995,272 T→C S524S (TCA→TCG uvrC ← excinuclease UvrABC, endonuclease subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,995,2720TC98.6% 245.5 / ‑4.2 69S524S (TCA→TCGuvrCexcinuclease UvrABC, endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (40/28);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

ACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTAC  >  W3110S.gb/1995206‑1995339
                                                                  |                                                                   
aCGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGaa                                                                  <  1:3460796/70‑1 (MQ=255)
      cccGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGCGAATCTggc                                                            <  1:1895154/70‑1 (MQ=255)
      cccGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGCGAATCTggc                                                            <  1:2797983/70‑1 (MQ=255)
                  gTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGCGaa                                                                  >  1:1519316/1‑52 (MQ=255)
                  gTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATc                                               >  1:2796160/1‑71 (MQ=255)
                  gTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATc                                               >  1:1396763/1‑71 (MQ=255)
                  gTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATc                                               >  1:704875/1‑71 (MQ=255)
                           tgATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCtt                                      >  1:1726018/1‑71 (MQ=255)
                           tgATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCtt                                      >  1:3677871/1‑71 (MQ=255)
                           tgATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCtt                                      >  1:3976495/1‑71 (MQ=255)
                             aTTCATCGCGAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTcc                                    >  1:1739506/1‑71 (MQ=255)
                             aTTCATCGCGAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTc                                     >  1:1161169/1‑70 (MQ=255)
                                      gAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATc                                               >  1:3395305/1‑51 (MQ=255)
                                      gAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATc                                               >  1:1376740/1‑51 (MQ=255)
                                      gAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATc                                               >  1:3184129/1‑51 (MQ=255)
                                      gAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCtt                                      >  1:3161071/1‑60 (MQ=255)
                                      gAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCtt                                      >  1:3296626/1‑60 (MQ=255)
                                      gAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCtt                                      >  1:1859358/1‑60 (MQ=255)
                                      gAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAg                           <  1:2386343/71‑1 (MQ=255)
                                      gAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAg                           <  1:795391/71‑1 (MQ=255)
                                      gAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAg                           <  1:2510734/71‑1 (MQ=255)
                                      gAATATGCTGGATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAg                           <  1:1364761/71‑1 (MQ=255)
                                         tatGCTGTATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAg                           <  1:1824313/68‑1 (MQ=255)
                                               ggATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGtt                  >  1:3162168/1‑71 (MQ=255)
                                               ggATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGtt                  >  1:569247/1‑71 (MQ=255)
                                                gATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGtt                  <  1:617067/70‑1 (MQ=255)
                                                gATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGtt                  <  1:2721017/70‑1 (MQ=255)
                                                gATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGtt                  <  1:2416099/70‑1 (MQ=255)
                                                gATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGtt                  <  1:356494/70‑1 (MQ=255)
                                                gATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGtt                  <  1:3583549/70‑1 (MQ=255)
                                                gATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGtt                  <  1:3620564/70‑1 (MQ=255)
                                                gATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGtt                  <  1:3886652/70‑1 (MQ=255)
                                                gATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGtt                  <  1:904355/70‑1 (MQ=255)
                                                gATAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGtt                  <  1:3839117/70‑1 (MQ=255)
                                                  tAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtc                >  1:922428/1‑70 (MQ=255)
                                                  tAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtc                >  1:2581769/1‑70 (MQ=255)
                                                  tAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtc                >  1:845316/1‑70 (MQ=255)
                                                  tAACATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtc                >  1:2255935/1‑70 (MQ=255)
                                                    aCATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGttta                >  1:1976630/1‑67 (MQ=255)
                                                    aCATGCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtc                >  1:2480811/1‑68 (MQ=255)
                                                        gCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGtt                  >  1:4082800/1‑62 (MQ=255)
                                                        gCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGtt                  >  1:3062720/1‑62 (MQ=255)
                                                        gCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtccagtc           >  1:3148894/1‑69 (MQ=255)
                                                        gCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtccagtc           >  1:2348827/1‑69 (MQ=255)
                                                        gCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtccagtc           >  1:2076821/1‑69 (MQ=255)
                                                        gCAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtccagtc           >  1:1722489/1‑69 (MQ=255)
                                                        gCAGCGAGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtccagtc           >  1:5621/1‑69 (MQ=255)
                                                         cAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtcca              >  1:2248382/1‑65 (MQ=255)
                                                         cAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtcca              >  1:1709822/1‑65 (MQ=255)
                                                         cAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtcca              >  1:231827/1‑65 (MQ=255)
                                                         cAGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtcca              >  1:2864977/1‑65 (MQ=255)
                                                          aGCGCGGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGc       <  1:1558867/71‑1 (MQ=255)
                                                             gcgGGCTAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTAAAAGAACAGCGTTTCCAGTCCAGCCtt    <  1:3802766/71‑1 (MQ=255)
                                                             gcgGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCtt    <  1:880130/71‑1 (MQ=255)
                                                             gcgGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCtt    <  1:1382140/71‑1 (MQ=255)
                                                             gcgGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTAAAAGAACAGCGTTTCCAGTCCAGCCtt    <  1:1058528/71‑1 (MQ=255)
                                                             gcgGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTAAAAGAACAGCGTTTCCAGTCCAGCCtt    <  1:3733678/71‑1 (MQ=255)
                                                             gcgGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTAAAAGAACAGCGTTTCCAGTCCAGCCtt    <  1:2792076/71‑1 (MQ=255)
                                                             gcgGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTAAAAGAACAGCGTTTCCAGTCCAGCCtt    <  1:2531389/71‑1 (MQ=255)
                                                             gcgGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTAAAAGAACAGCGTTTCCAGTCCAGCCtt    <  1:1457946/71‑1 (MQ=255)
                                                              cgGGCGAATCTGGCGGCAAACTAAATCCCTCACCtt                                      >  1:2686891/1‑36 (MQ=255)
                                                              cgGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtccagtc           >  1:32317/1‑63 (MQ=255)
                                                              cgGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtccagtc           >  1:2601467/1‑63 (MQ=255)
                                                              cgGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtccagtc           >  1:2459122/1‑63 (MQ=255)
                                                              cgGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtccagtc           >  1:988207/1‑63 (MQ=255)
                                                              cgGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtccagt            >  1:1263049/1‑62 (MQ=255)
                                                              cgGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCtt    <  1:2812295/70‑1 (MQ=255)
                                                              cgGGCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCtt    <  1:1014484/70‑1 (MQ=255)
                                                               gggCGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTAc  <  1:1655103/71‑1 (MQ=255)
                                                                  |                                                                   
ACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTAC  >  W3110S.gb/1995206‑1995339

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