Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,174,665 T→G L206W (TTG→TGG)  nagK → N‑acetylglucosamine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,174,6650TG98.5% 236.3 / ‑3.6 65L206W (TTG→TGG) nagKN‑acetylglucosamine kinase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (39/25);  total (40/25)
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.83e-01

GGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGC  >  NZ_CP009273/1174580‑1174749
                                                                                     |                                                                                    
ggCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAg                                                                                  <  2:215816/90‑1 (MQ=255)
   tgtgGTCAGCGTGGCTGCATTGAAAGTTATCTGTCTGGGCGGGGGTTTGCGTGGCGGTATCAACACTGTTCTGATCAACCGTGGCAgtcg                                                                               >  1:1806272/1‑87 (MQ=255)
      ggTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGc                                                                                >  2:712488/1‑86 (MQ=255)
      ggTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGc                                                                                <  1:712488/86‑1 (MQ=255)
      ggTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTccc                                                                            <  1:2225845/90‑1 (MQ=255)
      ggTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTccc                                                                            >  2:1804135/1‑90 (MQ=255)
              tGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAattat                                                                    <  2:1164794/90‑1 (MQ=255)
                gCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTg                                                                  >  2:1943564/1‑90 (MQ=255)
                gCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTg                                                                  <  1:204092/90‑1 (MQ=255)
                 cTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTgc                                                                 <  1:2169067/90‑1 (MQ=255)
                     aTTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGGGCtt                                                             >  2:3006538/1‑90 (MQ=255)
                     aTTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCtt                                                             >  1:421418/1‑90 (MQ=255)
                     aTTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCtt                                                             >  2:2569224/1‑90 (MQ=255)
                     aTTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCtt                                                             >  2:152189/1‑90 (MQ=255)
                     aTTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCtt                                                             >  2:1889966/1‑90 (MQ=255)
                         aaaaTTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCCGGCTCCCGAAATTATTGCGCTTTATg                                                         >  1:416657/1‑90 (MQ=255)
                         aaaaTTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATg                                                         <  2:2270418/90‑1 (MQ=255)
                             ttATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCa                                                     <  1:3006538/90‑1 (MQ=255)
                              tATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCaa                                                    <  1:36132/90‑1 (MQ=255)
                                   gtctgGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCa                                                     <  2:2711168/84‑1 (MQ=255)
                                   gtctgGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCa                                                     >  1:684663/1‑84 (MQ=255)
                                   gtctgGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCa                                                     <  2:684663/84‑1 (MQ=255)
                                   gtctgGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCa                                                     >  1:2711168/1‑84 (MQ=255)
                                   gtctgGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGa                                               <  2:620980/90‑1 (MQ=255)
                                     ctgGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATg                                             >  1:2396525/1‑90 (MQ=255)
                                     ctgGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATg                                             >  2:2913602/1‑90 (MQ=255)
                                     ctgGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATg                                             <  2:2982758/90‑1 (MQ=255)
                                     ctgGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATg                                             <  1:1759056/90‑1 (MQ=255)
                                      tgGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGa                                            >  2:158200/1‑90 (MQ=255)
                                        gTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGc                                          <  2:2956124/90‑1 (MQ=255)
                                        gTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGc                                          <  1:1943564/90‑1 (MQ=255)
                                        gTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGc                                          <  1:81675/90‑1 (MQ=255)
                                                  tGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGTAAGGgc                                >  2:850416/1‑90 (MQ=255)
                                                    cGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGgcgc                              <  2:1152907/90‑1 (MQ=255)
                                                    cGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGgcgc                              <  2:1990235/90‑1 (MQ=255)
                                                                  cacTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTAt                >  1:1315411/1‑90 (MQ=255)
                                                                   acTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATc               <  1:806612/90‑1 (MQ=255)
                                                                     tattatGATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTg             >  2:2274761/1‑90 (MQ=255)
                                                                     tattatCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTg             >  2:2438777/1‑90 (MQ=255)
                                                                     tattatCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTg             >  2:1900384/1‑90 (MQ=255)
                                                                     tattatCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTg             >  1:225222/1‑90 (MQ=255)
                                                                      attatCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTgg            >  2:606425/1‑90 (MQ=255)
                                                                      attatCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTgg            >  2:2693130/1‑90 (MQ=255)
                                                                       ttatCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGa           >  1:101576/1‑90 (MQ=255)
                                                                        tatCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGGTCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGAt          <  1:1073112/90‑1 (MQ=255)
                                                                        tatCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGGTCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGAt          <  1:1073113/90‑1 (MQ=255)
                                                                        tatCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGAt          <  1:1889966/90‑1 (MQ=255)
                                                                           catcaACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATtta       >  1:846839/1‑90 (MQ=255)
                                                                           catcaACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATtta       >  1:1903565/1‑90 (MQ=255)
                                                                            atcaACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGc                                      <  2:1537799/58‑1 (MQ=255)
                                                                            atcaACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGc                                      >  1:1537799/1‑58 (MQ=255)
                                                                            atcaACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGAtttgt      >  2:967997/1‑88 (MQ=255)
                                                                             tcaACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATttatt     >  2:542759/1‑90 (MQ=255)
                                                                              caACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGTGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATttatta    >  1:1281299/1‑90 (MQ=255)
                                                                              caACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATttatta    >  1:600995/1‑90 (MQ=255)
                                                                              caACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATttatta    >  2:2578730/1‑90 (MQ=255)
                                                                              caACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATttatta    >  1:1509391/1‑90 (MQ=255)
                                                                              caACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATttatta    >  2:877489/1‑90 (MQ=255)
                                                                                aCCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAg                                           <  2:2829541/49‑1 (MQ=255)
                                                                                aCCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAg                                           >  1:2829541/1‑49 (MQ=255)
                                                                                aCCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGc  >  1:225656/1‑90 (MQ=255)
                                                                                aCCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGc  >  2:377927/1‑90 (MQ=255)
                                                                                aCCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGc  >  1:1585203/1‑90 (MQ=255)
                                                                                aCCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGc  >  1:86093/1‑90 (MQ=255)
                                                                                aCCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGc  >  2:1565486/1‑90 (MQ=255)
                                                                                     |                                                                                    
GGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGC  >  NZ_CP009273/1174580‑1174749

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