Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,867,919 C→T 13.3% G73D (GGC→GAC)  mtlA ← fused mannitol‑specific PTS enzyme IIABC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,867,9190CT13.3% 120.8 / 9.3 45G73D (GGC→GAC) mtlAfused mannitol‑specific PTS enzyme IIABC components
Reads supporting (aligned to +/- strand):  ref base C (19/20);  new base T (3/3);  total (22/23)
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.96e-01

GGCATGTCTGCGCCGACGATAACGCCCATGGTGGTGATGGCACCGACTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGGAGATAAGTGATCATCGGCCCGA  >  W3110S.gb/3867853‑3867982
                                                                  |                                                               
ggCATGTCTGCGCCGACGATAACGCCCATGGTGGTGATGGCACCGACTACGCCGCCACGTTCGCCGCCTa                                                              >  1:899539/1‑70 (MQ=255)
ggCATGTCTGCGCCGACGATAACGCCCATGGTGGTGATGGCACCGACTACGCCGCCACGTTCGCCGCCTAc                                                             <  1:4415369/71‑1 (MQ=255)
ggCATGTCTGCGCCGACGATAACGCCCATGGTGGTGATGGCACCGACTACGCCGCCACGTTCGCCGCCTAc                                                             >  1:3619392/1‑71 (MQ=255)
    tGTCTGCGCCGACGATAACGCCCATGGTGGTGATGGCACCGACTACGCCGCCACGTTCGCCGCCTACCAGc                                                         <  1:1310315/71‑1 (MQ=255)
      tCTGCGCCGACGATAACGCCCATGGTGGTGATGGCACCGACTACGCCGCCACGTTCGCCGTCTACCAGCt                                                        >  1:2935710/1‑70 (MQ=255)
      tCTGCGCCGACGATAACGCCCATGGTGGTGATGGCACCGACTACGCCGCCACGTTCGCCGCCTACCAGCt                                                        >  1:3418846/1‑70 (MQ=255)
          cgcCGACGATAACGCCCATGGTGGTGATGGCACCGACTACGCCGCCACGTTCGCCGCCTACCAGCTTacca                                                   >  1:3261211/1‑71 (MQ=255)
                 gaTAACGCCCATGGTGGTGATGGCACCGACTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGta                                             >  1:488213/1‑70 (MQ=255)
                 gaTAACGCCCATGGTGGTGATGGCACCGACTACGCCGCCACGTTCGCCGCCTAACAGCTTACCACCGGta                                             >  1:1021574/1‑70 (MQ=255)
                  aTAACGCCCATGGTGGTGATGGCACCGACTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGtata                                           <  1:2614376/71‑1 (MQ=255)
                     aCGCCCATGGTGGTGATGGCACCGACTACGCCGCCACGTTCGCCTCCTACCAGCTTACCACCGGtat                                            <  1:2890891/67‑1 (MQ=255)
                     aCGCCCATGGTGGTGATGGCACCGACTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGtat                                            <  1:4153177/67‑1 (MQ=255)
                     aCGCCCATGGTGGTGATGGCACCGACTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAAcc                                        <  1:2217725/71‑1 (MQ=255)
                     aCGCCCATGGTGGTGATGGCACCGACTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAAcc                                        <  1:4511944/71‑1 (MQ=255)
                           aTGGTGGTGATGGCACCGACTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATcag                                  >  1:4301094/1‑71 (MQ=255)
                            tggtggTGATGGCACCGACTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATcagc                                 <  1:2381653/71‑1 (MQ=255)
                                 gtgATGGCACCGACTACGCCGCCACGTTCGCCGCCTACCa                                                           <  1:3981104/40‑1 (MQ=255)
                                   gATGGCACCGACTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCgg                                               >  1:2687956/1‑50 (MQ=255)
                                      ggCACCGACTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAg                         >  1:4344385/1‑69 (MQ=255)
                                      ggCACCGACTACGACGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAg                         >  1:2925455/1‑69 (MQ=255)
                                       gCACCGACTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGta                                             >  1:2039956/1‑48 (MQ=255)
                                       gCACCGACTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAg                         >  1:3138195/1‑68 (MQ=255)
                                       gCACCGACTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAg                         >  1:3668968/1‑68 (MQ=255)
                                        cACCGACTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGgaga                     <  1:4581809/71‑1 (MQ=255)
                                         aCCGACTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGgaga                     >  1:625069/1‑70 (MQ=255)
                                         aCCGACTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGGAGAt                    <  1:3104730/71‑1 (MQ=255)
                                          ccGACTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGC‑‑CAGGAGATAAg                 >  1:277217/1‑71 (MQ=255)
                                            gACTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGGAGATatg                 <  1:3946511/71‑3 (MQ=255)
                                            gACTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGGAGATAAg                 >  1:4318763/1‑71 (MQ=255)
                                             aCTACGCCGCCACGTTCGCCGTCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGGAGATAAGt                <  1:2964918/71‑1 (MQ=255)
                                             aCTACGCCGCCACGTTCGCCGTCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGGAGATAAGt                <  1:2970178/71‑1 (MQ=255)
                                             aCTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGGAGATAAGt                >  1:3157760/1‑71 (MQ=255)
                                              cTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGGAGATAAGt                >  1:2759005/1‑70 (MQ=255)
                                              cTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGGAGATAAGTg               <  1:1150121/71‑1 (MQ=255)
                                              cTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGGAGATAAGTg               <  1:3182901/71‑1 (MQ=255)
                                              cTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGGAGATAAGTg               <  1:164526/71‑1 (MQ=255)
                                                aCGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGtata                                           <  1:4256462/41‑1 (MQ=255)
                                                 cgccgcCACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGGAGATAAGTGa              <  1:555040/69‑1 (MQ=255)
                                                   ccgccACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGGAGATAAGTGatcat          <  1:4114823/71‑1 (MQ=255)
                                                      ccACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGGAGAt                    <  1:4401253/58‑1 (MQ=255)
                                                       cACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGGAGATAAGTg               <  1:954626/62‑1 (MQ=255)
                                                        aCGTTCGCCGTCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGGAGATAAGTGATCATCGGcc     >  1:3686063/1‑71 (MQ=255)
                                                        aCGTTCGCCGTCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGGAGATAAGTGATCATCGGcc     >  1:2686166/1‑71 (MQ=255)
                                                         cGTTCGCCGTCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGGAGATAAGTGa              <  1:3829543/61‑1 (MQ=255)
                                                           ttCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGGAGATAAGTGATCATCGGCCCGa  >  1:1501989/1‑71 (MQ=255)
                                                                  |                                                               
GGCATGTCTGCGCCGACGATAACGCCCATGGTGGTGATGGCACCGACTACGCCGCCACGTTCGCCGCCTACCAGCTTACCACCGGTATAACCGATCAGCAGCGGCAGGAGATAAGTGATCATCGGCCCGA  >  W3110S.gb/3867853‑3867982

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: