Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273294,6640AC56.2% 12.8 / 19.1 16V661G (GTG→GGG) paoCaldehyde oxidoreductase molybdenum‑binding subunit PaoC
Reads supporting (aligned to +/- strand):  ref base A (7/0);  new base C (0/9);  total (7/9)
Fisher's exact test for biased strand distribution p-value = 8.74e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-04
Rejected as consensus: Frequency below/above cutoff threshold.

ACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCG  >  NZ_CP009273/294577‑294739
                                                                                       |                                                                           
acacCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCAcc                                                                           >  2:294747/1‑90 (MQ=255)
             tCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCa                                                              >  1:275657/1‑90 (MQ=255)
                  ggggAGGATGGGGGGTCGTTATCACCCAGAAAAACCCCCTCTTTTTTTGGGATATCGGGATGACCCGCCCCCTCATCCCCCGCCATATCg                                                         <  1:76736/90‑1 (MQ=255)
                   gggAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCa                                                                             <  2:523127/69‑1 (MQ=255)
                   gggAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCa                                                                             >  1:523127/1‑69 (MQ=255)
                        gaTATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTa                                                   >  1:160430/1‑90 (MQ=255)
                        gaTATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTa                                                   >  1:498522/1‑90 (MQ=255)
                          tggggggTGGGTATCACCAAGAAAAACCCCCTCCTTTTTTGGGATATCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAAc                                                 <  1:107984/87‑1 (MQ=255)
                              gggTCGGTACCCCCAAGAAAAACCCCCCCTTTTTTTGGGATATCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAg                                             <  1:100665/90‑1 (MQ=255)
                               ggTGGGTACCACCCAGAAAAACCCCCCCCTGTTTTGGGATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGt                                            <  2:82320/90‑1 (MQ=255)
                                  cGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGc                                         >  1:324723/1‑90 (MQ=255)
                                      atcaccCAGAAAAACCCCCCCCTTTTTGGGGATACCCGAATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCaa                                     <  2:501121/90‑1 (MQ=255)
                                          tcCAGAAAAACCCCCTCCTTTTTTGGGATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACgg                                 <  2:403167/90‑1 (MQ=255)
                                              ggAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCa                                                                             <  1:458951/42‑1 (MQ=255)
                                              ggAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCa                                                                             >  2:458951/1‑42 (MQ=255)
                                                       ccccctTTTTGGGGATACCCGAATGACCCGCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCa                    <  2:512485/87‑1 (MQ=255)
                                                         tcctGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGc                  >  1:373820/1‑90 (MQ=255)
                                                                  gggATACCCGCATGACCCGGCCCCCCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc         <  1:103007/90‑1 (MQ=255)
                                                                      tatCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTg     >  2:527662/1‑90 (MQ=255)
                                                                         ccacATGACCCGCCCCCCCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcg  <  1:16512/87‑1 (MQ=255)
                                                                                       |                                                                           
ACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCG  >  NZ_CP009273/294577‑294739

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