Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,178,242 G→A 13.6% P48S (CCT→TCT)  gatA ← galactitol‑specific enzyme IIA component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,178,2420GA13.6% 114.9 / 9.6 44P48S (CCT→TCT) gatAgalactitol‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base G (19/19);  new base A (2/4);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

TCCTTAGCATGAATCGCCTCACAATGCGGTATTGCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTGGCAATTAATGCCTGTGGCCAGGTATCATGAACCACACCTTTGGCGAGCATC  >  W3110S.gb/2178175‑2178309
                                                                   |                                                                   
gacttagCATGAATCGCCTCACAATGCGGTATTGCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGGGaa                                                                  >  1:3442424‑M2/8‑71 (MQ=255)
gacttagCATGAATCGCCTCACAATGCGGTATTGCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGGGaa                                                                  <  1:3224510‑M2/64‑1 (MQ=255)
   ttagCATGAATCGCCTCACAATGCGGTATTGCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTc                                                               >  1:3296542‑M2/5‑71 (MQ=255)
     agCATGAATCGCCTCACAATGCGGTATTGCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTc                                                               <  1:1546101‑M2/67‑1 (MQ=255)
        aTGAATCGCCTCACAATGCGGTATTGCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGAGAATTCTGCtt                                                          <  1:3527343/71‑1 (MQ=255)
         tGAATCGCCTCACAATGCGGTATTGCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGc                                                            >  1:883851/1‑68 (MQ=255)
                 cTCACAATGCGGTATTGCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTGGCa                                                  <  1:2796721/70‑1 (MQ=255)
                  tCACAATGCGGTATTGCAATAGCGTGCTTCTCAAGCATTATCCCGGTAGGGAATTCTGCTtct                                                        >  1:2920102/1‑63 (MQ=255)
                  tCACAATGCGGTATTGCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTGGCaa                                                 >  1:2323923/1‑70 (MQ=255)
                  tCACAATGCGGTATTGCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGAGAATTCTGc                                                            <  1:898464/59‑1 (MQ=255)
                    acaATGCGGTATTGCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGAGAATTCTGCtt                                                          >  1:134606/1‑59 (MQ=255)
                       aTGCGGTATTGCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTGGCAATTAATg                                           >  1:421899/1‑71 (MQ=255)
                        tGCGGTATTGCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGGGaa                                                                  <  1:1397986/47‑1 (MQ=255)
                          cGGTATTGCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTGGCAATTAATGcc                                         <  1:2038443/70‑1 (MQ=255)
                             tATTGCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGAGAATTCTGCTtct                                                        <  1:3724299/52‑1 (MQ=255)
                               ttGCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCtt                                                          <  1:234397/48‑1 (MQ=255)
                               ttGCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCtt                                                          <  1:3712992/48‑1 (MQ=255)
                               ttGCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTGGCAATTAATGcc                                         <  1:2544071/65‑1 (MQ=255)
                               ttGCAATAGCGTGCTGCTAAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTGGCAATTAATGCCTGTGGc                                   >  1:3892524/1‑71 (MQ=255)
                                 gCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTGGCAATTAATGCCTGTGGc                                   >  1:2919021/1‑69 (MQ=255)
                                 gCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTGGCAATTAATGCCTGTGGc                                   >  1:3550280/1‑69 (MQ=255)
                                 gCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTGGCAATTAATGCCTGTGGc                                   >  1:1562789/1‑69 (MQ=255)
                                 gCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTGGCAATTAATGCCTGTGGc                                   >  1:317876/1‑69 (MQ=255)
                                 gCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTGGCAATTAATGACTGTGGc                                   >  1:3707822/1‑69 (MQ=255)
                                   aaTAGCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTGGCAATTAATGCCTGTGGCCAgg                               <  1:4370165/71‑1 (MQ=255)
                                    aTAGCGTGCTGCTCAAGCATTATCCCGGTAGGaaa                                                                  <  1:2063739/35‑1 (MQ=255)
                                    aTAGCGTGCTGCTCAAGCATTATCCCGGTAGGGaa                                                                  <  1:579840/35‑1 (MQ=255)
                                    aTAGCGTGCTGCTCAAGCATTATCCCGGTAGGGaa                                                                  <  1:4621229/35‑1 (MQ=255)
                                    aTAGCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGc                                                            <  1:1917283/41‑1 (MQ=255)
                                     tAGCGTGCTGCTCAAGCATTATCCCGGTAGGGAAtt                                                                >  1:3719508/1‑36 (MQ=255)
                                       gCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTg                                                     >  1:880760/1‑45 (MQ=255)
                                       gCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTg                                                     >  1:441830/1‑45 (MQ=255)
                                       gCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTGGCAATTAATGCCTGTGGCCAGGTAt                            <  1:231295/70‑1 (MQ=255)
                                       gCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTGGCAATTAATGCCTGTGGCCAGGTAt                            <  1:1123685/70‑1 (MQ=255)
                                       gCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTGGCAATTAATGCCTGTGGCCAGGTAt                            >  1:3950882/1‑70 (MQ=255)
                                       gCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTGGCAATTAATGCCTGTGGCCAGGTAt                            <  1:3139665/70‑1 (MQ=255)
                                       gCGTGCTGCTCAAGCATTATCCCGGTAGAGAATTCTGCTTCTCTGGCa                                                  >  1:4595038/1‑48 (MQ=255)
                                          tgctgcTCAAGCACTATCCCGGTAGGGAATTCTGCTTCTCTGGCAATTAATGCCTGTGGCCAGGTATCATg                        <  1:4164827/71‑1 (MQ=255)
                                           gctgctCAAGCATTATCCCGGTAGAGAATTCTGCTTCTCTGGCAATTAATGCCTGTGGCCAgg                               <  1:4379462/63‑1 (MQ=255)
                                            ctgctCAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTGGCAATTAATGCCTGTGGCCAGGTATCATGaa                      <  1:401057/71‑1 (MQ=255)
                                                     cATTATCCCGGTAGGGAATTCTGCTTCTCTGGCAATTaa                                             >  1:3045373/1‑39 (MQ=255)
                                                        tATCCCGGTAGGGAATTCTGCTTCTCTGGCAATTAATGCCTGTGGCCAGGTATc                           >  1:2510119/1‑54 (MQ=255)
                                                         aTCCCGGTAGGGAATTCTGCTTCTCTGGCAATTAATGCCTGTGGCCAGGTATCATGAACCACACCTTTGGc         >  1:1688626/1‑71 (MQ=255)
                                                                tAGGGAATTCTGCTTCTCTGGCAATTAATGCCTGTGGCCAGGTATCATGAACCACACCTTTGGCGAGCAtc  <  1:2789739/71‑1 (MQ=255)
                                                                   |                                                                   
TCCTTAGCATGAATCGCCTCACAATGCGGTATTGCAATAGCGTGCTGCTCAAGCATTATCCCGGTAGGGAATTCTGCTTCTCTGGCAATTAATGCCTGTGGCCAGGTATCATGAACCACACCTTTGGCGAGCATC  >  W3110S.gb/2178175‑2178309

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: