Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,413,422 A→G S254G (AGC→GGC)  ulaE → L‑ribulose‑5‑phosphate 3‑epimerase UlaE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,413,4220AG90.9% 60.2 / ‑0.7 22S254G (AGC→GGC) ulaEL‑ribulose‑5‑phosphate 3‑epimerase UlaE
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (15/5);  total (16/6)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

GGTGTAGTGGATTTCGAACGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGCCGTACCTGATTGAGATGTGGAGCGAAACGGCGGAAGACCCGGCGGCAGAAGTGGCGAAAGCGCGTGATTGGGTGAAAGCGCGCATGGCGAAAGCGGGCATGGTGGAG  >  NZ_CP009273/4413344‑4413508
                                                                              |                                                                                      
ggTGTAGTGGATTTCGAACGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGCCGTACCTGATTGAGATGTGGGGCGAAAcggcg                                                                             >  2:172846/1‑90 (MQ=255)
          aTTTCGAACGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGCCGTACCTGATTGAGATGTGGAGCGAAACGGCGGAAGACCcgg                                                                   >  1:32709/1‑90 (MQ=255)
           tttCGAACGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGCCGTTCCTGATTGAGATGGGGGGCGAAACGGCGGAAGACCcggc                                                                  >  2:291880/1‑90 (MQ=255)
           tttCGAACGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGCCGTACCTGATTGAGATGTGGGGCGAAACGGCGGAAGACCcggc                                                                  >  1:210165/1‑90 (MQ=255)
           tttCGAACGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGCCGTACCTGATTGAGATGTGGGGCGAAACGGCGGAAGACCcggc                                                                  >  1:3553/1‑90 (MQ=255)
              cGAACGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGCCGTACCTGATTGAGATGTGGGGCGAAACGGCGGAAGACCCggcggc                                                               <  2:280691/90‑1 (MQ=255)
                  cGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGCCGTACCTGATTGAGATGTGGGGCGAAACGGCGGAAGACCCGGCGGCAGaa                                                           >  1:465323/1‑90 (MQ=255)
                   gttgttTCGAAACGCTCAAACAGAGTGGCTATTGCGGGCCGTACCTGATTGAGATGTGGAGCGAAACGGCGGAAGACCCGGCGGCAGAAg                                                          <  2:32709/90‑1 (MQ=255)
                           gAAACGCTCAAACAGAGTGGCTATTGCGGGCCGTACCTGATTGAGATGTGGGGCGAAACGGCGGAAGACCCGGCGGCAGAAGTGGCGaaa                                                  <  1:457639/90‑1 (MQ=255)
                                gCTCAAACAGAGTGGCTATTGCGGGCCGTACCTGATTGAGATGTGGGGCGAAACGGCGGAAGACCCGGCGGCAGAAGTGGCGAAAGcgcg                                             >  2:491035/1‑90 (MQ=255)
                                  tCAAACAGAGTGGCTATTGCGGGCCGTACCTGATTGAGATGTGGGGCGAAACGGCGGAAGACCCGGCGGCAGAAGTGGCGAAAGCGCGTg                                           <  1:391828/90‑1 (MQ=255)
                                           gTGGCTATTGCGGGCCGTACCTGATTGAGATGTGGGGCGAAACGGCGGAAGACCCGGCGGCAGAAGTGGCGAAAGCGCGTGATTGGGTGa                                  >  1:280166/1‑90 (MQ=255)
                                                   tGCGGGCCGTACCTGATTGAGATGTGGGGCGAAACGGCGGAAGACCCGGCGGCAGAAGTGGCGAAAGCGCGGGATTGGGGGAAAgcgcgc                          >  2:212543/1‑90 (MQ=255)
                                                           gTACCTGATTGAGATGTGGGGCGAAACGGCGGAAGACCCGGCGGCAGAAGTGGCGAAAGCGCGTGATTGGGTGAAAGCGCGCATGGCGaa                  >  2:181406/1‑90 (MQ=255)
                                                             aCCTGATTGAGATGTGGGGCGAAACGGCGGAAGACCCGGCGGCAGAAGTGGCGAAAGCGCGTGATTGGGTGAAAGCGCGCATGGCGAAAg                <  1:187720/90‑1 (MQ=255)
                                                                 gATTGAGATGTGGGGCGAAACGGCGGAAGACCCGGCGGCAGAAGTGGCGAAAGCGCGTGATTGGGTGAAAGCGCGCATGGCGAAAGCggg            >  1:35656/1‑90 (MQ=255)
                                                                 gATTGAGATGTGGGGCGAAACGGCGGAAGACCCGGCGGCAGAAGTGGCGAAAGCGCGTGATTGGGTGAAAGCGCGCATGGCGAAAGCggg            >  2:407273/1‑90 (MQ=255)
                                                                     gagaTGTGGGGCGAAACGGCGGAAGACCCGGCGGCAGGAGTGGCGAAAGCGCGTGATTGGGGGAAAGGGCGCGTGGCGAAAGCGGGCAtg        >  1:262770/1‑90 (MQ=255)
                                                                     gagaTGTGGGGCGAAACGGCGGAAGACCCGGCGGCAGAAGTGGCGAAAGCGCGTGATTGGGTGAAAgcgcgc                          >  1:348491/1‑72 (MQ=255)
                                                                     gagaTGTGGGGCGAAACGGCGGAAGACCCGGCGGCAGAAGTGGCGAAAGCGCGTGATTGGGTGAAAgcgcgc                          <  2:348491/72‑1 (MQ=255)
                                                                     gagaTGTGGGGCGAAACGGCGGAAGACCCGGCGGCAGAAGTGGCGAAAGCGCGTGATTGGGTGAAAGCGCGCATGGCGAAAGCGGGCAtg        >  1:44001/1‑90 (MQ=255)
                                                                           tgGGGCGAAACGGCGGAAGACCCGGCGGCAGAAGTGGCGAAAGCGCGTGATTGGGTGAAAGCGCGCATGGCGAAAGCGGGCATGGTGGAg  >  1:263096/1‑90 (MQ=255)
                                                                              |                                                                                      
GGTGTAGTGGATTTCGAACGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGCCGTACCTGATTGAGATGTGGAGCGAAACGGCGGAAGACCCGGCGGCAGAAGTGGCGAAAGCGCGTGATTGGGTGAAAGCGCGCATGGCGAAAGCGGGCATGGTGGAG  >  NZ_CP009273/4413344‑4413508

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑