Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273700,0410TG55.6% 14.7 / 19.4 19G214G (GGT→GGGnagEPTS N‑acetyl glucosamine transporter subunit IIABC
Reads supporting (aligned to +/- strand):  ref base T (1/7);  new base G (10/0);  total (12/7)
Fisher's exact test for biased strand distribution p-value = 2.51e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-03
Rejected as consensus: Frequency below/above cutoff threshold.

ATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGATGTTCATGTCCGGCTTCTTCCCGATCATGATGTT  >  NZ_CP009273/699952‑700115
                                                                                          |                                                                          
aTCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGt                                                                            <  2:460439/90‑1 (MQ=255)
      aCCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGACGGGTTTTCCCGGGGACAtt                                                                      >  1:308988/1‑90 (MQ=255)
         ggTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGACGGTTTTTCACGGTGGCATTAAc                                                                   >  1:6691/1‑90 (MQ=255)
          gTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGTGACATTAAcc                                                                  <  1:261542/90‑1 (MQ=255)
               cATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTc                                                             <  1:406744/90‑1 (MQ=255)
                        cTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGCACGGGTTTCCCCGGGGACATTAACCGCTTCTATTCCggg                                                    >  1:414973/1‑89 (MQ=38)
                            acacCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGCTTTTCAACGGGGGATTAAACGGCTCTATGCCGGTGGCg                                                >  1:251440/1‑90 (MQ=255)
                              acCATCGCCTGTTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGc                                              <  2:152217/90‑1 (MQ=255)
                              acCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGACGGGTTTTCCCGGGGACATTAACCGCCTCTATGCCGGGggcggc                                              >  1:106170/1‑90 (MQ=255)
                                    gCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGGTTTTCACGGGGGGATTAAACGCTTCTATTCCCGTGACGGCGCcccg                                        >  1:369537/1‑90 (MQ=255)
                                    gCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTAAGGGTTTCCACGGGGGCCTTAACCGCTTCTATACCGGAGGCGGCACcgcg                                        >  2:372618/1‑90 (MQ=25)
                                         gTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGAt                                   <  1:150917/90‑1 (MQ=255)
                                         gTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGGTTTCCCCGGGGACAATAAACGCCTCTATTCCCGTGGCGGCACCCCcggggg                                   >  1:552837/1‑88 (MQ=255)
                                              aGATTGGTGAATTCACCAACGCGGCGGGTACGGGTTTTCCCGGGGACATTAACCGCTTCTATGCCGGGGGCGGCACCGCGGGGGTggtta                              >  1:69543/1‑86 (MQ=255)
                                                   ggTGAATTCACCAACGCGGCGGGGACGGGTTTTCACGGGGGCATTAACCGCTTCTATTCCGGGGACGGCACCGCGGGGGTGGTTATTTcc                         >  2:63460/1‑90 (MQ=255)
                                                   ggTGAATTCACCAACGCGGCGGGGACGGATTTCCACAAGGACAGTAAACACCTCCATTCCCCTGGCGGCACCCCGGGGGAGTTCCTGttc                         >  1:490191/1‑88 (MQ=255)
                                                              cAACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGATGTTCATGTCCGGCTTCTTccc              <  1:417264/90‑1 (MQ=255)
                                                                     gcggGTACGGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGATGTTCATGTCCGGCTTCTTCCCGATCatg       <  1:63460/90‑1 (MQ=255)
                                                                         ggaCGGGTTTCCACGGGGACATTAACCGCTTCTATGCCGGGGACGGCCCCGCGGGGGGGGTCATGTCCGGCTTCTTCCCGATCATGgtgt   >  1:74866/3‑90 (MQ=255)
                                                                          tACGGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGATGTTCATGTCCGGCTTCTTCCCGATCATGATGtt  <  2:308988/90‑1 (MQ=255)
                                                                                          |                                                                          
ATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGATGTTCATGTCCGGCTTCTTCCCGATCATGATGTT  >  NZ_CP009273/699952‑700115

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