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,6650TG94.1% 103.1 / ‑0.8 34L206W (TTG→TGG) nagKN‑acetylglucosamine kinase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (21/11);  total (22/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-01

GCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGC  >  NZ_CP009273/1174596‑1174749
                                                                     |                                                                                    
gCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTg                                                                  >  1:445067/1‑90 (MQ=255)
 cTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTgc                                                                 >  2:807471/1‑90 (MQ=255)
  tGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTgcg                                                                <  2:985806/90‑1 (MQ=255)
   gCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTgcgc                                                               >  2:80418/1‑90 (MQ=255)
   gCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTgcgc                                                               >  2:530656/1‑90 (MQ=255)
     aTTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCtt                                                             >  2:1357232/1‑90 (MQ=255)
        gAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTAt                                                          >  1:37738/1‑90 (MQ=255)
          aaaTTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATCATCATCAACCGTGGCAGGc                                                                                >  1:1561363/1‑66 (MQ=255)
          aaaTTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATCATCATCAACCGTGGCAGGc                                                                                <  2:1561363/66‑1 (MQ=255)
             ttATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCa                                                     >  1:1408058/1‑90 (MQ=255)
                  tgtctgGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCg                                                >  1:1545041/1‑90 (MQ=255)
                  tgtctgGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCg                                                >  2:206131/1‑90 (MQ=255)
                   gtctgGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGa                                               <  1:807471/90‑1 (MQ=255)
                     ctgGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGc                                                                                <  1:208067/55‑1 (MQ=255)
                     ctgGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGc                                                                                >  2:208067/1‑55 (MQ=255)
                       ggTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATTAAGGCGATGAg                                           >  2:1413892/1‑90 (MQ=255)
                       ggTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATTAAGGCGATGAg                                           <  1:1413892/90‑1 (MQ=255)
                       ggTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAg                                           <  2:1497222/90‑1 (MQ=255)
                       ggTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAg                                           <  2:51910/90‑1 (MQ=255)
                       ggTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAg                                           <  1:1544725/90‑1 (MQ=255)
                        gTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGc                                          <  2:1485889/90‑1 (MQ=255)
                              gTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCaa                                    <  2:200093/90‑1 (MQ=255)
                                tttGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAgg                                  <  2:901907/90‑1 (MQ=255)
                                  tGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGgc                                >  2:571522/1‑90 (MQ=255)
                                                  cacTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTAt                >  1:28726/1‑90 (MQ=255)
                                                     tattatCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTg             >  1:10178/1‑90 (MQ=255)
                                                     tattatCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTg             >  1:1119190/1‑90 (MQ=255)
                                                     tattatCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTg             >  1:744367/1‑90 (MQ=255)
                                                              caACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGgcg                               >  2:834972/1‑63 (MQ=255)
                                                              caACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGgcg                               <  1:834972/63‑1 (MQ=255)
                                                              caACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATttatta    >  2:594244/1‑90 (MQ=255)
                                                              caACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATttatta    >  1:1580097/1‑90 (MQ=255)
                                                              caACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATttatta    >  1:895174/1‑90 (MQ=255)
                                                                aCCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGc  >  2:619431/1‑90 (MQ=255)
                                                                     |                                                                                    
GCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGC  >  NZ_CP009273/1174596‑1174749

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