Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,422,887 G→A 10.4% A14A (GCG→GCA)  ulaE → L‑xylulose 5‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,422,8870GA10.4% 136.0 / 6.4 48A14A (GCG→GCA) ulaEL‑xylulose 5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (28/15);  new base A (2/3);  total (30/18)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

CACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGA  >  NC_000913/4422741‑4423035
                                                                                                                                                  |                                                                                                                                                    
cACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCt                                                                                                                                                    <  2:160501/149‑1 (MQ=255)
 aCGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTc                                                                                                                                                   >  1:3952/1‑149 (MQ=255)
      tttATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCccgc                                                                                                                                              >  1:571079/1‑149 (MQ=255)
        tATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCccgccg                                                                                                                                            >  4:117883/1‑149 (MQ=255)
                gCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCACTCCCCGGCGGGGAgcgc                                                                                                                                    <  4:172649/149‑1 (MQ=255)
                gCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCACTCCCCGGCGGGGAgcgc                                                                                                                                    <  4:86373/149‑1 (MQ=255)
                       taTCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGctggctgg                                                                                                                             >  3:117220/1‑149 (MQ=255)
                       taTCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCACTCCCCGGCGGGGAgcgc                                                                                                                                    >  4:435227/1‑139 (MQ=255)
                       taTCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCACTCCCCGGCGGGGAgcgc                                                                                                                                    <  3:435227/142‑1 (MQ=255)
                            gTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGCAGGGCTGGCTGGAACGc                                                                                                                        >  3:67699/1‑149 (MQ=255)
                                     cTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTg                                                                                                               <  4:171998/149‑1 (MQ=255)
                                         tcCGGTGGAAGCCGCGCGTCACTTTAAACGCTCCAACGCTGAACTGTGGGGCTAAGGAGCGGATATGTCGTCCACACAAGTCCCGCTTGTCATCTATGAAAAAGCGCTCCCCGCCGGGGTGTTATGGCTGGAACGCCTGCAACTGGCaa                                                                                                           >  4:351410/1‑149 (MQ=255)
                                         tcCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCaa                                                                                                           >  1:203539/1‑149 (MQ=255)
                                                gAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGtt                                                                                                    >  1:174787/1‑149 (MQ=255)
                                                 aaGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTa                                                                                                   <  3:532410/149‑1 (MQ=255)
                                                 aaGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTa                                                                                                   <  3:211539/149‑1 (MQ=255)
                                                  aGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAg                                                                                                  <  2:302475/149‑1 (MQ=255)
                                                               ttCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGt                                                                                     >  4:500835/1‑149 (MQ=255)
                                                               ttCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGt                                                                                     <  4:315021/149‑1 (MQ=255)
                                                                   aaCGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGaa                                                                                 >  3:442395/1‑149 (MQ=255)
                                                                    aCGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGaaa                                                                                <  3:46954/149‑1 (MQ=255)
                                                                          cATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCgg                                                                          <  3:477535/149‑1 (MQ=255)
                                                                                tGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATg                                                                    <  1:116814/149‑1 (MQ=255)
                                                                                  aaCTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGGAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGaa                                                                  >  2:351565/1‑149 (MQ=255)
                                                                                  aaCTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGaa                                                                  >  1:173720/1‑149 (MQ=255)
                                                                                             tAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCg                                                       >  1:269216/1‑149 (MQ=255)
                                                                                              aaGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGc                                                      <  1:6356/149‑1 (MQ=255)
                                                                                               aGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGcc                                                     >  1:49131/1‑149 (MQ=255)
                                                                                                   gCGGATATGTTGTCCAAACAAACCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGAGGATCGCCTGTc                                                 <  3:321087/149‑1 (MQ=255)
                                                                                                                  aaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGATGGCCTGCAACTGGCAAAAACGTTAGGCTTCTATTTTGTCGAAATGTAGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAg                                  >  2:469688/1‑149 (MQ=255)
                                                                                                                  aaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAg                                  >  2:356440/1‑149 (MQ=255)
                                                                                                                  aaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAg                                  >  2:196016/1‑149 (MQ=255)
                                                                                                                  aaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAg                                  >  4:248366/1‑149 (MQ=255)
                                                                                                                  aaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAg                                  >  1:3169/1‑149 (MQ=255)
                                                                                                                         tCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGag                           >  1:544115/1‑149 (MQ=255)
                                                                                                                          cccGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGagc                          >  3:527295/1‑149 (MQ=255)
                                                                                                                               ttGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGt                     >  4:548041/1‑149 (MQ=255)
                                                                                                                                tgtcATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTc                    <  3:500835/146‑1 (MQ=255)
                                                                                                                                tGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTc                    <  4:263789/149‑1 (MQ=255)
                                                                                                                                     tCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGgcg               >  1:513124/1‑149 (MQ=255)
                                                                                                                                      cTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGgcgc              >  3:462401/1‑149 (MQ=255)
                                                                                                                                       tATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCt             >  1:437889/1‑149 (MQ=255)
                                                                                                                                       tATGAAAAAGCGCTCCCCGCCGGGGAGGGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCt             >  4:275094/1‑149 (MQ=255)
                                                                                                                                         tGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTgg           <  1:351565/149‑1 (MQ=255)
                                                                                                                                           aaaaaGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTc                                                 >  1:12218/1‑109 (MQ=255)
                                                                                                                                           aaaaaGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTc                                                 <  2:12218/109‑1 (MQ=255)
                                                                                                                                            gccaGCACTCCCCGGCGGGGAGCGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCa        >  2:172697/4‑149 (MQ=255)
                                                                                                                                            aaaaGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCa        >  2:97920/1‑149 (MQ=255)
                                                                                                                                                gcgcTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGc    <  2:269216/149‑1 (MQ=255)
                                                                                                                                                  gcTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGa  >  2:79563/1‑149 (MQ=255)
                                                                                                                                                  gcTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGa  >  4:556628/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGA  >  NC_000913/4422741‑4423035

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑