Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,562,1230AT46.4% 29.5 / 12.3 28I189N (ATT→AAT) glgCglucose‑1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (10/5);  new base T (0/13);  total (10/18)
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.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGTAGATACCCATACTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTCAACGAATTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAA  >  NZ_CP009273/3562040‑3562203
                                                                                     |                                                                                
cGTAGATACCCATACTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGGTTTTCAACGAATTCGATAATTTTAt                                                                            >  1:162208/1‑90 (MQ=255)
   aGATACCCATACTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTCAACGAATTCGATAATTTTATCGt                                                                         <  2:590499/90‑1 (MQ=255)
     aTACCCATACTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTCAACGAATTCGATAATTTTATCGTtc                                                                       <  1:204875/90‑1 (MQ=255)
     aTACCCATACTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTCAACGAATTCGATAATTTTATCGTtc                                                                       <  1:623849/90‑1 (MQ=255)
          cATACTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTCAACGAATTCGATAATTTTATCGTTCtcatc                                                                  >  1:407209/1‑90 (MQ=255)
                          ttGCTGGGTTCGTTCGGCATTGGGGGGGGGTTAGCAGTTTTTTCAACGAATTCGATATTTTTATCTTTCTCATCAACCGCCATAACGCCa                                                  <  1:103575/90‑1 (MQ=255)
                          ttGCTCGGTTCGTTCGGCTTGGAGGGGGGGTTAGCAGTTTTTTCAACAAATTCGATATTTTTATCTTTCTCATCAACCGCCATAACGCCa                                                  <  1:69832/90‑1 (MQ=255)
                          ttGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTCAACGAATTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCa                                                  >  2:552035/1‑90 (MQ=255)
                            gCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTCAACGAATTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCaaa                                                >  2:519188/1‑90 (MQ=255)
                                 gATCGTTCGGCATTGAGGGGGGGTTAGCGGTTTTTAAAACAAATTCGATATTTTTATCCTTCTCATCAACCGCCATAACGCCAAATGCgg                                           <  2:479149/90‑1 (MQ=255)
                                 gATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTCAACGAATTCGATAATTTTATCGt                                                                         <  1:213203/60‑1 (MQ=255)
                                 gATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTCAACGAATTCGATAATTTTATCGt                                                                         >  2:213203/1‑60 (MQ=255)
                                 gATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTCAACGAATTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCgg                                           >  2:164175/1‑90 (MQ=255)
                                               ggcgggggTTTAGCAGTTTTTTAAACGAATTCGATATTTTTATCCTCCTCATCACCCGCCATAACGCCAAATGCGGAGGCTTCTTCAATc                             <  2:67871/88‑1 (MQ=255)
                                               ggcgggggGTTAGCAGTTTTTACAACGAATTCGATATTTTTATCTTTCTCACCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATc                             <  1:376135/88‑1 (MQ=255)
                                               gACGGCGGGTTAGCAGGTTTTTCAACGAATTCGATAATTTTATCGt                                                                         <  2:169006/46‑1 (MQ=255)
                                               gACGGCGGGTTAGCAGGTTTTTCAACGAATTCGATAATTTTATCGt                                                                         >  1:169006/1‑46 (MQ=255)
                                                gcgggggGTTAAGAGGTTTTTTCAACAAATTCAATATTTTTATCTTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATc                             <  2:539386/89‑1 (MQ=255)
                                                  ggggggTTAGCAGTTTTTTCAACGAATTCGATATTTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGt                          <  2:259269/87‑1 (MQ=255)
                                                   gcggGTTAGCAGGTTTTTCAACGAATTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTa                         >  2:353135/1‑90 (MQ=255)
                                                    ggggTTAGCAGTTTTTTAAACAAATTCGATATTTTTATCTTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTAc                        <  2:197416/89‑1 (MQ=255)
                                                    ggggTAAGCGGTTTTTTAAACAAATTCGATATTTTTATCTTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTAc                        <  1:98206/89‑1 (MQ=255)
                                                      ggTTAGCAGGTTTTTCAACGAATTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTg                      >  2:465099/1‑90 (MQ=255)
                                                               gttttttAAACAAATTCGATATTTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACaa              <  1:100852/88‑1 (MQ=255)
                                                                 tttttAAACGAATTGGATATTTTTATCGTTCTCATAAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGc            <  2:90178/90‑1 (MQ=255)
                                                                  ttttCAAACAAATTCAATATTTTTATCTTTCCCACCACCCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGc            <  1:149100/90‑1 (MQ=255)
                                                                  ttttCAACGAATTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCa           >  1:171327/1‑90 (MQ=255)
                                                                            aaTTCAATATTTTTATCGTTCTCATCACCCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCaa  <  2:407209/90‑1 (MQ=255)
                                                                                     |                                                                                
CGTAGATACCCATACTCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGGTTTTTCAACGAATTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAA  >  NZ_CP009273/3562040‑3562203

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