Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,268,313 T→C 9.2% I50T (ATC→ACC)  kdsA → 3‑deoxy‑D‑manno‑octulosonate 8‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,268,3130TC 9.2% 125.9 / 3.1 43I50T (ATC→ACC) kdsA3‑deoxy‑D‑manno‑octulosonate 8‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (20/19);  new base C (2/2);  total (22/21)
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

GAAACAAAAAGTGGTTAGCATTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTTCACGAAC  >  NC_000913/1268167‑1268459
                                                                                                                                                  |                                                                                                                                                  
gAAACAAAAAGTGGTTAGCATTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATc                                                                                                                                                   >  5:266137/1‑148 (MQ=255)
           tGGTTAGCATTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGtt                                                                                                                                        <  8:67249/148‑1 (MQ=255)
              ttAGCATTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCaa                                                                                                                                     <  4:113239/148‑1 (MQ=255)
                 gCATTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGc                                                                                                                                  <  5:140525/148‑1 (MQ=255)
                   aTTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCct                                                                                                                                <  7:222608/148‑1 (MQ=255)
                          aCATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGa                                                                                                                         >  6:4587/1‑148 (MQ=255)
                              cAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACa                                                                                                                       >  1:120334/1‑146 (MQ=255)
                              cAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACa                                                                                                                       >  7:262579/1‑146 (MQ=255)
                               aaCGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACa                                                                                                                       >  3:53658/1‑145 (MQ=255)
                               aaCGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAg                                                                                                                    >  4:110101/1‑148 (MQ=255)
                                aCGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACaa                                                                                                                      <  4:53658/145‑1 (MQ=255)
                                aCGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGc                                                                                                                   >  7:118084/1‑148 (MQ=255)
                                        aaTGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGct                                                                                                           <  8:262579/148‑1 (MQ=255)
                                         aTGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTCGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGctc                                                                                                          >  3:61711/1‑148 (MQ=255)
                                          tGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGctcc                                                                                                         <  3:110101/148‑1 (MQ=255)
                                                 ccGTTAGTACTGTTTGTCGGTATGAACGTGTTGTAATCTCGCGATCTTGCGATTCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCtccat                                                                                                    <  6:266137/146‑1 (MQ=255)
                                                                     tATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTACCCCTTACGTGTTCAAAGCctct                                                                                                                              <  8:256032/100‑1 (MQ=255)
                                                                     tATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTACCCCTTACGTGTTCAAAGCctct                                                                                                                              >  7:256032/1‑100 (MQ=255)
                                                                       tGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGcc                                                                             <  2:120334/147‑1 (MQ=255)
                                                                       tGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCt                                                                            <  6:118835/148‑1 (MQ=255)
                                                                        gAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCtccat                                                                                                    >  7:169331/1‑123 (MQ=255)
                                                                        gAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCtccat                                                                                                    <  8:169331/123‑1 (MQ=255)
                                                                        gAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTACCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTg                                                                           >  8:14306/1‑148 (MQ=255)
                                                                         aaCGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGCGTATCCCGTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGCCCTg                                                                           >  1:197002/1‑147 (MQ=255)
                                                                          aCGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTgg                                                                          >  1:320678/1‑147 (MQ=255)
                                                                          aCGTGTTGGAATCTCGCGATCTGGCGATGCGCATTAGCGAGCACTACGTAACTGTGACCCAGAAACTGTGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGga                                                                         <  2:101673/148‑1 (MQ=255)
                                                                            gtgtTGGAATCTCGCGATCTGGCGATGCGCATTCGCGAGCACTACGTACCTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGgaag                                                                       <  4:61711/148‑1 (MQ=255)
                                                                               ttGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCTAAGCCTCTTTTGACAAAGCCACCCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGGaagaag                                                                    >  5:95846/1‑148 (MQ=255)
                                                                                      ctcGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGt                                                                                       >  1:124217/1‑122 (MQ=255)
                                                                                      ctcGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGt                                                                                       <  2:124217/122‑1 (MQ=255)
                                                                                        cgcgATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGGAAGAAGGGATGaaaa                                                           >  5:14991/1‑148 (MQ=255)
                                                                                                 gCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTACCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTcc                                                     <  7:14306/145‑1 (MQ=255)
                                                                                                        gcATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGaa                                           <  5:268027/148‑1 (MQ=255)
                                                                                                         cATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAg                                          >  3:124667/1‑148 (MQ=255)
                                                                                                                 aGCACTAAGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCAACCCGCTCCTCCATCCACTCTTATCGTGAAACTGGCCTGGAAGAAGGGATGAAAATCATCCAGGAGTAGAAGCAGACttt                                  >  4:157005/1‑148 (MQ=255)
                                                                                                                    aCTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTgg                               >  8:180713/1‑148 (MQ=255)
                                                                                                                      tACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAAt                                                          >  4:143858/1‑119 (MQ=255)
                                                                                                                      tACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAAt                                                          <  3:143858/119‑1 (MQ=255)
                                                                                                                                     cAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAgg                                                  >  4:66602/1‑112 (MQ=255)
                                                                                                                                     cAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAgg                                                  <  3:66602/112‑1 (MQ=255)
                                                                                                                                         aaCTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGt          <  8:118084/148‑1 (MQ=255)
                                                                                                                                           cTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTTc        <  2:320678/148‑1 (MQ=255)
                                                                                                                                                tATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTTCACGaa   >  1:324968/1‑148 (MQ=255)
                                                                                                                                                  tCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTTCACGAAc  >  4:35186/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GAAACAAAAAGTGGTTAGCATTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTTCACGAAC  >  NC_000913/1268167‑1268459

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