Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,436,0580TG52.9% 14.4 / 18.1 17G272G (GGT→GGGmnmCbifunctional tRNA (5‑methylaminomethyl‑2‑thiouridine)(34)‑methyltransferase MnmD/FAD‑dependent 5‑carboxymethylaminomethyl‑2‑thiouridine(34) oxidoreductase MnmC
Reads supporting (aligned to +/- strand):  ref base T (2/6);  new base G (9/0);  total (11/6)
Fisher's exact test for biased strand distribution p-value = 2.26e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-03
Rejected as consensus: Frequency below/above cutoff threshold.

TACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTGG  >  NZ_CP009273/2435985‑2436146
                                                                         |                                                                                        
tACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTTCCcgcgcgg                                                                          >  1:531769/1‑89 (MQ=255)
    gCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGCCCGCGCgttgtt                                                                      >  1:771632/1‑90 (MQ=255)
     cTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGTATTGCCAGCGCggtgttt                                                                     >  1:50036/1‑89 (MQ=255)
     cTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGCCAGCGCggtgttt                                                                     >  1:29882/1‑89 (MQ=255)
       ccccTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGCCAGCGCGGTGTTTTc                                                                   >  1:968109/1‑90 (MQ=255)
            gCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGCCCGCGCGGTGTTTTCGCTgg                                                              >  2:57309/1‑90 (MQ=255)
                 gcgcCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTTCCCGCGCGGTGTTGTTGCTGGGGCTa                                                         >  1:635395/1‑90 (MQ=255)
                 gcgcCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTTCCCGCGCGGTGGTGTTGCTGGCGCTa                                                         >  2:902824/1‑90 (MQ=255)
                 gcgcCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGCCAGCGCGGTGTTGTCGCTGGCGCTa                                                         >  1:318400/1‑90 (MQ=255)
                                   aCGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTTCCAGCGCGTTGTTTTCGCTGGGGCTATTACGGCGCGGGTGGGAg                                       >  1:973490/1‑90 (MQ=255)
                                       gcagcagcaAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTaa                                   <  1:450060/90‑1 (MQ=255)
                                                aaCGGGAAGCGGCGATTATCGGCGGTGGGATTGCCCGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGGAGGGAAAGCCTtttt                          >  1:914096/1‑90 (MQ=255)
                                                    ggAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTgcgc                      <  1:243680/90‑1 (MQ=255)
                                                         cggcgATTATCGGCGGTGGTATTGCCAGTGCGTTGTTGTCGCTGGCGCTATTCCGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATg                 <  2:973490/90‑1 (MQ=255)
                                                           gcgATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAg               <  2:235086/90‑1 (MQ=255)
                                                                tATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGcccc          <  1:542503/90‑1 (MQ=255)
                                                                     gcggTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCAc     <  2:968109/90‑1 (MQ=255)
                                                                        gtggtATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTgg  <  2:976755/90‑1 (MQ=255)
                                                                         |                                                                                        
TACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTGG  >  NZ_CP009273/2435985‑2436146

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