Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,584,214 T→G 8.5% E457A (GAG→GCG)  glgE ← alpha‑1,4‑glucan:maltose‑1‑phosphate maltosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,584,2140TG 8.5% 136.8 / 4.3 47E457A (GAG→GCG) glgEalpha‑1,4‑glucan:maltose‑1‑phosphate maltosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (21/22);  new base G (2/2);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCGGCGAGGAAGATGACGTCGGGGTGCT  >  NC_002947/4584067‑4584352
                                                                                                                                                   |                                                                                                                                          
cacaAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAATACGGGTTGATGTCCGGTGTATTGACCATGAAGTTGGGACGGTAGCACTGGCTCCACGGCGGCTGATAGAGCTGCTc                                                                                                                                           >  4:4518/1‑149 (MQ=255)
cacaAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTc                                                                                                                                           >  2:49097/1‑149 (MQ=255)
  caAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCTGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGa                                                                                                                                         <  7:103248/149‑1 (MQ=255)
     gCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCGCGaata                                                                                                                                      >  3:328753/1‑147 (MQ=255)
          aGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTcg                                                                                                                                 <  5:207341/149‑1 (MQ=255)
           gCCCATGGTCGCCAGCGCGGCACGGATCAGGTAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTcgc                                                                                                                                <  3:592/149‑1 (MQ=255)
            cccATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAACCGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTcgcg                                                                                                                               <  4:206561/149‑1 (MQ=255)
            cccATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTcgcg                                                                                                                               >  4:88563/1‑149 (MQ=255)
               aTGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAg                                                                                                                            >  7:169889/1‑149 (MQ=255)
               aTGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCGCGAATATTTCGAGCAg                                                                                                                            >  8:90903/1‑149 (MQ=255)
                  gTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTc                                                                                                                         >  1:269341/1‑149 (MQ=255)
                   tCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCGCGAATATTTCGAGCAGCTca                                                                                                                        <  4:328753/149‑2 (MQ=255)
                    cGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCGCGAATATTTCGAGCAGCTcaa                                                                                                                       <  7:90903/149‑3 (MQ=255)
                      ccAGCTCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTTTTCAGGAATAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGc                                                                                                                     <  8:93442/149‑1 (MQ=255)
                      ccAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGc                                                                                                                     <  2:269341/149‑1 (MQ=255)
                         gcgcGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTtg                                                                                                                  >  1:141203/1‑149 (MQ=255)
                          cgcgGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTtgt                                                                                                                 <  4:145962/149‑1 (MQ=255)
                           gcgGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGc                                                                                                                     >  2:167470/1‑144 (MQ=255)
                           gcgGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGc                                                                                                                     <  1:167470/144‑1 (MQ=255)
                              gCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGtt                                                                                                             <  5:96143/149‑1 (MQ=255)
                                    aTCAGGAAACCTGCGAGGACTGAGGAGTGCAGGAAAAACGGGTTGATGTCCGGAGTAATGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCGTGCTTGTGGTTgcgc                                                                                                         >  6:333129/1‑147 (MQ=255)
                                     tCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAg                                                                                                      <  3:88563/149‑1 (MQ=255)
                                         gAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTg                                                                                                   >  5:30463/1‑148 (MQ=255)
                                           acccTGCGCGGCCTGTGGTGTGCACGTAACCCGGGTTCATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACGGGCTCCACGGCGGCGGATTGAGCTGCTCGAAATATTCGCGCAGCTGCGGCTTGTGGTTGCGCCAGGTGaaa                                                                                                <  6:277335/147‑1 (MQ=255)
                                              cTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAg                                                                                             >  8:29921/1‑149 (MQ=255)
                                              cTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAg                                                                                             <  8:169889/149‑1 (MQ=255)
                                                 cgcgGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGta                                                                                          <  6:30463/149‑1 (MQ=255)
                                                  gcgGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGta                                                                                          >  8:243830/1‑148 (MQ=255)
                                                        ggAGGTGTGCAGGAAAACCGGGGTGATGTCGGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTTAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGctct                                                                                   <  3:266588/148‑1 (MQ=255)
                                                          aGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCtgt                                                                                 <  1:49097/149‑1 (MQ=255)
                                                               tgCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGaaa                                                                                                                                       <  8:268660/90‑1 (MQ=255)
                                                               tgCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGaaa                                                                                                                                       >  7:268660/1‑90 (MQ=255)
                                                                                  tGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGATGCGCCAGGTGAAATAGGTATAGCACTGTGCGTACCCGACCTTGCCCAGgcgc                                                         >  7:54332/1‑149 (MQ=255)
                                                                                      cGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGcc                                                      >  6:29644/1‑148 (MQ=255)
                                                                                      cGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCca                                                     >  2:127767/1‑149 (MQ=255)
                                                                                           tATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCatcatc                                                <  7:29921/149‑1 (MQ=255)
                                                                                                  gaagaaGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCt                                         >  8:309818/1‑149 (MQ=255)
                                                                                                  gaagaaGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGATCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCAGGCt                                         >  5:158505/1‑149 (MQ=255)
                                                                                                  gaagaaGTTGGGCCGGTAGAACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCt                                         >  2:16700/1‑149 (MQ=255)
                                                                                                   aagaagTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCtt                                        <  7:243830/149‑1 (MQ=255)
                                                                                                   aagaagTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCtt                                        >  5:239209/1‑149 (MQ=255)
                                                                                                              ccGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCtt                             <  5:29644/149‑1 (MQ=255)
                                                                                                                 gTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCgg                          <  1:127767/149‑1 (MQ=255)
                                                                                                                        tGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCGGCGAgg                     >  8:238738/1‑147 (MQ=255)
                                                                                                                          gCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCGGCGAgg                     >  4:211953/1‑145 (MQ=255)
                                                                                                                                  gcggcTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCGGCGAGGAAGATGACGTCg         <  7:238738/149‑1 (MQ=255)
                                                                                                                                         aTTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCGGCGAGGAAGATGACGTCGGGGTgct  <  6:239209/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
CACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATCATCGCCGGCTTGGTAAAGGCTTCGGCGAGGAAGATGACGTCGGGGTGCT  >  NC_002947/4584067‑4584352

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: ‑