Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,698,3020AC55.6% 10.1 / 12.4 18G88G (GGT→GGGydgJoxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (7/1);  new base C (0/10);  total (7/11)
Fisher's exact test for biased strand distribution p-value = 2.51e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-03
Rejected as consensus: Frequency below/above cutoff threshold.

CCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGG  >  NZ_CP009273/1698221‑1698379
                                                                                  |                                                                             
ccccAGGCTTTTTGCCATCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTc                                                                       >  1:987004/1‑90 (MQ=255)
      gttttttGCCGGCCCACCCAGCTCCCGCGTTTGTGACGGTGCCAGGGTAAGGGTTTTATCAACGCCCCCTTTTTTCCCCGCCTCAAgcgc                                                                 <  1:408695/88‑1 (MQ=255)
      gttttttCCCGGCGCTCCCACCTCCCGCGTTTGTGACAGTGTCAGGGTAAGGGTTTTACAACCGCCCACTTTTTTACCCGCCTCAAgcgc                                                                 <  1:449407/88‑1 (MQ=255)
          tttGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCt                                                             >  2:934752/1‑90 (MQ=255)
                aGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCt                                                       >  2:816094/1‑90 (MQ=255)
                     aTCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACgg                                                  >  2:633607/1‑90 (MQ=255)
                        cgcctctCGCGTTTGTGACAGTGTCAGGGTAAGGGGTTTATCACCGCCCACTTTTTTCCCCGCCTCAAGCGCCGCTTTGGCTAACGGGaa                                               <  1:403237/87‑1 (MQ=255)
                                    ttGTGACAGTGCCCGGGTAAGGGTTTTAACAACCCCCACTTTTTTCCCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGtt                                   <  2:987004/90‑1 (MQ=255)
                                     tgtgACGGTGTCAGGGAAAGGGTTTTACAAACGACCACTTTTTTCCCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTg                                  <  1:197230/90‑1 (MQ=255)
                                     tgggACGGTGTCAGGGAAAGGGTTTTACCAACGCCCAATTTTTTCCCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTg                                  <  1:591447/87‑1 (MQ=255)
                                     tgggACAGTGCCGGGGAAAGGGTTTTAACAACGACCATTTTTTTCCCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTg                                  <  2:993081/87‑1 (MQ=255)
                                       cgACGGTGCCACGATAAGGGTTTTACCAACGTCCATATTTTTCCCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTggg                                <  1:633607/89‑1 (MQ=255)
                                            gtgtCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAg                           >  1:590897/1‑90 (MQ=255)
                                                   ggaaaaGGTTTTACAAACGACCATTTTTTTCCCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGAc                    <  2:1029134/87‑1 (MQ=255)
                                                       aaGGGTTTTACCAACAACCATTTTTTTCCCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAAt                 <  2:590897/90‑1 (MQ=255)
                                                             tttATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTc           >  1:340901/1‑90 (MQ=255)
                                                               tATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTa         >  1:914499/1‑90 (MQ=255)
                                                                      ccccAAATTTTTCCCCGCCTCAAGCGCCGCTTTGGATAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTggg  <  1:402787/88‑1 (MQ=255)
                                                                                  |                                                                             
CCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGG  >  NZ_CP009273/1698221‑1698379

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