Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,122,101 C→G 11.5% S65S (TCG→TCC wzxC ← colanic acid exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,122,1010CG11.5% 152.7 / 9.6 61S65S (TCG→TCCwzxCcolanic acid exporter
Reads supporting (aligned to +/- strand):  ref base C (29/24);  new base G (4/3);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

TCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAG  >  W3110S.gb/2122031‑2122169
                                                                      |                                                                    
tcgccAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATc                                                                      <  1:142436/68‑1 (MQ=255)
tCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATc                                                                      <  1:65065/71‑1 (MQ=255)
  cccAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATc                                                                      >  1:2280276/1‑69 (MQ=255)
  cccAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCg                                                                     <  1:1123963/70‑1 (MQ=255)
  cccAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGa                                                                    <  1:26886/71‑1 (MQ=255)
  cccAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGa                                                                    <  1:3493484/71‑1 (MQ=255)
    cAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGt                                                                  >  1:1847693/1‑71 (MQ=255)
    cAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGt                                                                  >  1:2539589/1‑71 (MQ=255)
    cAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGt                                                                  <  1:1190425/71‑1 (MQ=255)
        ccGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGt                                                                  <  1:845200/67‑1 (MQ=255)
        ccGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGc                                                              <  1:2027896/71‑1 (MQ=255)
           aCGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTATCGCTGAATAATCGAGATAGCGa                                                            >  1:2331842/1‑70 (MQ=255)
                cAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCg                                                       >  1:3225178/1‑70 (MQ=255)
                   ccAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGt                                                   <  1:148691/71‑1 (MQ=255)
                      gTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGt                                                   <  1:3328990/68‑1 (MQ=255)
                      gTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAg                                                <  1:1293850/71‑1 (MQ=255)
                         cAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGa                                                            >  1:1673860/1‑56 (MQ=255)
                         cAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGa                                               >  1:638354/1‑69 (MQ=255)
                         cAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGa                                               >  1:3259939/1‑69 (MQ=255)
                         cAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGa                                               >  1:3634674/1‑69 (MQ=255)
                           aCGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGt                                                   >  1:418606/1‑63 (MQ=255)
                           aCGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGt                                                   >  1:1837616/1‑63 (MQ=255)
                             gtggtgAGTTCAAGGTGAGTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGc                                           >  1:3034094/1‑69 (MQ=255)
                             gtggtgAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGc                                                              >  1:2806837/1‑50 (MQ=255)
                             gtggtgAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGt                                         <  1:184242/71‑1 (MQ=255)
                                 tgAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCt                                     <  1:2521729/71‑1 (MQ=255)
                                   aGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGc                                   >  1:710953/1‑71 (MQ=255)
                                     ttCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAg                                                                   <  1:2789549/37‑1 (MQ=255)
                                     ttCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAg                                                                   <  1:2112839/37‑1 (MQ=255)
                                     ttCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAg                                                                   <  1:872092/37‑1 (MQ=255)
                                     ttCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAg                                                                   <  1:3801762/37‑1 (MQ=255)
                                     ttCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGc                                   >  1:2793990/1‑69 (MQ=255)
                                       cAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGa                                                            >  1:2954240/1‑42 (MQ=255)
                                       cAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAgc                               >  1:1443742/1‑71 (MQ=255)
                                       cAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGAGATACCGAAGTCAGAAAGCGTATCTGCCAgc                               >  1:1943205/1‑71 (MQ=255)
                                         aGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGc                                                              <  1:164067/38‑1 (MQ=255)
                                          ggTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGcgcg                            <  1:3532694/71‑1 (MQ=255)
                                            tGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGa                           >  1:673010/1‑70 (MQ=255)
                                             gACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGaa                                                     >  1:2176708/1‑43 (MQ=255)
                                              aCTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATaa                        >  1:2924625/1‑71 (MQ=255)
                                                 gATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCa                     >  1:3233737/1‑71 (MQ=255)
                                                   tttCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTa                                        <  1:2043577/50‑1 (MQ=255)
                                                     tCTTTTCGCTGAATAATCGAGTTAGCGATACCGaa                                                     >  1:2492029/1‑35 (MQ=255)
                                                     tCTTTTCGCTGAATAATCGAGTTAGCGATACCGaa                                                     >  1:1957891/1‑35 (MQ=255)
                                                         ttCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTAt                                       >  1:2832664/1‑45 (MQ=255)
                                                         ttCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTAt                                       >  1:1264458/1‑45 (MQ=255)
                                                         ttCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATATTCACCAGCGa              >  1:27246/1‑70 (MQ=255)
                                                         ttCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGa              >  1:2581279/1‑70 (MQ=255)
                                                          tCGCTGAATAATAGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGaca            >  1:3199720/1‑71 (MQ=255)
                                                             cTGAATAATGGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCAATGATCACCAGCGacacag         <  1:197859/71‑3 (MQ=255)
                                                             cTGAATAATGGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCAATGATCACCAGCGacacag         <  1:977366/71‑3 (MQ=255)
                                                              tGAATAATGGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCAATGATCACCAGCGacacag         <  1:2686239/70‑3 (MQ=255)
                                                                aataatCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGt        <  1:2141850/69‑1 (MQ=255)
                                                                aataatCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGt        <  1:1847007/69‑1 (MQ=255)
                                                                 ataatCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAg     <  1:2640425/71‑1 (MQ=255)
                                                                 ataatCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAg     <  1:2950139/71‑1 (MQ=255)
                                                                 ataatCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAg     >  1:1382862/1‑71 (MQ=255)
                                                                   aatGGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCAATGATc                      >  1:3657794/1‑52 (MQ=38)
                                                                   aatGGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCAATGATc                      >  1:2500120/1‑52 (MQ=38)
                                                                   aatGGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCAATGATc                      >  1:152949/1‑52 (MQ=38)
                                                                   aatGGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCAATGATc                      >  1:1486312/1‑52 (MQ=38)
                                                                    atCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAg  >  1:745882/1‑71 (MQ=255)
                                                                    atCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAg  <  1:1152296/71‑1 (MQ=255)
                                                                     tCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCAc                    <  1:951094/52‑1 (MQ=255)
                                                                      |                                                                    
TCCCCAGCCCGACGTTCAGCCAGTACAACGTGGTGAGTTCAAGGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAG  >  W3110S.gb/2122031‑2122169

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: