Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,507,250 A→G 21.9% I182T (ATC→ACC)  ydcO ← predicted benzoate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,507,2500AG21.9% 165.4 / 45.6 82I182T (ATC→ACC) ydcOpredicted benzoate transporter
Reads supporting (aligned to +/- strand):  ref base A (44/20);  new base G (14/4);  total (58/24)
Fisher's exact test for biased strand distribution p-value = 5.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTT  >  W3110S.gb/1507181‑1507315
                                                                     |                                                                 
tCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGAt                                                                  <  1:343455/71‑1 (MQ=255)
         atatAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGGTGacgatcac                                                         <  1:3896686/71‑1 (MQ=255)
         atatAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGGTGacgatcac                                                         <  1:3646752/71‑1 (MQ=255)
         atatAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGacgatcac                                                         >  1:2539158/1‑71 (MQ=255)
         atatAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGacgatcac                                                         >  1:1278281/1‑71 (MQ=255)
         atatAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGacgatcac                                                         >  1:260303/1‑71 (MQ=255)
         atatAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGacgatcac                                                         >  1:3606781/1‑71 (MQ=255)
          tataAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGGTGacg                                                              <  1:1118899/65‑1 (MQ=255)
           ataAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGacgatcacga                                                       >  1:160476/1‑71 (MQ=255)
           ataAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGacgatcacga                                                       >  1:3005231/1‑71 (MQ=255)
           ataAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGacgatcacga                                                       >  1:1371220/1‑71 (MQ=255)
            taAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGacgatcacgat                                                      >  1:2435965/1‑71 (MQ=255)
                    gagaACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCaataa                                              <  1:3511245/71‑1 (MQ=255)
                    gagaACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCaataa                                              <  1:389204/71‑1 (MQ=255)
                    gagaACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCaataa                                              <  1:2370622/71‑1 (MQ=255)
                         cGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAatcatc                                         >  1:1789898/1‑71 (MQ=255)
                              ttAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCgg                                     >  1:2047840/1‑70 (MQ=255)
                               tAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGGTGACGATCACGATCCCAATAATCATCGCGGc                                    <  1:82978/70‑1 (MQ=255)
                               tAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCaataat                                             >  1:950562/1‑61 (MQ=255)
                               tAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCaataat                                             >  1:346023/1‑61 (MQ=255)
                               tAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCgg                                     >  1:2579881/1‑69 (MQ=255)
                               tAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCgg                                     >  1:2593402/1‑69 (MQ=255)
                               tAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGc                                    <  1:2718855/70‑1 (MQ=255)
                               tAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGc                                    <  1:579093/70‑1 (MQ=255)
                               tAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGc                                    >  1:1628420/1‑70 (MQ=255)
                                aaaGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCa                                   <  1:2786816/70‑1 (MQ=255)
                                 aaGACAACATCAGTTGTGACAACGTCACCTTGCGCGGTGACGATCACGATCCCAATAATCATCGCGGCaa                                  >  1:2608898/1‑70 (MQ=255)
                                    acaacaTCAGTTGTGACAACGTCACCTTGCGCGGTGACGATCACGATCCCAATAATCATCGCGGCAATTAc                              >  1:3877240/1‑71 (MQ=255)
                                    acaacaTCAGTTGTGACAACGTCACCTTGCGCGGTGACGATCACGATCCCAATAATCATCGCGGCAATTAc                              >  1:364686/1‑71 (MQ=255)
                                      aacaTCAGTTGTGACAACGTCACCTTGCGCGGTGACGATCACGATCCCAATAATCATCGCGGCaa                                  >  1:1091012/1‑65 (MQ=255)
                                      aacaTCAGTTGTGACAACGTCACCTTGCGCGGTGACGATCACGATCCCAATAATCATCGCGGCaa                                  >  1:2332547/1‑65 (MQ=255)
                                      aacaTCAGTTGTGACAACGTCACCTTGCGCGGTGACGATCACGATCCCAATAATCATCGCGGCaa                                  >  1:3996331/1‑65 (MQ=255)
                                      aacaTCAGTTGTGACAACGTCACCTTGCGCGGTGACGATCACGATCCCAATAATCATCGCGGCaa                                  >  1:1204942/1‑65 (MQ=255)
                                      aacaTCAGTTGTGACAACGTCACCTTGCGCGGTGACGATCACGATCCAATAATCATCGCGGCaa                                  >  1:2916143/1‑64 (MQ=255)
                                      aacaTCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATcc                                                    <  1:2687425/47‑1 (MQ=255)
                                      aacaTCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATcc                                                    <  1:3814340/47‑1 (MQ=255)
                                      aacaTCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCg                            >  1:4021915/1‑71 (MQ=255)
                                      aacaTCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCg                            <  1:425641/71‑1 (MQ=255)
                                      aacaTCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCg                            >  1:660937/1‑71 (MQ=255)
                                      aacaTCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCg                            >  1:1735137/1‑71 (MQ=255)
                                      aacaTCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCg                            >  1:2551298/1‑71 (MQ=255)
                                       acaTCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGc                           <  1:2079300/71‑1 (MQ=255)
                                         aTCAGTTGTGACAACGTCACCTTGCGCTATGACGATCACGATCCCAATAATCATCGCGGCaa                                  >  1:3472277/1‑62 (MQ=255)
                                         aTCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCaa                                  >  1:1949732/1‑62 (MQ=255)
                                         aTCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCaa                                  >  1:1294286/1‑62 (MQ=255)
                                          tCAGTTGTGACAACGTCACCTTGCGCGATGacgatca                                                          >  1:2064556/1‑37 (MQ=255)
                                          tCAGTTGTGACAACGTCACCTTGCGCGATGacgatca                                                          >  1:2358851/1‑37 (MQ=255)
                                          tCAGTTGTGACAACGTCACCTTGCGCGATGacgatca                                                          >  1:3030749/1‑37 (MQ=255)
                                          tCAGTTGTGACAACGTCACCTTGCGCGATGacgatc                                                           >  1:3765904/1‑36 (MQ=255)
                                          tCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCaata                                               <  1:2682096/48‑1 (MQ=255)
                                          tCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCaata                                               <  1:931475/48‑1 (MQ=255)
                                          tCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCaata                                               <  1:331247/48‑1 (MQ=255)
                                          tCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAatca                                           <  1:3534107/52‑1 (MQ=255)
                                          tCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAatca                                           <  1:4037962/52‑1 (MQ=255)
                                          tCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAatca                                           <  1:3529444/52‑1 (MQ=255)
                                          tCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAatca                                           <  1:2818216/52‑1 (MQ=255)
                                           cAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAatcat                                          >  1:514296/1‑52 (MQ=255)
                                           cAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAatcat                                          >  1:1314246/1‑52 (MQ=255)
                                            aGTTGTGACAACGTCACCTTGCGCGGTGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCAt                         >  1:512715/1‑68 (MQ=255)
                                            aGTTGTGACAACGTCACCTTGCGCGGTGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCAt                         >  1:3361664/1‑68 (MQ=255)
                                            aGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAatcatc                                         >  1:2651722/1‑52 (MQ=255)
                                            aGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAatcatc                                         >  1:1878728/1‑52 (MQ=255)
                                            aGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAatcatc                                         >  1:1846759/1‑52 (MQ=255)
                                              ttGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCaataat                                             <  1:1569425/46‑1 (MQ=255)
                                               tgtgACAACGTCACCTTGCCCGATGACGATCACGATCCCAATAatca                                           >  1:3853361/1‑47 (MQ=255)
                                                gtgACAACGTCACCTTGCGCGGTGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGcgcg                   >  1:2598085/1‑70 (MQ=255)
                                                gtgACAACGTCACCTTGCGCGGTGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCgcgg                  >  1:1476006/1‑71 (MQ=255)
                                                gtgACAACGTCACCTTGCGCGGTGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCgcgg                  >  1:3157997/1‑71 (MQ=255)
                                                gtgACAACGTCACCTTGCGCGGTGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCgcgg                  >  1:3924443/1‑71 (MQ=255)
                                                gtgACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCgcgg                  >  1:405176/1‑71 (MQ=255)
                                                gtgACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCgcgg                  >  1:2132210/1‑71 (MQ=255)
                                                           aCCTTGCGCGATGACGATCACGATCCCAATAATCATCGCgg                                     >  1:247002/1‑41 (MQ=255)
                                                           aCCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGc                           >  1:1168447/1‑51 (MQ=255)
                                                           aCCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCa                          >  1:2086870/1‑52 (MQ=255)
                                                           aCCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCa                          >  1:3507716/1‑52 (MQ=255)
                                                           aCCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCa                          >  1:3213885/1‑52 (MQ=255)
                                                           aCCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCa                          >  1:1774738/1‑52 (MQ=255)
                                                           aCCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCt       >  1:2388825/1‑71 (MQ=255)
                                                           aCCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCt       >  1:3570402/1‑71 (MQ=255)
                                                           aCCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCt       >  1:2216671/1‑71 (MQ=255)
                                                           aCCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCt       >  1:2665109/1‑71 (MQ=255)
                                                                 cgcgATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGtt  <  1:3626362/70‑1 (MQ=255)
                                                                     |                                                                 
TCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTT  >  W3110S.gb/1507181‑1507315

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