Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273832,0420CG57.9% 5.1 / 18.8 20A79G (GCG→GGG) hcxBhydroxycarboxylate dehydrogenase HcXB
Reads supporting (aligned to +/- strand):  ref base C (2/6);  new base G (11/0);  total (14/6)
Fisher's exact test for biased strand distribution p-value = 1.03e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.17e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAG  >  NZ_CP009273/831962‑832115
                                                                                |                                                                         
gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGggggggggggggac                                                                  >  1:18374/1‑86 (MQ=255)
gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGgggggggggggccc                                                                  >  1:762055/1‑86 (MQ=255)
gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGgggcggggac                                                                  >  2:50917/1‑86 (MQ=255)
gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGcggcggggcc                                                                  >  1:167911/1‑86 (MQ=255)
gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc                                                                  >  2:793014/1‑90 (MQ=255)
gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc                                                                  >  2:329436/1‑90 (MQ=255)
gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCAc                                                                  >  1:19863/1‑90 (MQ=255)
gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCAc                                                                  <  1:793013/90‑1 (MQ=255)
gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCAc                                                                  >  2:656976/1‑90 (MQ=255)
           aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACCCTCGATGGCg                                                       >  2:812932/1‑90 (MQ=255)
           aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCTCGATgggg                                                       >  2:166046/1‑88 (MQ=255)
                cGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATcgc                                                  <  1:168571/90‑1 (MQ=255)
                  cTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCCCCCTGGCGGTCgcgc                                                >  1:680537/1‑90 (MQ=255)
                            cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGTGTGGGGGTTACTCTCGAGGGGGATTTCGCATTTTGTctg                                      >  2:541647/1‑88 (MQ=255)
                            cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCTCGATGGCGATCGCGCATTTGGTCAg                                      >  2:45801/1‑90 (MQ=255)
                            cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCt                                                               <  1:64750/65‑1 (MQ=255)
                            cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCt                                                               >  2:64750/1‑65 (MQ=255)
                                    ccTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGc                              <  1:166046/90‑1 (MQ=255)
                                      tGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCTCGGTGGGGGTCGCGCATTTGGGCAGGGCGCggggc                            >  2:836908/1‑87 (MQ=255)
                                        cAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACAt                          <  1:234420/90‑1 (MQ=255)
                                              aaCCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCCCTGGCGGTCGCGCATTTTGGTAGGGCGCGGGCCATGAAGCg                    >  2:457204/1‑90 (MQ=255)
                                                      tGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCCCGCCCCCTGGGGGGCGCGCCTTTTGTCAGGGCGGGGCACAAGAAGCGATGGCGCt            >  2:445100/1‑90 (MQ=255)
                                                       gCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTg           <  1:591212/90‑1 (MQ=255)
                                                                gTTAAAGAGGCGGGGGGGGCGGGCACGCTCGATGGGGATCGCGCATTTTGTCAGGGCGCGGGACATGAAGCGGTGGGGCTGGGGGTTgag  >  1:181363/1‑90 (MQ=255)
                                                                                |                                                                         
GATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAG  >  NZ_CP009273/831962‑832115

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