Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,216,1540TG56.2% 8.0 / 14.5 16V164G (GTG→GGG) dldD‑lactate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (3/4);  new base G (9/0);  total (12/4)
Fisher's exact test for biased strand distribution p-value = 1.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.25e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCGAACCGCACTCAGTGATTGGATCATCGTGTATAGGCGCATCGGTCATCGGCGGTATTTGTAACAACTCCGGCGGCTCGCTGGTGCAACGTGGCCCGGCGTATACCGAAATGTCGTTATTCGCGCGTATAAATGAAGACGGCAAACTGACGCTGGTGAACCATCTG  >  NZ_CP009273/2216070‑2216236
                                                                                    |                                                                                  
gcgAACCGCACTCAGTGATTGGATCATCGTGTATAGGCGCATCGGTCATCGGCGGTATTTTTAACAACTCCCGCGGCTCGCTGGGGCAAc                                                                               >  2:90839/1‑90 (MQ=255)
   aaCCGCACTCAGTGATTGGATCATCGTGTATAGGCGCATCGGTCATCGGCGGTATTTTTAACAACTCCGGCGGCTCGCTGGGGCAACGTg                                                                            >  2:372700/1‑90 (MQ=255)
           caggTGATTGGATCATCGTGTATAGGCGCATCGGTCATCGGCGGTATTTGTAACAACTCCGGCGGCTCGCTGGTGCAACGTGGCCCGGCg                                                                    <  2:862501/87‑1 (MQ=255)
                        catcGTGTATAGGCGCATCGGTCATCGGCGGTATTTTTTAAAAATCCCGCGGGCCGCTGGGGGAACGGGGCCCGGCGGATACCGaaacag                                                       >  1:862501/1‑87 (MQ=255)
                          tcGTGTATAGGCGCATCGGTCATCGGCGGAATTTGTAAAAACTCCGGCGGGCCGCTGGGGCAACGGGGCCCCGCGGATACCGAAATGTCg                                                     >  1:309196/1‑90 (MQ=255)
                              gtATAGGCGCATCGGTCATCGGCGGTATTTGTAACAACTCCGGCGGCTCGCTGGGGCAACGTGGCCCCGCGGATACCGAAATGTCGtttt                                                 >  2:634426/1‑88 (MQ=255)
                                   ggCGCATCGGTCATCGGCGGTATTTGTAACAACTCCGGCGGCTCGCTGGTGCAACGTGGCCCGGCGTATACCGAAATGTCGTTATTcgcg                                            <  1:634426/90‑1 (MQ=255)
                                      gcATCGGTCATCGGCGGTATTTGTAACAACTCCGGCGGCTCGCTGGGGCAACGTGGCCCGGCGTATACCGAAATGTCGTTATTCGCGCGt                                         >  2:430529/1‑90 (MQ=255)
                                      gcATCGGTCATCGGCGGTATTTGTAACAACTCCGGCGGCTCGCTGGGGCAACGTGGCCCGGCGTATACCCGAAAGTCGTTTTTTGCGGGt                                         >  2:524843/1‑90 (MQ=255)
                                       cATCGGTCATCGGCGGTATTTGTAACAACTCCGGCGGCTCGCTGGGGCAACGAGGCCCGGCGTATACCGAAATGTCGTTATTCGCGCGta                                        >  2:311560/1‑90 (MQ=255)
                                        aTCGGTCATCGGCGGTATTTTTAACAACTCCGGCGGCTCGCTGGTGCAACGTGGCCCCGCGTATACCGAAATGTCGTTATTTGCGCGgat                                       >  2:750268/1‑90 (MQ=255)
                                        aTCGGTCATCGGCGGTATTTGTAACAACTCCGGCGGCTCGCTGGGGCAACGTGGCCCCGCGGATACCGAAATGTCGTTATTTGCGCGtat                                       >  2:750269/1‑90 (MQ=255)
                                        aTCGGTCATCGGCGGGATTTTTAAAAACTCCCGCGGCTCGCTGGTGCAAACGGGCCCCGCGGATACCGAAAAGTCGTTATTCGCGCGtat                                       >  1:377134/1‑90 (MQ=255)
                                                                caacTCCGGCGGCTCGCTGGTGCAACGTGGCCCGGCGTATACCGAAATGTCGTTATTCGCGCGTATAAATGAAGACGGCAAACTGACGCt               >  1:536809/1‑90 (MQ=255)
                                                                           gcTCGCTGGTGCAACGTGGCCCGGCGTATACCGAAATGTCGTTATTCGCGCGTATAAATGAAGACGGCAAACTGACGCTGGTGAACCATc    <  1:90839/90‑1 (MQ=255)
                                                                             tCGCTGGTGCAACGTGGCCCGGCGTATACCGAAATGTCGTTATTCGCGCGTATAAATGAAGACGGCAAACTGACGCTGGTGAACCATCTg  <  2:848257/90‑1 (MQ=255)
                                                                                    |                                                                                  
GCGAACCGCACTCAGTGATTGGATCATCGTGTATAGGCGCATCGGTCATCGGCGGTATTTGTAACAACTCCGGCGGCTCGCTGGTGCAACGTGGCCCGGCGTATACCGAAATGTCGTTATTCGCGCGTATAAATGAAGACGGCAAACTGACGCTGGTGAACCATCTG  >  NZ_CP009273/2216070‑2216236

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