Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273296,4020TG56.2% 6.6 / 13.1 16T82P (ACC→CCC) paoCaldehyde oxidoreductase molybdenum‑binding subunit PaoC
Reads supporting (aligned to +/- strand):  ref base T (3/4);  new base G (9/0);  total (12/4)
Fisher's exact test for biased strand distribution p-value = 1.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATAT  >  NZ_CP009273/296313‑296491
                                                                                         |                                                                                         
ggCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGt                                                                                           <  1:256184/90‑1 (MQ=255)
   cTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGttag                                                                                        >  2:314676/1‑87 (MQ=255)
    tgatgaTAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATg                                                                                       <  1:41907/90‑1 (MQ=255)
     gatgatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAtta                                                                                      >  1:326995/1‑88 (MQ=255)
               cTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGc                                                                                  <  2:322804/84‑1 (MQ=255)
               cTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGc                                                                                  >  1:322804/1‑84 (MQ=255)
                        ggggCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGAAAGCCAGTACGCCCCgcgc                                                                   >  2:395055/1‑90 (MQ=255)
                        ggggCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAAGACACGCCGTACGCCCCgcgc                                                                   >  1:186097/1‑90 (MQ=255)
                           gccgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAACCCGGACCCCCCGCGCctt                                                                >  2:322135/1‑90 (MQ=255)
                           gccgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGAAAACCCGGACGCCCCGCGCctt                                                                >  1:407450/1‑90 (MQ=255)
                           gccgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAAGAAAGCCCGTACGCCCCGCGCctt                                                                >  1:158317/1‑90 (MQ=255)
                             cgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCAGTACCCCCCGCGCttttt                                                              >  1:91681/1‑90 (MQ=255)
                             cgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAATACCACCCCGACGCCCCGCGCctttt                                                              >  1:197768/1‑90 (MQ=255)
                              gccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGt                                                                                           <  2:329405/60‑1 (MQ=255)
                                      gCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGGAAGCCCCGCGCTTTTTTCgcggcgg                                                     >  2:143310/1‑89 (MQ=255)
                                         tGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGGATGACAGCCAGGACCCCCCGCGCGTTTTTCGCCGCGTCCg                                                  >  1:212444/1‑90 (MQ=255)
                                                gtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAg                                           <  2:427304/90‑1 (MQ=255)
                                                        cGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAAAGAAAGCCCCTACGCCCGGGGCTTTTTTGGCGGCGGCCGGATCAAggggggggg                                   >  2:16441/1‑87 (MQ=255)
                                                           cTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGc                                <  1:322135/90‑1 (MQ=255)
                                                                                         tAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGatat  <  1:249946/90‑1 (MQ=255)
                                                                                         tAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGatat  <  1:46549/90‑1 (MQ=255)
                                                                                         |                                                                                         
GGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATAT  >  NZ_CP009273/296313‑296491

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