Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,585,553 T→C 12.0% S11G (AGC→GGC)  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,585,5530TC12.0% 127.8 / 7.8 50S11G (AGC→GGC) glgEalpha‑1,4‑glucan:maltose‑1‑phosphate maltosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (22/22);  new base C (4/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CCAGGCGGTCGTGGCCGTCGCTGTAGGCCTTGGTGCTGACGGCCACAGGCTGGCCGCTGACGGCCTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGTTCAGAGCCGGGCACAGGCGGGGGGTTCAAAATAATTGAGCGAACGGCAGC  >  NC_002947/4585406‑4585701
                                                                                                                                                   |                                                                                                                                                     
ccAGGCGGTCGTGGCCGTCGCTGTAGGCCTTGGTGCTGACGGCCACAGGCTGGCCGCTGACGGCCTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCt                                                                                                                                                       >  5:218873/1‑148 (MQ=255)
               cGTCGCTGTAGGCCTTGGTGCTGACGGCCACAGGCTGGCCGCTGACGGCCTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCggg                                                                                                                                       >  2:113549/1‑149 (MQ=255)
               cGTCGCTGTAGGCCTTGGTGCTGACGGCCACAGGCTGGCCGCTGACGGCCTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGACCAGGCCATTAGCCTGGcccag                                                                                                                                        >  2:311632/1‑146 (MQ=255)
               cGTCGCTGTAGGCCTTGGTGCTGACGGCCACAGGCTGGCCGCTGACGGCCTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGACCAGGCCATTAGCCTGGcccag                                                                                                                                        <  1:311632/148‑3 (MQ=255)
                            cTTGGTGCTGACGGCCACAGGCTGGCCGCTGACGGCCTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCa                                                                                                                          >  1:113705/1‑149 (MQ=255)
                                gTGCTGACGGCCACAGGCTGGCCGCTGACGGCCTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATgg                                                                                                                       <  8:41097/148‑1 (MQ=255)
                                          ccACAGGCTGGCCGCTGACGGCCTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTc                                                                                                             >  1:25132/1‑148 (MQ=255)
                                           cacaGGCTGGCCGCTGACGGCCTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCaa                                                                                                           <  4:59134/149‑1 (MQ=255)
                                            acaGGCTGGCCTCTGACGGCCTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCaaa                                                                                                          <  8:261745/149‑1 (MQ=255)
                                               ggCTGGCCGCTGACGGCCTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAggg                                                                                                       <  1:113549/149‑1 (MQ=255)
                                                      cGCTGACGGCCTTGGCGGTGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGtt                                                                                                 >  4:226997/1‑148 (MQ=255)
                                                       gCTGACGGCCTTGGCGGTGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTc                                                                                               <  3:226997/149‑1 (MQ=255)
                                                       gCTGACGGCCTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTc                                                                                               >  3:262611/1‑149 (MQ=255)
                                                        cTGACGGCCTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCg                                                                                              <  3:183014/149‑1 (MQ=255)
                                                        cTGACGGCCTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCg                                                                                              <  1:157856/149‑1 (MQ=255)
                                                         tGACGGCCTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGt                                                                                             >  7:121183/1‑149 (MQ=255)
                                                           gcGGCCTTGGCGGCGTATCTGCCGGCGTCGAGCACGGGTTGGGGGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCCGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGa                                                                                           <  6:200668/148‑1 (MQ=255)
                                                             ggCCTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGa                                                                                           >  5:200668/1‑147 (MQ=255)
                                                              gCCTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACAt                                                                                        >  2:181954/1‑149 (MQ=255)
                                                               gcTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTc                                                                                                             <  3:228414/126‑1 (MQ=255)
                                                               ccTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTc                                                                                                             >  4:228414/1‑127 (MQ=255)
                                                                cTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCg                                                                                      >  7:132983/1‑149 (MQ=255)
                                                                    gcggcgAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATTGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCt                                                                                  <  8:121183/149‑1 (MQ=255)
                                                                                ccGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAg                                                                      <  3:308830/148‑1 (MQ=255)
                                                                                    cGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCtgg                                                                 <  2:113705/149‑1 (MQ=255)
                                                                                       cGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCtggtgg                                                              <  8:132983/149‑1 (MQ=255)
                                                                                       cGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCtggtgg                                                              >  7:188884/1‑149 (MQ=255)
                                                                                         aGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCg                                                            <  4:262611/149‑1 (MQ=255)
                                                                                                      tgtCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGTTCa                                                 >  5:224559/1‑147 (MQ=255)
                                                                                                            cGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAggg                                                     <  4:3037/137‑1 (MQ=255)
                                                                                                            cGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAggg                                                     >  3:3037/1‑137 (MQ=255)
                                                                                                            cGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCCCCCTCAGGCTGGTGGCGGTAAGGGTTCAGAGcc                                            <  7:177102/146‑1 (MQ=255)
                                                                                                            agATCACTATGCGTGGCGCCAACAGCGACTGGGCCAGGCTAATGGACTGGTACGGGTGATCGTTCGCCATGTGCACGCTTTCAAAGGTCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGTTCAGAGcc                                            <  2:25132/145‑1 (MQ=255)
                                                                                                              aTCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGTTCAGAGCCGGGca                                       >  8:20607/1‑149 (MQ=255)
                                                                                                                  cGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGTTCAGAGCCGGGCACAgg                                   >  4:48707/1‑149 (MQ=255)
                                                                                                                   gATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGTTCAGAGCCGGGCACAGGc                                  >  5:27947/1‑149 (MQ=255)
                                                                                                                       cGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGTTCAGAGCCGGGCACAGGCgggg                              <  6:27947/149‑1 (MQ=255)
                                                                                                                        gTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGTTCAGAGCCGGGCACAGGCggggg                             >  3:275450/1‑149 (MQ=255)
                                                                                                                           gcgcCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGTTCAGAGCCGGGCACAGGCGGGGGGtt                          <  7:20607/149‑1 (MQ=255)
                                                                                                                            cgcCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCg                                                            <  8:158572/114‑1 (MQ=255)
                                                                                                                              ccAACAGCGCCTGGGCCAGGCTAATGGCCTCGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGTTCAGAGCCGGGCACAGGCGGGGGGTTCa                        >  2:221233/1‑148 (MQ=255)
                                                                                                                                      gcCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGTTCAGAGCCGGGCACAGGCGGGGGGTTCAAAATAAtt                >  5:306088/1‑148 (MQ=255)
                                                                                                                                        ccggaccAGGCCATTAGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTCGCTCCCTCAGGCTGGTGGCGGTAAGGGtt                                                   <  1:24768/107‑1 (MQ=255)
                                                                                                                                        ccggaccAGGCCATTAGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGtt                                                   >  2:24768/6‑111 (MQ=255)
                                                                                                                                          gggCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGCCCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGTTCAGAGCCGGGCACAGGCGGGGGGTTCAAAATAATTGAGc            >  4:143113/1‑148 (MQ=255)
                                                                                                                                          ggaccAGGCCATTAGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGTTCAGAGCCGGGCACAGGCGGGGGGTTCAAAATAATTGAGc            >  6:186665/4‑148 (MQ=255)
                                                                                                                                          ggaccAGGCCATTAGCCTGGCCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGTTCAGAGCCGGGCACAGGCGGGGGGTTCAAAATAATTGAGc            >  3:88987/4‑148 (MQ=255)
                                                                                                                                              cAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGTTCAGAGCCGGGCACAGGCGGGGGGTTCAAAATAATTGAGCGAACg       <  6:306088/149‑1 (MQ=255)
                                                                                                                                                ggCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGTTCAGAGCCGGGCACAGGCGGGGGGTTCAAAATAATTGAGCGAACgg      <  4:12826/148‑1 (MQ=255)
                                                                                                                                                ggCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGTTCAGAGCCGGGCACAGGCGGGGGGTTCAAAATAATTGAGCGAACgg      >  3:12826/1‑148 (MQ=255)
                                                                                                                                                 gCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGTTCAGAGCCGGGCACAGGCGGGGGGTTCAAAATAAATGAGCGAACgg      >  2:132476/1‑147 (MQ=255)
                                                                                                                                                   tAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGTTCAGAGCCGGGCACAGGCGGGGGGTTCAAAATAATTGAGCGAACGGCAGc  >  7:41095/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCAGGCGGTCGTGGCCGTCGCTGTAGGCCTTGGTGCTGACGGCCACAGGCTGGCCGCTGACGGCCTTGGCGGCGAAGCTGCCGGCGTCGAGCACGGGTTGGGTGTCTTCGATCACGATGCGTGGCGCCAACAGCGCCTGGGCCAGGCTAATGGCCTGGTCCGGGTGATCGTTCGCCATGGGCACGCTTTCAAAGGGCTCGTTTCGTGACATCGACCTTGCTCCCTCAGGCTGGTGGCGGTAAGGGTTCAGAGCCGGGCACAGGCGGGGGGTTCAAAATAATTGAGCGAACGGCAGC  >  NC_002947/4585406‑4585701

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