Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273832,0420CG56.5% 10.3 / 26.9 23A79G (GCG→GGG) hcxBhydroxycarboxylate dehydrogenase HcXB
Reads supporting (aligned to +/- strand):  ref base C (2/8);  new base G (13/0);  total (15/8)
Fisher's exact test for biased strand distribution p-value = 9.18e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCAC  >  NZ_CP009273/831962‑832130
                                                                                |                                                                                        
gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGccgcggggac                                                                                 >  2:437941/1‑86 (MQ=255)
gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc                                                                                 >  1:312686/1‑90 (MQ=255)
gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc                                                                                 >  2:585062/1‑90 (MQ=255)
gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc                                                                                 >  1:61341/1‑90 (MQ=255)
      aaGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGTCACGCTCGa                                                                           >  1:120766/1‑90 (MQ=255)
           aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGGACGCCCCATgggg                                                                      >  2:527456/1‑88 (MQ=255)
                     cTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACCCTCGACTGCGGACGCGCact                                                            >  2:368149/1‑88 (MQ=255)
                            cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGTCACg                                                                                >  1:593826/1‑63 (MQ=255)
                            cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGTCACGCTCGATGGCGGTCGCGCATTTGGTCAg                                                     >  1:19376/1‑90 (MQ=255)
                            cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCTCGATGGCGATCGCGCATTTGGTCAg                                                     >  1:12107/1‑90 (MQ=255)
                            cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCCATGGGGATCGCGCCTTTTGGCAg                                                     >  2:90183/1‑90 (MQ=255)
                            cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACg                                                                                <  2:593826/63‑1 (MQ=255)
                            cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCACGCCCCCTGGGGGTCGCGCCTTTTGGCAg                                                     >  2:19818/1‑90 (MQ=255)
                                             tAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCCCCATGGCGGACGCGCCTTTTGGCAGGGCGCGGGACAAGaaac                                    >  2:368393/1‑88 (MQ=255)
                                              aaCCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCACCCTCCATGGCGATCGCGCATTTGGGCAGGGCGCCGGACATGAAGCg                                   >  1:134858/1‑90 (MQ=255)
                                                      tGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCt                           <  2:61341/90‑1 (MQ=255)
                                                                gTTAAAGAGGCGGGGGGGGCGGGCACGCTCGATGGCGATCGCGCCTTTGGTCAGGGGGCGGCCCATGAAGCGATGGCGCTGGGGGTTgag                 >  2:90897/1‑90 (MQ=255)
                                                                     agagGCGGGGGGGGCGGGCACGCCCGATGGCGATCGCGCATTTTGGCAGGGCGCGGCACATGAAACGATGGGGCTGGGGGTTGAGAAAgc            >  2:202826/1‑90 (MQ=255)
                                                                     agagGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAgc            <  2:134858/90‑1 (MQ=255)
                                                                          cGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATc       <  2:228121/90‑1 (MQ=255)
                                                                          cGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATc       <  2:120766/90‑1 (MQ=255)
                                                                          cGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATc       <  1:585062/90‑1 (MQ=255)
                                                                          cGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATc       <  1:527456/90‑1 (MQ=255)
                                                                             gggCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGc    <  1:368149/90‑1 (MQ=255)
                                                                               gcggcgGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCAc  <  1:605378/90‑1 (MQ=255)
                                                                                |                                                                                        
GATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCAC  >  NZ_CP009273/831962‑832130

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