Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,152,7450AG54.5% 5.6 / 11.0 11G171G (GGA→GGG) ycfHmetal‑dependent hydrolase
Reads supporting (aligned to +/- strand):  ref base A (1/4);  new base G (6/0);  total (7/4)
Fisher's exact test for biased strand distribution p-value = 1.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGT  >  NZ_CP009273/1152687‑1152830
                                                          |                                                                                     
gCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCt                                                        <  2:208363/90‑1 (MQ=255)
        cacTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAAt                                                <  1:218374/90‑1 (MQ=255)
         acTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGGTCTCGGGTTTTAAATCTCCTTTTCCGGCATTGTGACCTTCCGTAATg                                               >  1:58780/1‑90 (MQ=255)
            gTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTGGGTCTCGGGTTTTAAATCTCCTTTTTCCGCCTTTTGGACTTCCCTTATGCgg                                            >  2:110363/1‑90 (MQ=255)
               ttACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGc                                         <  2:113519/90‑1 (MQ=255)
                      ggACAGAGAAACGGCGGGGAAAATATTGGGTTTCGGGTTTTTTATTTTCTTTTTCGGCCTTTTGGCCTTCCCTAATGCGGGGCAACTgcg                                  >  2:337857/1‑90 (MQ=255)
                      ggACAGAGAAACGGCGGGGAAAATACTGGGTCTCGGGTTTTTTATCTCCCTTTTCCGCATTGTGGCCTTCCCTAAAGCGGGGGAACTgcg                                  >  2:176758/1‑90 (MQ=255)
                      ggACAGAGAAACGGCGGGGAAAATACTGGGTCTCGGGTTTTAAATCTCCTTTTTCGGGATTGTGGCCTTCCCTAATGCGGGGGAAATgcg                                  >  1:22480/1‑90 (MQ=255)
                                aCGGCGGGGAAAATAATGGATCTCGGGTTTTTAATCTCCTTTTTCGGCATTGTGGGCTTCCCTAATGCGGGGGAACTGCGCGATGCTgcg                        >  2:48994/1‑90 (MQ=255)
                                                  gATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGAt      <  2:17506/90‑1 (MQ=255)
                                                      tcGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGt  >  2:320715/1‑90 (MQ=255)
                                                          |                                                                                     
GCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGT  >  NZ_CP009273/1152687‑1152830

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