Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,273,8300AC58.8% 12.3 / 19.0 17T85T (ACA→ACCagaVPTS N‑acetylgalactosamine transporter subunit IIB
Reads supporting (aligned to +/- strand):  ref base A (7/0);  new base C (0/10);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 5.14e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.65e-04
Rejected as consensus: Frequency below/above cutoff threshold.

ATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAA  >  NZ_CP009273/3273741‑3273913
                                                                                         |                                                                                   
aTCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAaca                                                                                     >  1:558455/1‑90 (MQ=255)
                 gACGCTAAAAAAATTTATCAAAAACATTCCCCCCGCCGCCGATCGAAAAAAAATCCTGCTGGTTTGAAAAACCCCCGCCGATTTCCTGAc                                                                    <  2:524152/90‑1 (MQ=255)
                   cGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGc                                                                  >  1:412760/1‑90 (MQ=255)
                     cTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTg                                                                >  1:541186/1‑90 (MQ=255)
                                     acaactTTCCTCGCGCCGCCAATCGACAGAAAATCCTGCTGGTTTGTAAACCCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTc                                                <  1:593336/90‑1 (MQ=255)
                                        acaTTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCgg                                             >  1:182779/1‑90 (MQ=255)
                                        acaTTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCgg                                             >  2:406015/1‑90 (MQ=255)
                                             ccccgcgCCGCCGATCGAAAAAAAATCCTGCTGGTTTGAAAACCCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAAt                                        <  1:286688/87‑1 (MQ=255)
                                               tCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCg                                      >  2:603708/1‑90 (MQ=255)
                                                  ccccgccgACCGAAAAAAAATCCTGTTGTTTTTAAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCAt                                   <  2:179138/88‑1 (MQ=255)
                                                   gccgccGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCAtt                                  >  2:174014/1‑90 (MQ=255)
                                                           ccGACAAAAAACCCTGCTGTTTTGTAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTgg                          <  1:557881/89‑1 (MQ=255)
                                                             gACAGAAAACCCTGCTGGTTTGAAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCa                        <  1:227697/90‑1 (MQ=255)
                                                                 gAAAATCCTGCTGGTTTGTAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAatat                    <  1:115587/90‑1 (MQ=255)
                                                                 aaaaaCCCTGCTGTTTTGAAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAatat                    <  1:573997/89‑1 (MQ=255)
                                                                     aTCCTGTTGTTTTGAAAACCCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCAc                <  1:453480/90‑1 (MQ=255)
                                                                        ctgctgGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTAc             >  2:363181/1‑90 (MQ=255)
                                                                                   aaaaaCCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCaa  <  2:541186/89‑1 (MQ=255)
                                                                                         |                                                                                   
ATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAA  >  NZ_CP009273/3273741‑3273913

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