Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,122,143 G→A 14.3% I51I (ATC→ATT wzxC ← colanic acid exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,122,1430GA14.3% 117.5 / 10.7 49I51I (ATC→ATTwzxCcolanic acid exporter
Reads supporting (aligned to +/- strand):  ref base G (24/18);  new base A (4/3);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

GGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATC  >  W3110S.gb/2122073‑2122212
                                                                      |                                                                     
ggTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGcgcg                                                                       <  1:3532694/71‑1 (MQ=255)
  tGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGa                                                                      >  1:673010/1‑70 (MQ=255)
    aCTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATaa                                                                   >  1:2924625/1‑71 (MQ=255)
       gATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCa                                                                >  1:3233737/1‑71 (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:977366/71‑3 (MQ=255)
                   cTGAATAATGGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCAATGATCACCAGCGacacag                                                    <  1:197859/71‑3 (MQ=255)
                    tGAATAATGGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCAATGATCACCAGCGacacag                                                    <  1:2686239/70‑3 (MQ=255)
                      aataatCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGt                                                   <  1:1847007/69‑1 (MQ=255)
                      aataatCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGt                                                   <  1:2141850/69‑1 (MQ=255)
                       ataatCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAg                                                <  1:2950139/71‑1 (MQ=255)
                       ataatCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAg                                                >  1:1382862/1‑71 (MQ=255)
                       ataatCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAg                                                <  1:2640425/71‑1 (MQ=255)
                         aatGGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCAATGATc                                                                 >  1:152949/1‑52 (MQ=38)
                         aatGGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCAATGATc                                                                 >  1:3657794/1‑52 (MQ=38)
                         aatGGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCAATGATc                                                                 >  1:1486312/1‑52 (MQ=38)
                         aatGGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCAATGATc                                                                 >  1:2500120/1‑52 (MQ=38)
                          atCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAg                                             <  1:1152296/71‑1 (MQ=255)
                          atCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAg                                             >  1:745882/1‑71 (MQ=255)
                           tCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCAc                                                               <  1:951094/52‑1 (MQ=255)
                             gAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGcc                                          >  1:2243459/1‑71 (MQ=255)
                             gAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGcc                                          >  1:998497/1‑71 (MQ=255)
                              aGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCg                                         <  1:2741337/71‑1 (MQ=255)
                              aGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCg                                         <  1:311692/71‑1 (MQ=255)
                                  aGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACt                                     >  1:746782/1‑71 (MQ=255)
                                  aGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACt                                     >  1:3203799/1‑71 (MQ=255)
                                   gCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCa                                                                <  1:2106873/43‑1 (MQ=255)
                                   gCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACtg                                    >  1:2506016/1‑71 (MQ=255)
                                              gTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTTGTTGTCGa                         >  1:2934084/1‑71 (MQ=255)
                                              gTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGa                         <  1:3091656/71‑1 (MQ=255)
                                                 aGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATaa                      >  1:2175301/1‑71 (MQ=255)
                                                  gAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCATGCCg                                         <  1:228143/51‑1 (MQ=255)
                                                   aaaGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATc                    >  1:2128846/1‑71 (MQ=255)
                                                   aaaGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATc                    >  1:150599/1‑71 (MQ=255)
                                                    aaGCGTATCTGCCAGCGCGATAATCACCA‑‑GACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCcgc                 >  1:1071938/1‑71 (MQ=255)
                                                     aGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCcg                  >  1:3197547/1‑71 (MQ=255)
                                                       cGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCcgcg                >  1:954711/1‑71 (MQ=255)
                                                       cGTAACTGGCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACtgg                                   >  1:3306015/1‑52 (MQ=39)
                                                           tCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGa                                        <  1:3522085/43‑1 (MQ=255)
                                                           tCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGa                                        <  1:46704/43‑1 (MQ=255)
                                                           tCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAg            >  1:191276/1‑71 (MQ=255)
                                                             tGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCa          >  1:1842792/1‑71 (MQ=255)
                                                              gcaaGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCa          >  1:1479658/4‑70 (MQ=255)
                                                               ccAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAg            <  1:3618746/67‑1 (MQ=255)
                                                                  gcgcGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGa                         <  1:1378472/51‑1 (MQ=255)
                                                                  gcgcGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGa                         <  1:92666/51‑1 (MQ=255)
                                                                  gcgcGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGa                         <  1:2641477/51‑1 (MQ=255)
                                                                  gcgcGATAATCAACAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGa                         <  1:1820334/51‑1 (MQ=255)
                                                                    gcgATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCAt   >  1:943476/1‑71 (MQ=255)
                                                                     cgATAATCACCAGCGACACGGTAAGCAGGCCGAACtggtg                                 >  1:2818573/1‑40 (MQ=255)
                                                                      gATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATc  <  1:2118646/70‑1 (MQ=255)
                                                                      |                                                                     
GGTGACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATC  >  W3110S.gb/2122073‑2122212

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