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

GGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATAT  >  NZ_CP009273/296313‑296491
                                                                                         |                                                                                         
ggCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGt                                                                                           <  2:21979/90‑1 (MQ=255)
        gatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAg                                                                                   <  1:18024/90‑1 (MQ=255)
        gatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAg                                                                                   <  2:168274/90‑1 (MQ=255)
          tAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGcc                                                                                 >  2:165687/1‑90 (MQ=255)
                      gTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCAGGACCCCCCgc                                                                     >  2:235315/1‑90 (MQ=255)
                      gTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAATACAGCCCCTACCCCCCgc                                                                     >  1:383634/1‑90 (MQ=255)
                        ggggCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGGAGGACAACCCGGACCCCCCgcgc                                                                   >  1:454179/1‑90 (MQ=255)
                        ggggCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCAGTAAGCCCcgcgg                                                                   >  2:69488/1‑89 (MQ=255)
                        ggggCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAACCAGGACGCCCCgcgc                                                                   >  2:398690/1‑90 (MQ=255)
                           gccgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAAGACAGCCAGGAAGCCCCGCGCctt                                                                >  2:241190/1‑90 (MQ=255)
                                                gtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAg                                           <  1:83189/90‑1 (MQ=255)
                                                                    gTGCCCCGGCGTCACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGGGAGGCGGCCTTTgg                       >  2:192940/1‑90 (MQ=255)
                                                                      gCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCa                     <  2:199228/90‑1 (MQ=255)
                                                                             cGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCgg              <  1:167661/90‑1 (MQ=255)
                                                                                         tAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGatat  <  2:441342/90‑1 (MQ=255)
                                                                                         |                                                                                         
GGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATAT  >  NZ_CP009273/296313‑296491

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