Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273984,3600TG52.6% 21.8 / 25.1 19T28P (ACC→CCC) asnSasparagine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (1/8);  new base G (10/0);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 1.19e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCCCTGGAGTACGTCGGCTACAGGCACAACGCTCATA  >  NZ_CP009273/984291‑984442
                                                                     |                                                                                  
gTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGGACGTACCCCTCCCCGCACgg                                                                >  2:14234/1‑90 (MQ=255)
    aGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGGACGCAACCATGCACCCACgggggc                                                            >  1:447757/1‑86 (MQ=255)
      gATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGGACGGACCCATCCGCGCCCCGGGACtt                                                          >  2:676137/1‑90 (MQ=255)
         cAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGc                                                       <  1:854550/90‑1 (MQ=255)
          aaaGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGGACGTACCCATCCCCGCACGGGGGCTTTGCt                                                      >  2:773017/1‑90 (MQ=255)
          aaaGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGGACGTACCCATCCCCGCACGGGGACTTCGCt                                                      >  1:155188/1‑90 (MQ=255)
          aaaGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGGACGGACCCATCCCCGCACGGGGACTTCGCt                                                      >  1:368105/1‑90 (MQ=255)
          aaaGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGGACCGACCCCCCCCCCCACGGGGGCCTCGCt                                                      >  1:622396/1‑90 (MQ=255)
                ggAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGGACGGACCCAACCGCCCACGGGGACTTCGCCGTCAAc                                                >  2:864685/1‑90 (MQ=255)
                     cGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACCGACCCCCCCCCGCCCGGGGGCTTATCTGTCAAAGGcga                                           >  1:726780/1‑88 (MQ=255)
                      gTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGGACGGACCCCTCCGCGCACGGGGACTTCGCTGTCAACGGCTac                                          >  2:901003/1‑90 (MQ=255)
                        cATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTacac                                        <  2:169494/90‑1 (MQ=255)
                              cGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCCCt                                  <  1:68294/90‑1 (MQ=255)
                                gCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGGACGTACCCCTCCCCGCACGGGGACTTCGCTGTCAACGGGCACACCGCCCCgg                                >  1:266994/1‑90 (MQ=255)
                                                     tttGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCCCTGGAGTACGTCGGCTACAGGcaca           <  1:901003/90‑1 (MQ=255)
                                                       tGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCCCTGGAGTACGTCGGCTACAGGCACAAc         <  2:447757/90‑1 (MQ=255)
                                                        gAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCCCTGGAGTACGTCGGCTACAGGCACAACg        <  2:89324/90‑1 (MQ=255)
                                                          aTCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCCCTGGAGTACGTCGGCTACAGGCACAACGCt      <  1:951373/90‑1 (MQ=255)
                                                              cggcggGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCCCTGGAGTACGTCGGCTACAGGCACAACGTTCata  <  1:31177/90‑1 (MQ=255)
                                                                     |                                                                                  
GTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAGCTTTTGAATCTCGGCGGGTACGTACCCATCCGCGCACGGTGACTTCGCTGTCAACGGCTACACGGCCCTGGAGTACGTCGGCTACAGGCACAACGCTCATA  >  NZ_CP009273/984291‑984442

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