Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,175,9570AT51.8% 13.9 / 32.3 27I17F (ATT→TTT) zupTzinc transporter ZupT
Reads supporting (aligned to +/- strand):  ref base A (5/8);  new base T (14/0);  total (19/8)
Fisher's exact test for biased strand distribution p-value = 5.80e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-03
Rejected as consensus: Frequency below/above cutoff threshold.

GATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTC  >  NZ_CP009273/3175872‑3176046
                                                                                     |                                                                                         
gATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTgg                                                                                       <  2:634173/90‑1 (MQ=255)
  tCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGAGCCACGGTTTTTTgcg                                                                                     >  2:369641/1‑90 (MQ=255)
               gAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTTTTTGCGCGTTTTTCGGCGt                                                                        >  2:846032/1‑90 (MQ=255)
               gAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGACGGCACGGTTTTTGGTGCTTTTCTCgtggg                                                                        >  2:904677/1‑86 (MQ=255)
                   aTTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACCTTTTTTGGGGCGGTTTTCCGCGGtctc                                                                    >  2:561173/1‑90 (MQ=255)
                            gAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGCCCCGTTTTTTTGGGCGGTTTTCTGCCGTCTCCGGGgaaaa                                                           >  1:564610/1‑90 (MQ=255)
                             aGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGAGCCACGGTTTTTTGCGCGTTTTTCCGCGTTCTCCGGCaaaaaa                                                          >  2:926506/1‑89 (MQ=255)
                              gTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGCCCCGGTTTTTTGCGCGTTTTTTTGCGTTCTCTTGCGAAAAcc                                                         >  1:757529/1‑90 (MQ=255)
                              gTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGAGCCACGTTTATTGGCGCGTTTTTCGGCGTTCTCGGGCaaaaaac                                                         >  2:332729/1‑88 (MQ=255)
                              gTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCCCGTTTATTGGCGCGTTTTTCGGCGTTCTCGGGGaaaaaac                                                         >  2:360944/1‑88 (MQ=255)
                              gTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTTTTGGCGCGTTTCTCGGCGTTCTCGGGCaaaaaac                                                         >  2:755974/1‑88 (MQ=255)
                                  ttgATGTCAGTACCTCTCATTCTGGCCATACTGGCGGGGGCAGCCACGTTTATTGGCGGGTTTCTCGGCGGTCGCGGGCCAAAACCCTCg                                                     <  1:340460/88‑1 (MQ=255)
                                     aTGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAAc                                                  <  2:603127/90‑1 (MQ=255)
                                           gTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTTTTGGGGCGTTTTTCGGCGTTCTCGGGCAAAAAACCTCGAACCCCTTa                                            >  2:485366/1‑90 (MQ=255)
                                               ctctCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTgg                                        <  1:360944/90‑1 (MQ=255)
                                                tctcATTCTGACCATACTGGCGGGGGGAGCCACGTTTATTTGGGCGGTTTTCGGCGTTCTCGGGCAAAAAACCCCCAACCGCCTACTGGc                                       >  2:862138/1‑90 (MQ=255)
                                                tctcATTCTGACCATACTGGCGGGGGCAGCCCCGTTTTTTTGCGCGGTTCTCGGCGGTCTCCGGCAAAAAACCCCCAACCCCTTTCTggg                                       >  1:196719/1‑89 (MQ=255)
                                                tctcATTCTGACCATACTGGCGGGGGCAGCCACGTTTTTTGGGGCGTTTTTCGGCGTTTTCGGGGAAAAACCCTCGAACCCCTTACTTGc                                       >  2:564245/1‑90 (MQ=255)
                                                         gACCATACTGGCGGGGGGAGCCCCGGTTATTTGCGCGGTTTTCGCCGTTCTCGGGGAAAAAACCCCGAACCCCATCTTCGCGTTTTCGCt                              >  1:920240/1‑90 (MQ=255)
                                                         gACCATACTGGCGGGGGGAGCCACGTTTTTTTGCGCGGTTTTCGGCGTTCTCGGGCAAAAAACCTCGAACCGCTTTCTGGGGGTTTTGCt                              >  2:487379/1‑90 (MQ=255)
                                                             aTACTGGCGGGGGCAGCCACGTTTTTTTGCGCGTTTTTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTTCTGGCGTTTTTGCTAggg                          >  2:396562/1‑89 (MQ=255)
                                                                   gCGGGGGCAGCCCCGTTTTTTTGCGCGTTTCTCCGCGTTCTCGGGCAAAAAACCCCCAACCCCCTTCTGGCGGTTTTGCTAGGGTTTtcg                    >  1:66187/1‑90 (MQ=255)
                                                                      ggggCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcggcg                 <  2:674960/90‑1 (MQ=255)
                                                                      ggggCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcggcg                 <  1:926506/90‑1 (MQ=255)
                                                                      ggggCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcggcg                 <  1:872835/90‑1 (MQ=255)
                                                                      ggggCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcggcg                 <  1:786987/90‑1 (MQ=255)
                                                                          cAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCggggg             >  2:415732/1‑89 (MQ=255)
                                                                                     aTTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGctc  <  1:561173/90‑1 (MQ=255)
                                                                                     aTTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGctc  <  1:755974/90‑1 (MQ=255)
                                                                                     |                                                                                         
GATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTC  >  NZ_CP009273/3175872‑3176046

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