Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273296,4220GC52.9% 12.0 / 26.7 17P75R (CCG→CGG) paoCaldehyde oxidoreductase molybdenum‑binding subunit PaoC
Reads supporting (aligned to +/- strand):  ref base G (3/5);  new base C (9/0);  total (12/5)
Fisher's exact test for biased strand distribution p-value = 9.05e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGG  >  NZ_CP009273/296335‑296511
                                                                                       |                                                                                         
gTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAAAGACAACCCGTACGCCCCgc                                                                                         >  1:260433/1‑90 (MQ=255)
gTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGTGAGCCCCGACGCCCgcc                                                                                         >  1:790670/1‑88 (MQ=255)
  ggggCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGGAATACAACCCGGACGCCCCgcgc                                                                                       >  2:796524/1‑90 (MQ=255)
     gccgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCttt                                                                                    >  1:175101/1‑90 (MQ=255)
                gCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAAGCCGGACCCCCCGCGCCTTTTTCtcggcgg                                                                         >  2:387724/1‑89 (MQ=255)
                                  cGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAAGGACAGCCAGTAAGCCCCGCGCTTTTTTGGCGGGGTCCGGCTCGAGGGGGGTgg                                                       >  2:652934/1‑89 (MQ=255)
                                     cTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGc                                                    <  1:652934/90‑1 (MQ=255)
                                       ttGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTTCGCGGCGTCCCTATCAAAGGGGGGGGGGCgg                                                  >  1:185072/1‑89 (MQ=255)
                                           cGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCtt                                              <  2:704525/90‑1 (MQ=255)
                                            gAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCCGTACGCCCCGCGCTTTTTTCGCGGGGTCCGTATCAAAGGGGGGGGGGGGGCCttt                                             >  2:272950/1‑90 (MQ=255)
                                              gTGCCCCGGCGTTACTGGCGGGAAAGACAACCAGTACGCCCCGCGCCTTTTTTGCGGCGGCCGTATCAAGGGGGGGGAGGCGGCCctttg                                           >  2:785704/1‑88 (MQ=255)
                                                                cGGTAATGACAGCCAGTACGCCCCGCGCTTTTTGCGCGGCGGCCGTATCAAGGGGGGGGGGGGGGCCTTTGGCAAAGGCGGGAAccacga                         >  2:147994/1‑90 (MQ=255)
                                                                          aGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACCATATAGCCATa               >  2:535440/1‑90 (MQ=255)
                                                                           gCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAg              >  2:691110/1‑90 (MQ=255)
                                                                            ccAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCg                                                   <  1:74363/52‑1 (MQ=255)
                                                                            ccAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCg                                                   >  2:74363/1‑52 (MQ=255)
                                                                             cAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAggc            <  1:204186/90‑1 (MQ=255)
                                                                               gTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAggcgg          <  1:691110/90‑1 (MQ=255)
                                                                                       ggCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTgggg  >  2:450603/1‑90 (MQ=255)
                                                                                       |                                                                                         
GTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGG  >  NZ_CP009273/296335‑296511

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