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,6650TG95.5% 67.6 / ‑3.2 22L206W (TTG→TGG) nagKN‑acetylglucosamine kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (11/10);  minor base A (1/0);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

GGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGC  >  NZ_CP009273/1174586‑1174749
                                                                               |                                                                                    
ggTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTccc                                                                            <  1:37337/90‑1 (MQ=255)
            tGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTgcg                                                                <  1:599787/90‑1 (MQ=255)
               aTTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTTCCGAAATTATTGGGCtt                                                             >  1:428462/1‑90 (MQ=255)
                         aTCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAg                                                   >  1:483544/1‑90 (MQ=255)
                            tgtctgGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCg                                                >  1:169877/1‑90 (MQ=255)
                                 ggTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAg                                           <  2:613/90‑1 (MQ=255)
                                 ggTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAg                                           <  2:171604/90‑1 (MQ=255)
                                 ggTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAg                                           <  2:169877/90‑1 (MQ=255)
                                 ggTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAg                                           >  2:501893/1‑90 (MQ=255)
                                      cgGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGc                                      <  2:352406/90‑1 (MQ=255)
                                       ggTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCa                                     >  1:262071/1‑90 (MQ=255)
                                        gTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCaa                                    <  2:370931/90‑1 (MQ=255)
                                                         cAACACTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGt                   <  1:537742/90‑1 (MQ=255)
                                                            cacTATTATCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTAt                >  1:465219/1‑90 (MQ=255)
                                                             acTATTATCATCAACCGTAGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATc               >  2:155561/1‑90 (MQ=255)
                                                               tattatCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTg             >  2:44523/1‑90 (MQ=255)
                                                                attatCATCAACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTgg            >  2:119377/1‑90 (MQ=255)
                                                                        caACCGTGGCAGGCTCCCGAAATTATTGCGCTTTa                                                           >  1:663584/1‑35 (MQ=38)
                                                                        caACCGTGGCAGGCTCCCGAAATTATTGCGCTTTa                                                           <  2:663584/35‑1 (MQ=38)
                                                                        caACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCGAGGGCGCACGTTGAGCGTTATCTGGATttatta    >  2:419241/1‑90 (MQ=255)
                                                                        caACCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATttatta    <  1:740530/90‑1 (MQ=255)
                                                                          aCCGTGGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGc  >  1:385329/1‑90 (MQ=255)
                                                                               |                                                                                    
GGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGC  >  NZ_CP009273/1174586‑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: