Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,350,3810GC57.1% 2.0 / 16.1 14G457G (GGC→GGGcadClysine decarboxylation/transport transcriptional activator CadC
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base C (7/1);  total (10/4)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGGTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTGCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGTTTTACCTTTTACCA  >  NZ_CP009273/4350296‑4350469
                                                                                     |                                                                                        
gggTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCaa                                                                                      >  2:68826/1‑90 (MQ=255)
   tgtTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCa                                                                                   >  2:78407/1‑90 (MQ=255)
    gtTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCaa                                                                                  >  1:65979/1‑90 (MQ=255)
                     aaTTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTGCCAAGCAACACATAATTTAGCCAgg                                                                 <  2:152115/90‑1 (MQ=255)
                               ggTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTGCCAAGCAACACATAATTTAGCCAGGACATTTCAAg                                                       >  1:200188/1‑90 (MQ=255)
                                gTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTGCCAAGCAAcaca                                                                              <  1:45146/66‑1 (MQ=255)
                                gTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTGCCAAGCAAcaca                                                                              >  2:45146/1‑66 (MQ=255)
                                        tGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTGCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGcc                                              <  1:123355/90‑1 (MQ=255)
                                          cATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAg                                            >  1:245626/1‑90 (MQ=255)
                                            tCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTa                                          >  2:322639/1‑90 (MQ=255)
                                               gctgctTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTAttt                                       >  2:98367/1‑90 (MQ=255)
                                                       ccGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCaa                                                                                  >  1:243663/1‑39 (MQ=39)
                                                       ccGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCaa                                                                                  <  2:243663/39‑1 (MQ=39)
                                                                                    tGCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGGATTTATCGCCTGGTAAGACTCATCTTTTTTACCTTTTACCa  >  2:128723/1‑90 (MQ=255)
                                                                                     |                                                                                        
GGGTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTGCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGTTTTACCTTTTACCA  >  NZ_CP009273/4350296‑4350469

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