Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273371,0340CG60.0% 6.4 / 10.7 10G40G (GGC→GGGmhpT3‑(3‑hydroxy‑phenyl)propionate transporter
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base G (6/0);  total (7/3)
Fisher's exact test for biased strand distribution p-value = 3.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GACCATCGGGCTTTGTTTTTTGGTCGCTCTGATGGAAGGGCTGGATCTTCAGGCGGCTGGCATTGCGGCGGGTGGCATCGCCCAGGCTTTCGCACTCGATAAAATGCAAATGGGCTGGATATTTAGCGCCGGAATACTCGGTTT  >  NZ_CP009273/370959‑371102
                                                                           |                                                                    
gACCATCGGGCTTTGTTTTTTGGTCGCTCTGATGGAAGGGCTGGATCTTCAGGCGGCTGGCATTGCGGCGGGTGGCa                                                                     <  2:87703/77‑1 (MQ=255)
gACCATCGGGCTTTGTTTTTTGGTCGCTCTGATGGAAGGGCTGGATCTTCAGGCGGCTGGCATTGCGGCGGGGGGCa                                                                     >  1:87703/1‑77 (MQ=255)
     tCGGGCTTTGTTTTTTGGTCGCTCTGATGGAAGGGCTGGATCTTCAGGCGGCTGGCATTGCGGCGGGGGGGGTCGCCCCGGGGTTTGCAc                                                   >  2:47549/1‑90 (MQ=255)
            ttGTTTTTTGGTCGCTCTGATGGAAGGGCTGGATCTTCAGGCGGCTGGCATTGCGGCGGGGGGGATCGCCCCGGCTTTTGCACTCGATaa                                            >  1:359439/1‑90 (MQ=255)
                          ctctGATGGAAGGGCTGGATCTTCAGGCGGCTGGCATTGCGGCGGGTGGCATCGCCCAGGCTTTCGCACTCGATAAAATGCAAATGGGCt                              <  1:302932/90‑1 (MQ=255)
                                    aGGGCTGGATCTTCAGGCGGCTGGCATTGCGGCGGGGGGGGTCGCCCCCGCTTTTCCCCCCCCTAAAAAACAAAAGGGGTGGGTAttttt                    >  2:302932/1‑88 (MQ=255)
                                       gCTGGATCTTCAGGCGGCTGGCATTGCGGCGGGTGGCATCGCCCAGGCTTTCGCACTCGATAAAATGCAAATGGGCTGGATATTTAgcgc                 <  1:212205/90‑1 (MQ=255)
                                                tCAGGCGGCTGGCATTGCGGCGGGGGGGATCGCCCAGGGTTTCGCACTCCGTAAAAAGCAAATAGGGGGGGTATTTAGCGGCCGAATaat        >  1:396727/1‑88 (MQ=255)
                                                    gcggcTGGCATTGCGGCGGGGGGGATCGCCCCGGGTTTTGCACTCGGTAAAAAGCAAAAGGGGTGGGTATTTTGCGCCGGGATACTCggg    >  2:54760/1‑89 (MQ=255)
                                                      ggcTGGCATTGCGGCGGGGGGGATCGCCCCGGCTTTTGCCCTCGATAAAAAGCAAAAGGGGTGGGTATTTTGCGCCGGGATACTCGGgtt  >  2:113364/1‑90 (MQ=255)
                                                                           |                                                                    
GACCATCGGGCTTTGTTTTTTGGTCGCTCTGATGGAAGGGCTGGATCTTCAGGCGGCTGGCATTGCGGCGGGTGGCATCGCCCAGGCTTTCGCACTCGATAAAATGCAAATGGGCTGGATATTTAGCGCCGGAATACTCGGTTT  >  NZ_CP009273/370959‑371102

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