Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,896,6400GC54.2% 18.7 / 35.8 24A141G (GCT→GGT) bglB6‑phospho‑beta‑glucosidase BglB
Reads supporting (aligned to +/- strand):  ref base G (10/1);  new base C (0/13);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 5.61e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-06
Rejected as consensus: Frequency below/above cutoff threshold.

CAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTTGATCCCCGCCTGCGCCAT  >  NZ_CP009273/3896558‑3896725
                                                                                  |                                                                                     
cAGCCATAACGCCATTTTTTGTTGGTAGCGAGTAAAGCCCGTGGGGGCGTAATGCTCGAAGTGGTCGATGCCCGCTCGATTACCCCAAcc                                                                                <  1:394338/90‑1 (MQ=255)
               ttttTGTTGGGGGCGAGAAAAGCCCGGGGGGGGCTAACGATCGAAGTAATCGATTCCCGCTAAATTACCCAACCCGCCGTAGTTTTTCAc                                                                 <  2:318649/90‑1 (MQ=255)
               ttttTGTGGGTAGGGAGAAAAGCCCGGGGGGGCGTAATGCTCGAAGGGGCCGATGCCCGCTCGATTACCCCACCCGCCGTAGTTTTTCAc                                                                 <  2:323065/90‑1 (MQ=255)
               tttATTTGGGTGGGGAGAAAAGCCCGGGGGGGCGTAATGCTGGAAGTGATCGATGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCAc                                                                 <  2:134900/90‑1 (MQ=255)
                    gTTGGTAGGGAAAAAAGCCCGGGGGGGCGTAATGCTGGAAGTGGTCGATCCCCGCTCGATTACCCCACCCGCCGTAGTTTTTCACCAGcc                                                            <  2:285115/90‑1 (MQ=255)
                    gTTGGTAGCGAGAAAAGCCCGGGGGGGGGTAATGCTGGAAGTGATCGATGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCACCAGcc                                                            <  2:269866/90‑1 (MQ=255)
                    gTTGGTAGCGAGAAAAGCCCGGGGGGGCGTAATGCTCGAGGGGGTCGATGCCCGCTCGATTACCCAAACCGCCGTAGTTTTTCACCAGcc                                                            <  2:437893/90‑1 (MQ=255)
                        gTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCata                                                        >  1:358675/1‑90 (MQ=255)
                                  aaGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCATCAGCCCATATGGCATTTCg                                              >  2:380162/1‑90 (MQ=255)
                                        gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAGTgg                                        >  1:89204/1‑90 (MQ=255)
                                        gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                                        >  2:226916/1‑90 (MQ=255)
                                        gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                                        >  2:424826/1‑90 (MQ=255)
                                          gggggCGTAATGCTGGAGGGGGTCGATGCCCGCTCGTTTACCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGa                                      <  2:149330/88‑1 (MQ=255)
                                          gCGGGCGTAATGCTCGAGGTGGTCGATGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGa                                      <  1:142547/90‑1 (MQ=255)
                                           ggggCATAATGCTCGAGGGGGTCGATCCCCCCTCGATTACCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGAt                                     <  2:444577/89‑1 (MQ=255)
                                                gTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGt                                >  1:289170/1‑90 (MQ=255)
                                                       tCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGc                         >  1:180930/1‑90 (MQ=255)
                                                        ggaggTGGTCGATGCCCGCTATTTTACCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCg                        <  2:259379/86‑1 (MQ=255)
                                                            ggggtCGATGCCCGCTCTATTACCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCtt                    <  2:423998/86‑1 (MQ=255)
                                                               gtCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTTGAt                 >  2:78853/2‑90 (MQ=255)
                                                                         gCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTTGATCCCCGCCTgc       >  1:62384/1‑90 (MQ=255)
                                                                              tttACCCCACCCGCCGTAGTTTTTCCCCAGCCCATATGGCATTTCGTAATGGGATAACGTTAACAGCGGCTTGATCCCCGCCTGCGCCAt  <  2:342886/89‑1 (MQ=255)
                                                                                 aGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCtt                    <  1:105519/69‑1 (MQ=255)
                                                                                 aGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCtt                    >  2:105519/1‑69 (MQ=255)
                                                                                  |                                                                                     
CAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTTGATCCCCGCCTGCGCCAT  >  NZ_CP009273/3896558‑3896725

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