Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 4,071,423 +A coding (335/1296 nt) glgC → glucose‑1‑phosphate adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,071,4191.A96.2% 238.8 / 3.8 80coding (331/1296 nt)glgCglucose‑1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base A (51/26);  total (51/29)
Fisher's exact test for biased strand distribution p-value = 4.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ATGAAGAAATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTCGACATTATCCGCCGTTATAAAGC  >  W3110S.gb/4071363‑4071486
                                                         |                                                                   
aTGAAGAAATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGaaaa                                                                 <  1:2240045/61‑1 (MQ=255)
aTGAAGAAATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGaaaa                                                                 <  1:2271367/61‑1 (MQ=255)
aTGAAGAAATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGaaaa                                                                 <  1:2483123/61‑1 (MQ=255)
aTGAAGAAATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGaaaa                                                                 <  1:3356088/61‑1 (MQ=255)
  gaagaaATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGaaaaac                                                               >  1:1029366/1‑59 (MQ=255)
  gaagaaATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGaaaaac                                                               >  1:1791354/1‑59 (MQ=255)
  gaagaaATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGaaaaac                                                               >  1:2943866/1‑59 (MQ=255)
     gaaATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCa                                                   >  1:3650380/1‑71 (MQ=255)
     gaaATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCa                                                   >  1:332496/1‑71 (MQ=255)
     gaaATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCa                                                   >  1:3173948/1‑71 (MQ=255)
     gaaATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCa                                                   >  1:3087400/1‑71 (MQ=255)
     gaaATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCa                                                   >  1:3731116/1‑71 (MQ=255)
     gaaATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCa                                                   >  1:2535132/1‑71 (MQ=255)
     gaaATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCa                                                   >  1:519533/1‑71 (MQ=255)
     gaaATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCa                                                   >  1:117614/1‑71 (MQ=255)
     gaaATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCa                                                   >  1:1903038/1‑71 (MQ=255)
        aTGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCg                                                >  1:3538104/1‑71 (MQ=255)
        aTGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCg                                                >  1:36686/1‑71 (MQ=255)
        aTGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCg                                                >  1:2527648/1‑71 (MQ=255)
        aTGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCg                                                >  1:988303/1‑71 (MQ=255)
        aTGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCg                                                >  1:2112400/1‑71 (MQ=255)
        aTGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCg                                                >  1:2082153/1‑71 (MQ=255)
        aTGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCg                                                >  1:1093478/1‑71 (MQ=255)
        aTGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCg                                                >  1:1127960/1‑71 (MQ=255)
        aTGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCg                                                >  1:1667825/1‑71 (MQ=255)
        aTGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCg                                                >  1:1654283/1‑71 (MQ=255)
         tGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCg                                                >  1:1470808/1‑70 (MQ=255)
         tGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCg                                                >  1:4053087/1‑70 (MQ=255)
         tGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCg                                                >  1:3895563/1‑70 (MQ=255)
         tGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCg                                                >  1:1457918/1‑70 (MQ=255)
         tGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCg                                                >  1:316080/1‑70 (MQ=255)
         tGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCg                                                >  1:294274/1‑70 (MQ=255)
         tGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCg                                                >  1:2261225/1‑70 (MQ=255)
         tGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCg                                                >  1:1674234/1‑70 (MQ=255)
             cGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGAt                                           >  1:1810627/1‑71 (MQ=255)
               aGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGc                                         <  1:3079418/71‑1 (MQ=255)
               aGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGc                                         <  1:4056162/71‑1 (MQ=255)
               aGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGc                                         <  1:1444146/71‑1 (MQ=255)
                   tGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGt                                      <  1:59809/70‑1 (MQ=255)
                   tGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGt                                      <  1:2209186/70‑1 (MQ=255)
                     tCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCa                                    <  1:2912824/70‑1 (MQ=255)
                     tCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCa                                    <  1:3994764/70‑1 (MQ=255)
                     tCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCa                                    <  1:3619012/70‑1 (MQ=255)
                     tCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCAc                                   <  1:160641/71‑1 (MQ=255)
                     tCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCAc                                   <  1:1893732/71‑1 (MQ=255)
                     tCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCAc                                   <  1:4106502/71‑1 (MQ=255)
                     tCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCAc                                   <  1:425482/71‑1 (MQ=255)
                      cGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCAcc                                  >  1:2824899/1‑71 (MQ=255)
                      cGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCAcc                                  >  1:3138044/1‑71 (MQ=255)
                      cGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCAcc                                  >  1:1176554/1‑71 (MQ=255)
                         tCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCa                                >  1:131611/1‑70 (MQ=255)
                             ctgcCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGGACCGCAGATGCGGTCAcc                                  >  1:1774655/1‑64 (MQ=255)
                             ctgcCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCAcc                                  >  1:1368491/1‑64 (MQ=255)
                             ctgcCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCAcc                                  >  1:4077866/1‑64 (MQ=255)
                             ctgcCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCAcc                                  >  1:3281753/1‑64 (MQ=255)
                             ctgcCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCAcc                                  >  1:2389908/1‑64 (MQ=255)
                             ctgcCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCAcc                                  >  1:1567870/1‑64 (MQ=255)
                             ctgcCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCAcc                                  >  1:156775/1‑64 (MQ=255)
                             ctgcCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCaaaa                             <  1:625217/69‑1 (MQ=255)
                             ctgcCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCaaaa                             <  1:3602405/69‑1 (MQ=255)
                             ctgcCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAAc                            <  1:2958158/70‑1 (MQ=255)
                             ctgcCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAAc                            <  1:573778/70‑1 (MQ=255)
                             ctgcCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAAc                            <  1:1222953/70‑1 (MQ=255)
                             ctgcCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAAc                            <  1:2872119/70‑1 (MQ=255)
                             ctgcCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAAc                            <  1:3041993/70‑1 (MQ=255)
                             ctgcCAGCACAGCAGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAAc                            <  1:1506966/70‑1 (MQ=255)
                             ctgcCAGCACAGCAGAGAATGAAAGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAAc                            <  1:2607929/68‑1 (MQ=255)
                                  agcacagcaGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTCGa                       >  1:2940584/1‑70 (MQ=255)
                                    cacagcaGAGAATGAAAGGGGAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTCGACAtt                   <  1:2351778/71‑1 (MQ=255)
                                    cacagcaGAGAATGAAAGGGGAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTCGACAtt                   <  1:3151213/71‑1 (MQ=255)
                                     acagcaGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTCGACAt                    >  1:274580/1‑70 (MQ=255)
                                     acagcaGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTCGACAt                    >  1:1941999/1‑70 (MQ=255)
                                     acagcaGAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTCGACAt                    >  1:255746/1‑70 (MQ=255)
                                      cagcagAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTCGACATTa                  >  1:2899812/1‑71 (MQ=255)
                                      cagcagAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTCGACATTa                  >  1:2762389/1‑71 (MQ=255)
                                      cagaagAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTCGACAt                    >  1:3706764/1‑69 (MQ=255)
                                       agcagAGAATGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCa                                    <  1:3318096/52‑1 (MQ=38)
                                              actGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTc                         <  1:44312/54‑1 (MQ=38)
                                              aaTGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTc                         <  1:157455/56‑1 (MQ=39)
                                              aaTGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTc                         <  1:1838947/56‑1 (MQ=39)
                                              aaTGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTCGACATTATCCGCCGt          >  1:1296444/1‑71 (MQ=255)
                                              aaTGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTCGACATTATCCGCCGt          >  1:3283671/1‑71 (MQ=255)
                                              aaTGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTCGACATTATCCGCCGt          >  1:1366670/1‑71 (MQ=255)
                                              aaTGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTCGACATTATCCGCCGt          >  1:705285/1‑71 (MQ=255)
                                              aaTGAAAGGGGAAAAACTGGTATCGCGGAACCGCAGATGCGGTCACACAAAACCTCGACATTATCCGCCGt          >  1:2391568/1‑71 (MQ=255)
                                               aTGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTc                         <  1:3880645/55‑1 (MQ=39)
                                                tGAAAGGGGAAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAAc                            <  1:3830877/51‑1 (MQ=38)
                                                   aggggaAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTCGACATTATCCGCCGTtata      >  1:78307/3‑69 (MQ=255)
                                                     gggaAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTCGACATTGTCCGCCGTTATAAAGc  <  1:1965690/71‑1 (MQ=255)
                                                     gggaAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTCGACATTATCCGCCGTTATAAAGc  <  1:2600499/71‑1 (MQ=255)
                                                         |                                                                   
ATGAAGAAATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAGAATGAAAGGGGAAAACTGGTATCGCGGCACCGCAGATGCGGTCACCCAAAACCTCGACATTATCCGCCGTTATAAAGC  >  W3110S.gb/4071363‑4071486

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: