Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273832,0480TG56.2% 13.4 / 20.2 16V81G (GTC→GGC) hcxBhydroxycarboxylate dehydrogenase HcXB
Reads supporting (aligned to +/- strand):  ref base T (0/7);  new base G (9/0);  total (9/7)
Fisher's exact test for biased strand distribution p-value = 8.74e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACG  >  NZ_CP009273/831962‑832131
                                                                                      |                                                                                   
gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc                                                                                  >  1:493235/1‑90 (MQ=255)
gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc                                                                                  >  1:89880/1‑90 (MQ=255)
gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc                                                                                  >  2:152436/1‑90 (MQ=255)
           aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCGATGGCg                                                                       >  1:122917/1‑90 (MQ=255)
                  cTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCgcgc                                                                <  1:360033/90‑1 (MQ=255)
                     cTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACCCTCGATGGCGATCGCGCAtt                                                             <  2:263699/90‑1 (MQ=255)
                            cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCCCGCCCGATGGCGGCCGCGCCATTggggag                                                      >  2:234722/1‑86 (MQ=255)
                                                   tcatGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCACGCCCGATGGCGGTCGCGCATTTTGTTAGGGCGCGGCCCATGAAGCGATggg                               >  1:86360/1‑89 (MQ=255)
                                                     atGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGgcgc                             <  1:507187/90‑1 (MQ=255)
                                                                gTTAAAGAGGCGGGGGGGGCGGGCACGCCCGATGGCGATCGCGCATTTTGTCAGGGCGCGGCACATGAAGCGATGGGGCTGGGGGTTgag                  >  2:159888/1‑90 (MQ=255)
                                                                gTTAAAGAGGCGGGGGGGGCGGGCACGCCCGATGGCGATCGCGCATTTTGGCAGGGCGCGGCCCATGAAGCGATGGCGCTGGGGGTTgag                  >  2:210719/1‑90 (MQ=255)
                                                                gTTAAAGAGGCGGGGGCGGCGGGCACTCGCTGTGGCGGACACGGGTTTTGTCTGGGCGCGGGGCCTTAAGCGCCGGGGCTGGGGGTTgag                  >  2:191549/1‑90 (MQ=255)
                                                                     agagGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAgc             <  1:328371/90‑1 (MQ=255)
                                                                          cGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATc        <  1:194884/90‑1 (MQ=255)
                                                                             gggCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGc     <  2:86360/90‑1 (MQ=255)
                                                                                cggcggTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACg  <  2:494554/90‑1 (MQ=255)
                                                                                      |                                                                                   
GATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACG  >  NZ_CP009273/831962‑832131

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: