Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273296,4220GC55.6% 8.0 / 24.8 18P75R (CCG→CGG) paoCaldehyde oxidoreductase molybdenum‑binding subunit PaoC
Reads supporting (aligned to +/- strand):  ref base G (1/7);  new base C (10/0);  total (11/7)
Fisher's exact test for biased strand distribution p-value = 2.51e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGG  >  NZ_CP009273/296337‑296503
                                                                                     |                                                                                 
ggggCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGGATGACAGCCCGCACGCCCCgcgc                                                                               >  2:377838/1‑90 (MQ=255)
     cgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAAAGACAACCCGTACGCCCCGCGCttttt                                                                          >  1:58240/1‑90 (MQ=255)
     cgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGGCAGTACGCCCCGGGCttttt                                                                          >  1:97818/1‑90 (MQ=255)
       ccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAAAACAGCCAATACGCCCCGCCCctttttg                                                                        >  1:301771/1‑88 (MQ=255)
         tAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTgcgc                                                                      <  1:48708/90‑1 (MQ=255)
                 tGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACCCCCCGCGCTTTTTTCGCGGCGTCCg                                                              >  1:141445/1‑90 (MQ=255)
                        gtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTTCGCGGCGTCCGTATCAAg                                                       >  1:80971/1‑90 (MQ=255)
                                    tttGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCg                                           <  1:287825/90‑1 (MQ=255)
                                            gTGCCCCGGCGTTACTGGCGGTAATTACAACCCCTACGCCCCGCGCCTTTTGCGCGGCGTCCGGATCAAGGGGGGGGAGGCGGCCttttg                                   >  1:149055/1‑88 (MQ=255)
                                            gTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTGCGCGGCGTCCGTATCAAGGGGGGGGGGGCGGCCTTTgg                                   >  1:387523/1‑90 (MQ=255)
                                               ccccGGCGTTACTGGCGGTAATGACAACCAGGACGCCCCGCGCTTTTTTCGCGGCGGCCGGATCAAGGGCGGGGAGGCGGCCTTTGGCaa                                >  2:327495/1‑90 (MQ=255)
                                                      gTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCCTTTTTCGCGGCGGCCGGATCAAGGGGGGGGGGGGGGCCCTTGGCCATGGCggg                         >  2:368493/1‑89 (MQ=255)
                                                          cTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcg                     <  2:97818/90‑1 (MQ=255)
                                                            ggCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcgac                   <  1:378118/90‑1 (MQ=255)
                                                            ggCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcgac                   <  1:259204/90‑1 (MQ=255)
                                                             ggggTAAGGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcgacg                  <  2:220317/88‑1 (MQ=255)
                                                                 tAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGatat              <  2:387866/90‑1 (MQ=255)
                                                                             gTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAggcgg  >  2:93338/1‑90 (MQ=255)
                                                                                     |                                                                                 
GGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGG  >  NZ_CP009273/296337‑296503

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