Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,436,0610TG55.6% 14.6 / 24.0 18G273G (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 (0/8);  new base G (10/0);  total (10/8)
Fisher's exact test for biased strand distribution p-value = 2.29e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTGGG  >  NZ_CP009273/2435974‑2436147
                                                                                       |                                                                                      
ggAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTAt                                                                                      <  1:256054/90‑1 (MQ=255)
    cAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGcc                                                                                  <  1:582440/90‑1 (MQ=255)
     aGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGGTTGccc                                                                                 >  2:380636/1‑89 (MQ=255)
     aGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGCCa                                                                                 >  1:119723/1‑90 (MQ=255)
           tACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTTCCcgcgcgg                                                                           >  1:50818/1‑89 (MQ=255)
                cTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCCGCGCGttggtg                                                                      >  2:379293/1‑90 (MQ=255)
                    ccTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTTCCCGCGCGGTGTTGTCGc                                                                  >  1:73870/1‑90 (MQ=255)
                          ccGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGgcgc                                                            <  1:449054/90‑1 (MQ=255)
                          ccGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGgcgc                                                            <  1:467261/90‑1 (MQ=255)
                          ccGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTTGCCGCGCGGGGTTTTCGCCTgcgc                                                            >  1:570062/1‑90 (MQ=255)
                            gcgcCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTTCCCGCGCGGTGTTGTTGCTGgcgcct                                                          >  1:133255/1‑88 (MQ=255)
                            gcgcCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGGGTGCCAGCGCGGTGTTGTTGGTGGCGCTa                                                          >  2:467509/1‑90 (MQ=255)
                            gcgcCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGCCAGCGCGTTGTTGTCGCTGGCGCTa                                                          >  2:279344/1‑90 (MQ=255)
                                         aCCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGCCAGCGCGTTGTTGTCGCTGGGGCTATTACGGCGCGGCt                                             >  2:478848/1‑90 (MQ=255)
                                                  gcagcagcaAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTaa                                    <  1:126995/90‑1 (MQ=255)
                                                               ggAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTgcgc                       <  1:699523/90‑1 (MQ=255)
                                                               ggAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTgcgc                       <  1:498016/90‑1 (MQ=255)
                                                                           tATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGcccc           <  2:424257/90‑1 (MQ=255)
                                                                                    tggtATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTggg  <  2:694218/90‑1 (MQ=255)
                                                                                       |                                                                                      
GGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTGGG  >  NZ_CP009273/2435974‑2436147

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