Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,422,887 G→A 9.3% 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,8870GA 9.3% 120.9 / 3.5 43A14A (GCG→GCA) ulaEL‑xylulose 5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (22/17);  new base A (2/2);  total (24/19)
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.27e-01

TCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGT  >  NC_000913/4422751‑4423027
                                                                                                                                        |                                                                                                                                            
tCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCggg                                                                                                                                  <  1:25253/149‑1 (MQ=255)
          tagtaTCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCACTCCCCGGCGGGGAGCgcttt                                                                                                                         <  4:325892/148‑3 (MQ=255)
           agtaTCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTgctggc                                                                                                                        >  3:183435/1‑148 (MQ=255)
            gtaTCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTgctggctg                                                                                                                      >  2:276704/1‑149 (MQ=255)
             taTCCGTGATGCCGCTTCGCCGGGGGAAGCCGCACTCCAGGTCAAACGTGCCATCGCTTAACTGGGGGGCTAGGGCGCGGTTGGGCTGTGCAAACAAATCCCTCCTGGCATCGAGCGAAAAGTGCTGCCCGCCGGGGAGCgctggctg                                                                                                                      >  4:117991/1‑148 (MQ=255)
                ccGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCACTCCCCGGCGGGGAGCgctttttcat                                                                                                                    <  4:157088/147‑8 (MQ=255)
                ccGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCACTCCCCGGCGGGGAGCgctttttcat                                                                                                                    >  3:157088/1‑140 (MQ=255)
                  gTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACg                                                                                                                 >  4:167185/1‑148 (MQ=255)
                     aTGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCccgccg                                                                                                                                    <  1:318378/126‑1 (MQ=255)
                     aTGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCccgccg                                                                                                                                    >  2:318378/1‑126 (MQ=255)
                      tGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTg                                                                                                             <  3:138914/148‑1 (MQ=255)
                      tGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTg                                                                                                             >  3:401255/1‑148 (MQ=255)
                                  ggTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCaaaa                                                                                                 <  3:64489/148‑1 (MQ=255)
                                  ggTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCaaaa                                                                                                 <  3:172139/148‑1 (MQ=255)
                                      gAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGt                                                                                             >  3:237997/1‑148 (MQ=255)
                                             cACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCtt                                                                                      >  4:359449/1‑148 (MQ=255)
                                                    gTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTg                                                                              <  2:393112/149‑1 (MQ=255)
                                                      tCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGAtttt                                                                               <  4:309887/146‑1 (MQ=255)
                                                      tCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGAtttt                                                                               >  3:309887/1‑146 (MQ=255)
                                                      tCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGt                                                                             <  4:237997/148‑1 (MQ=255)
                                                      tCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGt                                                                             >  3:187642/1‑148 (MQ=255)
                                                      tCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTc                                                                            <  2:16333/149‑1 (MQ=255)
                                                                 aTCGCTGAACTGTGGGGCTAAGGAGCGGATATGGTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGt                                                                 >  1:285557/1‑149 (MQ=255)
                                                                  tCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGt                                                                 >  3:313680/1‑148 (MQ=255)
                                                                  tCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGt                                                                 >  3:261499/1‑148 (MQ=255)
                                                                  tCGCTGAACTGTGGGGCTAAGGAGCGCATATGTTGTCCAAGCAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTa                                                                <  1:38680/149‑1 (MQ=255)
                                                                     cTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGa                                                              <  3:210605/148‑1 (MQ=255)
                                                                      tGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATg                                                            >  1:204570/1‑149 (MQ=255)
                                                                      tGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCACTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATg                                                            >  2:278097/1‑149 (MQ=255)
                                                                              ggggCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGa                                                     <  3:167185/148‑1 (MQ=255)
                                                                                    aaGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCg                                               <  4:313680/148‑1 (MQ=255)
                                                                                                tgttgtCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCt                                   <  3:309544/148‑1 (MQ=255)
                                                                                                   tgtCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCt                                                                                       >  4:379108/1‑93 (MQ=255)
                                                                                                   tgtCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCt                                                                                       <  3:379108/93‑1 (MQ=255)
                                                                                                       caaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGa                           >  1:391924/1‑149 (MQ=255)
                                                                                                        aaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAg                          >  1:250038/1‑149 (MQ=255)
                                                                                                         aacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAg                          >  4:6411/1‑148 (MQ=255)
                                                                                                              aTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGa                    <  2:285557/149‑1 (MQ=255)
                                                                                                                cccGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTc                                                                            <  2:43421/91‑1 (MQ=255)
                                                                                                                cccGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTc                                                                            >  1:43421/1‑91 (MQ=255)
                                                                                                                     ttGGCATCTATGAAACAGCGCTCCCCGCCGGGTAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACTTTCGGCTTCGATTTTGTTGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGGGCAGCGt             >  2:310647/1‑149 (MQ=255)
                                                                                                                             tATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCt     >  2:320744/1‑149 (MQ=255)
                                                                                                                                gAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGt  >  1:401277/1‑149 (MQ=255)
                                                                                                                                        |                                                                                                                                            
TCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGT  >  NC_000913/4422751‑4423027

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