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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,122,1460AG13.2% 128.6 / 11.1 53I50I (ATT→ATCwzxCcolanic acid exporter
Reads supporting (aligned to +/- strand):  ref base A (24/22);  new base G (4/3);  total (28/25)
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.94e-01

ACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGC  >  W3110S.gb/2122077‑2122215
                                                                     |                                                                     
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:2141850/69‑1 (MQ=255)
                  aataatCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGt                                                      <  1:1847007/69‑1 (MQ=255)
                   ataatCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAg                                                   <  1:2950139/71‑1 (MQ=255)
                   ataatCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAg                                                   <  1:2640425/71‑1 (MQ=255)
                   ataatCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAg                                                   >  1:1382862/1‑71 (MQ=255)
                     aatGGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCAATGATc                                                                    >  1:152949/1‑52 (MQ=38)
                     aatGGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCAATGATc                                                                    >  1:1486312/1‑52 (MQ=38)
                     aatGGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCAATGATc                                                                    >  1:3657794/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:311692/71‑1 (MQ=255)
                          aGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCg                                            <  1:2741337/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:150599/1‑71 (MQ=255)
                                               aaaGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATc                       >  1:2128846/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:2641477/51‑1 (MQ=255)
                                                              gcgcGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGa                            <  1:1378472/51‑1 (MQ=255)
                                                              gcgcGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGa                            <  1:92666/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)
                                                                   aTAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTc                              <  1:1376958/44‑1 (MQ=255)
                                                                    tAATCACCAGCGACACGGTAAGCATGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGc  <  1:2428753/71‑1 (MQ=255)
                                                                    tAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGc  <  1:3261619/71‑1 (MQ=255)
                                                                     |                                                                     
ACTGATTTCTTTTCGCTGAATAATCGAGTTAGCGATACCGAAGTCAGAAAGCGTATCTGCCAGCGCGATAATCACCAGCGACACGGTAAGCAGGCCGAACTGGTGGTTGTCGATAATCCGCGCCAGCACGGTCATCTGC  >  W3110S.gb/2122077‑2122215

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: