Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 4,557,768 G→A 8.1% A472T (GCC→ACC)  glgA → Glycogen synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,557,7680GA 8.1% 153.4 / 2.7 51A472T (GCC→ACC) glgAGlycogen synthase
Reads supporting (aligned to +/- strand):  ref base G (18/28);  new base A (2/2);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01

TACGCGCAACGCTTCGGCTCGCTGCCGGTAGCACGTAATACCGGTGGCCTGGCCGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGAACCCTATGCCCGCCTGTACGAGGACCTGGTCAAGCAGACACAACTG  >  ACB143/4557636‑4557902
                                                                                                                                    |                                                                                                                                      
tACGCGCAACGCTTCGGCTCGCTGCCGGTAGCACGTAATACCGGTGGCCTGGCCGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCt                                                                                                                                     >  4:211657/1‑136 (MQ=255)
tACGCGCAACGCTTCGGCTCGCTGCCGGTAGCACGTAATACCGGTGGCCTGGCCGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCt                                                                                                                                     >  4:169632/1‑136 (MQ=255)
tACGCGCAACGCTTCGGCTCGCTGCCGGTAGCACGTAATACCGGTGGCCTGGCCGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCt                                                                                                                                     <  3:99237/136‑1 (MQ=255)
         cggttCGGCTCGCTGCCGGTAGCACGTAATACCGGTGGCCTGGCCGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGttt                                                                                                                            <  2:210908/133‑1 (MQ=255)
          gCTTCGGCTCGCTGCCGGTAGCACGTAATACCGGTGGCCTGGCCGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTa                                                                                                                           <  4:240406/136‑1 (MQ=255)
          gCTTCGGCTCGCTGCCGGTAGCACGTAATACCGGTGGCCTGGCCGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCGTCTATGTTTa                                                                                                                           >  1:213354/1‑136 (MQ=255)
          gCTTCGGCTCGCTGCCGGTAGCACGTAATACCGGTGGCCTGGCCGACACCATAGACTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGACGAGAGCTACCGCGACGCACTCAGACGCCCCATCTATGttt                                                                                                                            >  3:144948/1‑135 (MQ=255)
               ggCTCGCTGCCGGTAGCACGTAATACCGGTGGCCTGGCCGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCa                                                                                                                      <  5:162419/136‑1 (MQ=255)
                 cTCGCTGCCGGTAGCACGTAATACCGGTGGCCTGGCCGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGc                                                                                                                       >  7:226929/1‑133 (MQ=255)
                 cTCGCTGCCGGTAGCACGTAATACCGGTGGCCTGGCCGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGc                                                                                                                       >  7:2122/1‑133 (MQ=255)
                      tGCCGGTAGCACGTAATACCGGTGGCCTGGCCGCCACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACACAGCCGCGCCTTCTATGTTTACGGCAAGAAGGa                                                                                                               >  4:241778/1‑136 (MQ=255)
                         cGGTAGCACGTAATACCGGTGGCCTGGCCGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACct                                                                                                            <  1:35206/136‑1 (MQ=255)
                         cGGTAGCACGTAATACCGGTGGCCTGGCCGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCACCTTCTATGTTTACGGCAAGAAGGACct                                                                                                            >  5:134021/1‑136 (MQ=255)
                         cGGTAGCACGTAATACCGGTGGCCTGGCCGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCACCTTCTATGTTTACGGCAAGAAGGACct                                                                                                            >  5:129482/1‑136 (MQ=255)
                           gTAGCACGTAATACCGGTGGCCTGGCCGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACctgc                                                                                                          <  1:68828/136‑1 (MQ=255)
                              gCACGTAATACCGGTGGCCTGGCCGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACctgctg                                                                                                        >  6:250012/1‑135 (MQ=255)
                              gCACGTAATACCGGTGGCCTGGCCGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCATCTATGTTTACGGCAAGAAGGACctgctg                                                                                                        >  6:146512/1‑135 (MQ=255)
                                       aCCGGTGGCCTGGCCGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATgc                                                                                              >  2:136011/1‑136 (MQ=255)
                                         cGGTGGCCTGGCCGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCATCGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATgcgc                                                                                            <  3:84350/136‑1 (MQ=255)
                                              gcctggccGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCt                                                                                       >  2:174528/1‑136 (MQ=255)
                                              gcctggccGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCt                                                                                       >  7:223105/1‑136 (MQ=255)
                                              gcctggccGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCATCGATGTTTACGCCAAGAACGACCTGCAGAACGCCATGCGCTGCCt                                                                                       >  5:85623/1‑136 (MQ=255)
                                               cctggccGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTg                                                                                      >  6:136751/1‑136 (MQ=255)
                                                        cacCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGAcc                                                                             <  4:19021/136‑1 (MQ=255)
                                                                aGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGtt                                                                     >  7:14563/1‑136 (MQ=255)
                                                                 gTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTc                                                                    <  3:73049/136‑1 (MQ=255)
                                                                        gTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGt                                                             >  3:221209/1‑136 (MQ=255)
                                                                           aagggCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGcc                                                          <  5:172574/134‑1 (MQ=255)
                                                                           aCGGGCTTCCTGTTCAACGACTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGcc                                                          <  1:283917/136‑1 (MQ=255)
                                                                                   ccTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCgg                                                    <  5:39041/134‑1 (MQ=255)
                                                                                   ccTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCgg                                                    <  7:261908/134‑1 (MQ=255)
                                                                                   ccTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCgg                                                    <  4:240655/134‑1 (MQ=255)
                                                                                          aaCGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGAACCCt                                           <  5:47540/136‑1 (MQ=255)
                                                                                          aaCGAGGCCACCGTCGAGAGCTCCCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGAACCCt                                           <  7:176629/136‑1 (MQ=255)
                                                                                           aCGCGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGAACCCTa                                          <  8:140188/136‑1 (MQ=255)
                                                                                           aCGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCTATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGAACCCTa                                          <  4:22258/136‑1 (MQ=255)
                                                                                            cGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTTGTGCCAGGCGGTTGAACCCTAt                                         <  5:112781/136‑1 (MQ=255)
                                                                                            cGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGAACCCTAt                                         <  5:1211/136‑1 (MQ=255)
                                                                                                tCCAGCGTCGAGATCTTCCGCGAGGGACTCAGCCGCGGCTTCTATGTTTTCTGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGTCGGTTGAACCCTATGccc                                     <  4:144948/136‑1 (MQ=255)
                                                                                                tCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGAACCCTATGccc                                     <  8:77173/136‑1 (MQ=255)
                                                                                                  cACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGAACCCTATGCCCg                                    <  1:212008/135‑1 (MQ=255)
                                                                                                    ccGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGAACCCTATGCCCGCCt                                 >  8:163573/1‑136 (MQ=255)
                                                                                                       tCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGCACCCTATGCCCGCCTGTa                              >  1:282581/1‑136 (MQ=255)
                                                                                                                 ccGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGAACCCTATGCCCGCCTGTACGAGGACCt                     <  8:158192/135‑1 (MQ=255)
                                                                                                                   gcgAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGAACCCTATGCCCGCCTGTACGAGGACCTGGt                  <  3:24891/136‑1 (MQ=255)
                                                                                                                    cgAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGAACCCTATGCCCGCCTGTACGAGGACCTGGt                  >  3:258853/1‑135 (MQ=255)
                                                                                                                        gCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGAACCCTATGCCCGCCTGTACGAGGACCTGGTCAAGc             <  3:241778/136‑1 (MQ=255)
                                                                                                                              aGCCGCACCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGAACCCTATGCCCGCCTGTACGAGGACCTGGTCAAGCAGacac       <  6:134021/136‑1 (MQ=255)
                                                                                                                              aGCCGCACCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGAACCCTATGCCCGCCTGTACGAGGACCTGGTCAAGCAGacac       <  6:129482/136‑1 (MQ=255)
                                                                                                                                ccGCGCCTTCTATGTTTACGTCAAGAAGGACCTGCTGAACGCCAAGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGAACCCTATGCCCGCCTGTACGAGGACCTGGTCAAGCAGACACaa     <  5:41048/136‑1 (MQ=255)
                                                                                                                                ccGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGAACCCTATGCCCGCCTGTACGAGGACCTGGTCAAGCAGACACaa     <  1:2524/136‑1 (MQ=255)
                                                                                                                                   cgcCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGAACCCTATGCCCGCCTGTACGAGGACCTGGTCAAGCAGACACAACTg  <  3:74303/136‑1 (MQ=255)
                                                                                                                                   cgcCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGAACCCTATGCCCGCCTGTACGAGGACCTGGTCAAGCAGACACAACTg  <  3:133613/136‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                      
TACGCGCAACGCTTCGGCTCGCTGCCGGTAGCACGTAATACCGGTGGCCTGGCCGACACCATCGAGTGCGGCGTGACGGGCTTCCTGTTCAACGAGTCCACCGTCGAGAGCTACCGCGAGGCACTCAGCCGCGCCTTCTATGTTTACGGCAAGAAGGACCTGCTGAACGCCATGCGCTGCCTGTCGATGACCCAGCCGTTCAACTGGTGCCAGGCGGTTGAACCCTATGCCCGCCTGTACGAGGACCTGGTCAAGCAGACACAACTG  >  ACB143/4557636‑4557902

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