Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,036,522 A→G 7.1% E351G (GAG→GGG)  algC → phosphomannomutase/phosphoglucomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,036,5220AG 7.1% 162.6 / 2.0 56E351G (GAG→GGG) algCphosphomannomutase/phosphoglucomutase
Reads supporting (aligned to +/- strand):  ref base A (28/24);  new base G (2/2);  total (30/26)
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.76e-01

GATGTGGAAAACCGGCCATTCGCTGATCAAGAAGAAGATGAAGCAGACCGGTTCGCTGCTGGCAGGCGAGATGAGCGGTCACATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAACGCGACGCCGACTGGGGC  >  NC_002947/6036376‑6036667
                                                                                                                                                   |                                                                                                                                                   
gATGTGGAAAACCGGCCATTCGCTGATCAAGAAGAAGATGAAGCAGACCGGTTCGCTGCTGGCAGGCGAGATGAGCGGTCACATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGgag                                                                                                                                                    >  3:103069/1‑148 (MQ=255)
    tgGAAAACCGGCCATTCGCTGATCAAGAAGAAGATGAAGCAGACCGGTTCGCTGCTGGCAGGCGAGATGAGCGGTCACATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATcc                                                                                                                                                <  7:277707/148‑1 (MQ=255)
      gAAAACCGGCCATTCGCTTATCAAGAAGAAGATGAATCAGACCGGTTCGGTGCTGGCAGTCGAGATGAGCGGTCACATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCa                                                                                                                                             <  5:56844/149‑1 (MQ=255)
           ccGGCCATTCGCTGATCAAGAAGAAGATGAAGCAGACCGGTTCGCTGCTGGCAGGCGAGATGAGCGGTCACATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAg                                                                                                                                        <  2:41715/149‑1 (MQ=255)
                 aTTCGCTGATCAAGAAGAAGATGAAGCAGACCGGTTCGCTGCTGGCAGGCGAGATGAGCGGTCACATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAg                                                                                                                                  >  1:129250/1‑149 (MQ=255)
                         aTCAAGAAGAAGATGAAGCAGACCGGTTCGCTGCTGGCAGGCGAGCTGAGCGGTCACATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAgcgc                                                                                                                         >  5:89550/1‑149 (MQ=255)
                          tCAAGAAGAAGATGAAGCAGACCGGTTCGCTGCTGGCAGGCGAGATGAGCGGTCACATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAgcgc                                                                                                                         <  4:225592/148‑1 (MQ=255)
                             agaagaagaTGAAGCAGACCGGTTCGCTGCTGGCAGGCGAGATGAGCGGTCACATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGa                                                                                                                     >  6:309580/1‑148 (MQ=255)
                                 gaagaTGAAGCAGACCGGTTCGCTGCTGGCAGGCGAGATGAGCGGTCACATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAAcc                                                                                                                >  2:254133/1‑149 (MQ=255)
                                               ccGGTTCGCTGCTGGCAGGCGAGATGAGCGGTCACATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCtt                                                                                                   <  3:110350/148‑1 (MQ=255)
                                                      gctgctGGCAGGCGAGATGAGCGGTCACTATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAAc                                                                                            >  1:18680/1‑149 (MQ=255)
                                                      gctgctGGCAGGCGAGATGAGCGGTCACATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAATACCGAGCAGAGCGCCGAAAACATGTTTGCCGCCTTCCCGAAc                                                                                            >  7:199160/1‑148 (MQ=255)
                                                             gCAGGCGAGATGAGCGGTCACATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTc                                                                                    >  7:198809/1‑149 (MQ=255)
                                                                 gCGAGATGAGCGGTCACATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACg                                                                                >  6:244996/1‑149 (MQ=255)
                                                                       tGAGCGGTCACATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTCCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGa                                                                           >  7:246736/1‑148 (MQ=255)
                                                                          gCGGTCACATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGAATGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCATTCCCGAACGACATTTCCACGCCGGAAAt                                                                        >  8:39836/1‑148 (MQ=255)
                                                                           aggTCACATCTTCATCAATGAACGCTGGTACGGTTTTTACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCTAAAACCTGTTTGCCGCATTCCCTAACGACATTTCCACGCCGGAAATc                                                                       <  4:103069/147‑1 (MQ=255)
                                                                             gTCACATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCaa                                                                     >  8:170352/1‑148 (MQ=255)
                                                                              tCACATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAata                                                                   <  7:170352/149‑1 (MQ=255)
                                                                                acaTCTTAATCAAGGTACGCTGGTACGGGTTTGACGACGGCATCTATATCGCCGCGCGCCTGCTGGAGATCCTCATCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAatat                                                                  <  8:199160/148‑1 (MQ=255)
                                                                                      ttCATCAACGATCGCTGGTTCGGTTTTGTCGACGGGATCTATATCCCCGCGCGGCTGCTGGAGATCCTCACCAAGACCGATCAGAGCGGCTCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGAt                                                              <  1:254133/149‑1 (MQ=255)
                                                                                        catcaAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGAtgtg                                                           <  5:244996/148‑1 (MQ=255)
                                                                                           caAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAATTCAATATTGATGTGAcc                                                        >  5:201064/1‑148 (MQ=255)
                                                                                           caAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGAcc                                                        >  1:314034/1‑148 (MQ=255)
                                                                                           caAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGAGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACAGAGCAGAGCGCCGAAAACCTGTTTGCAGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATTTGa                                                          >  1:318848/1‑146 (MQ=255)
                                                                                            aaGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCa                                                                      <  7:239525/133‑1 (MQ=255)
                                                                                            aaGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCa                                                                      >  8:239525/1‑133 (MQ=255)
                                                                                            aaGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTGTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGGTACcga                                                      >  3:225970/1‑149 (MQ=255)
                                                                                             aGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACcgac                                                     >  1:178172/1‑149 (MQ=255)
                                                                                             aGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACcgac                                                     >  2:39292/1‑149 (MQ=255)
                                                                                             aGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACcgac                                                     >  8:108959/1‑149 (MQ=255)
                                                                                              ggAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACcgac                                                     >  3:64051/1‑148 (MQ=255)
                                                                                                     tGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTa                                              >  2:245581/1‑148 (MQ=255)
                                                                                                     tGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTa                                              >  6:285186/1‑148 (MQ=255)
                                                                                                       gTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTaaa                                            <  5:309580/148‑1 (MQ=255)
                                                                                                       gTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAAt                                           >  2:106305/1‑149 (MQ=255)
                                                                                                       gTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAAt                                           >  8:225008/1‑149 (MQ=255)
                                                                                                          cGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTc                                         <  7:108959/148‑1 (MQ=255)
                                                                                                          cGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTc                                         <  8:198809/148‑1 (MQ=255)
                                                                                                          cGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCa                                        <  1:245581/149‑1 (MQ=255)
                                                                                                          cGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCa                                        <  2:18680/149‑1 (MQ=255)
                                                                                                           ggTTTTGACGACGGCATCTATAGCTCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAg                                       <  6:89550/149‑1 (MQ=255)
                                                                                                                   cgacgGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTg                               <  6:78979/149‑1 (MQ=255)
                                                                                                                    gacgGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGa                              >  4:75915/1‑149 (MQ=255)
                                                                                                                              ataGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCaa                    >  1:29028/1‑149 (MQ=255)
                                                                                                                               taGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGCCATTTCCACGCCGGAAATCAATATTGATGTGCCCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAAc                   >  5:324171/1‑149 (MQ=255)
                                                                                                                                 gctccgcGCGCCTGCTGGAGATCCTAAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTTCCGCATTCCCTAACGACAATTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAAcgc                 <  2:318848/146‑1 (MQ=255)
                                                                                                                                 gcgcCGCGCGCCTGCTGGAGATCCTCAGCAATACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAATTCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAAcgc                 <  6:201064/149‑1 (MQ=255)
                                                                                                                                 gcgcCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAAcgc                 <  2:314034/149‑1 (MQ=255)
                                                                                                                                         cgcCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAACGCGACGCCGa         <  7:286545/149‑1 (MQ=255)
                                                                                                                                          tcttgctgAGGATCTCCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCGCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAACGCGAc              <  2:98792/140‑1 (MQ=255)
                                                                                                                                          tcttgctgAGGATCTCCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAACGCGAc              >  1:98792/4‑143 (MQ=255)
                                                                                                                                          tcttgctgAGGATCTCCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAACGCGAc              <  2:75680/140‑1 (MQ=255)
                                                                                                                                          tcttgctgAGGATCTCCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTGGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAACGCGAc              >  1:75680/4‑143 (MQ=255)
                                                                                                                                               ctgGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCCGCATCATTGATGCACTGCAACGCGACGCCGACTggg    <  3:75915/148‑1 (MQ=255)
                                                                                                                                               ctgGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAACGCGACGCCGACTgggg   <  2:178172/149‑1 (MQ=255)
                                                                                                                                                tgGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCAGTGATGCACTGCAACGCGACGCCGACTGGGGc  >  2:273023/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GATGTGGAAAACCGGCCATTCGCTGATCAAGAAGAAGATGAAGCAGACCGGTTCGCTGCTGGCAGGCGAGATGAGCGGTCACATCTTCATCAAGGAACGCTGGTACGGTTTTGACGACGGCATCTATAGCGCCGCGCGCCTGCTGGAGATCCTCAGCAAGACCGAGCAGAGCGCCGAAAACCTGTTTGCCGCCTTCCCGAACGACATTTCCACGCCGGAAATCAATATTGATGTGACCGACGAGGGTAAATTCAGCATCATTGATGCACTGCAACGCGACGCCGACTGGGGC  >  NC_002947/6036376‑6036667

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: