Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,521,572 G→A 37.5% L157L (CTG→CTA glpX → fructose 1,6‑bisphosphatase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,521,5720GA37.5% 51.5 / 58.2 56L157L (CTG→CTAglpXfructose 1,6‑bisphosphatase II
Reads supporting (aligned to +/- strand):  ref base G (16/19);  new base A (8/13);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 7.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

CCGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCGTACGC  >  W3110S.gb/3521513‑3521641
                                                           |                                                                     
ccGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACg                                                              <  1:3712222/69‑1 (MQ=255)
   cTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTAACGGTAACGATTCt                                                         <  1:2360532/71‑1 (MQ=255)
        ggATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTa                                                    <  1:1363008/71‑1 (MQ=255)
        ggATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTa                                                    <  1:2501082/71‑1 (MQ=255)
         gATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTaa                                                   <  1:3592669/71‑1 (MQ=255)
          aTAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTaaa                                                  >  1:3215265/1‑71 (MQ=255)
          aTAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTaaa                                                  >  1:2161743/1‑71 (MQ=255)
          aTAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTaaa                                                  >  1:514395/1‑71 (MQ=255)
           tAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAc                                                 >  1:1400806/1‑71 (MQ=255)
           tAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAAc                                                 <  1:4818247/71‑1 (MQ=255)
              ccTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAccac                                              <  1:3741495/71‑1 (MQ=255)
              ccTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAccac                                              <  1:4253426/71‑1 (MQ=255)
              ccTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAccac                                              <  1:1789239/71‑1 (MQ=255)
               cTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAAccacg                                             >  1:2558777/1‑71 (MQ=255)
                   gcAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAccacgccac                                         <  1:1559694/71‑1 (MQ=255)
                      aTGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAccac                                              <  1:2787724/63‑1 (MQ=255)
                      aTGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAccac                                              <  1:515216/63‑1 (MQ=255)
                       tGTAGCGGCCGCGCTCGGCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTaa                                                   <  1:157507/57‑1 (MQ=255)
                         tAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTg                                                                       >  1:3395190/1‑35 (MQ=255)
                         tAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGa                                                                      >  1:1058638/1‑36 (MQ=255)
                         tAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGcc                                   >  1:86956/1‑71 (MQ=255)
                         tAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGcc                                   >  1:4219587/1‑71 (MQ=255)
                             ggcggcgCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTa                                                    >  1:4256550/1‑50 (MQ=255)
                              gcggcgCTCGGCAAACCGTTGAGCGAACTGACGGTAAc                                                               <  1:295007/38‑1 (MQ=255)
                               cggcgCTCGGCAAACCGTTGAGCGAACTGACTGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATc                             <  1:3092228/71‑1 (MQ=255)
                                    cTCGGCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAAcc                                                <  1:4959025/47‑1 (MQ=255)
                                     tCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGa                        <  1:1998634/70‑1 (MQ=255)
                                      cGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGaaa                      >  1:1435282/1‑71 (MQ=255)
                                       ggCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAAt                     >  1:3722475/1‑71 (MQ=255)
                                       ggCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAAt                     >  1:3674026/1‑71 (MQ=255)
                                        gCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATg                    <  1:4118789/71‑1 (MQ=255)
                                        gCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATg                    <  1:964238/71‑1 (MQ=255)
                                         cAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATg                    >  1:3575324/1‑70 (MQ=255)
                                         cAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATg                    >  1:729648/1‑70 (MQ=255)
                                          aaaCCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATg                    >  1:2701391/1‑69 (MQ=255)
                                           aaCCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCt                          <  1:205764/62‑1 (MQ=255)
                                            aCCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATgcagc                <  1:2687301/71‑1 (MQ=255)
                                            aCCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATgcagc                <  1:1258290/71‑1 (MQ=255)
                                            aCCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATgcagc                <  1:3377123/71‑1 (MQ=255)
                                            aCCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATgcag                 >  1:4877483/1‑70 (MQ=255)
                                                 tGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGt                                 <  1:2850684/49‑1 (MQ=255)
                                                 tGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTc           >  1:593538/1‑71 (MQ=255)
                                                  gAGCGAACTAACGGTAACGATTCTGGCTAAAccacgccac                                         <  1:4159574/40‑1 (MQ=255)
                                                  gAGCGAACTAACGGTAACGATTCTGGCTAAAccacgccac                                         <  1:4539824/40‑1 (MQ=255)
                                                  gAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGa                                       >  1:2070174/1‑42 (MQ=255)
                                                  gAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGa                                       >  1:2064291/1‑42 (MQ=255)
                                                     cGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCg       >  1:189774/1‑71 (MQ=255)
                                                     cGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCg       <  1:2860542/71‑1 (MQ=255)
                                                     cGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCg       <  1:3531641/71‑1 (MQ=255)
                                                      gAACTGACGGTAACGATTCTGGCTAAACCATGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGc        >  1:2210830/1‑69 (MQ=255)
                                                      gAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGc        >  1:1193420/1‑69 (MQ=255)
                                                      gAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGc        >  1:4829160/1‑69 (MQ=255)
                                                      gAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCGt      <  1:2283205/71‑1 (MQ=255)
                                                         cTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGc                           <  1:1960028/47‑1 (MQ=255)
                                                         cTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGc                           <  1:3323221/47‑1 (MQ=255)
                                                          tGACGGTAACGATTCTTGCTAAACCACGCCACGATGACGTTATCGc                           <  1:132216/46‑1 (MQ=39)
                                                          tGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCGTAcgc  >  1:2141466/1‑71 (MQ=255)
                                                           |                                                                     
CCGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCGTACGC  >  W3110S.gb/3521513‑3521641

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: