Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,175,9460CG55.6% 13.3 / 26.6 18A13G (GCA→GGA) zupTzinc transporter ZupT
Reads supporting (aligned to +/- strand):  ref base C (0/8);  new base G (10/0);  total (10/8)
Fisher's exact test for biased strand distribution p-value = 2.29e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-04
Rejected as consensus: Frequency below/above cutoff threshold.

AGTTATTTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCG  >  NZ_CP009273/3175860‑3176028
                                                                                      |                                                                                    
aGTTATTTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGc                                                                                   <  2:42906/90‑1 (MQ=255)
  ttATTTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCa                                                                                 <  1:150465/90‑1 (MQ=255)
  ttATTTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCa                                                                                 <  1:438142/90‑1 (MQ=255)
           cGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGAGGCACGGTttttt                                                                       >  1:328929/1‑89 (MQ=255)
               cTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGCCACGTTTTTTTgcgg                                                                   >  2:567331/1‑89 (MQ=255)
               cTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGgcgc                                                                   <  1:449917/90‑1 (MQ=255)
               cTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGgcgc                                                                   <  1:91559/90‑1 (MQ=255)
                       aGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGAGCCCCGCTTTTTTGGGCGTTtcttg                                                           >  2:377230/1‑88 (MQ=255)
                           gAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGGCACGGTTTTTTGGGGGGTTTTTCgcgt                                                        >  1:68054/1‑87 (MQ=255)
                                       ggAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCaaaa                                          <  2:48517/90‑1 (MQ=255)
                                        gAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGAGCCACGCTTATTTGCGCGGTTTTCCCCGTTCTCCGGCaaaaa                                         >  2:500391/1‑90 (MQ=255)
                                          gTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGCCCCGGTTTTTTGGGCGGTTTTTCGCGTTCTCCGGGCAAAAAc                                        >  1:71686/1‑90 (MQ=255)
                                          gTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGCCACGTTTTTTTGCGCGGTTTTCGGCGGTCTCGGGGgaaaaac                                       >  1:5344/1‑88 (MQ=255)
                                           ttGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGAGCCACGTTTTTTTGCGCGGTTTTCGGCGTTCTCGGGGAAAAAacc                                      >  1:127068/1‑90 (MQ=255)
                                               cGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGa                                  <  2:258026/90‑1 (MQ=255)
                                               cGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGa                                  <  1:299124/90‑1 (MQ=255)
                                                             ctcATTCTGACCATACTGGCGGGGGGCGCCCCGTTTATTTGGGCGGTTTTCTGCGTTCTCCGGGAAAAAACCTCCAAACGCTTACTgggg                    >  1:537865/1‑88 (MQ=255)
                                                                               gCGGGGGGAGCCACGTTTTTTTGGGCGTTTTTCGGCGTTCTCGGGCAAAAAACCCCGAACCGCTTTCTGGGGGTTTTGCTAGGTTTTtcg  >  1:45518/1‑90 (MQ=255)
                                                                                      |                                                                                    
AGTTATTTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCG  >  NZ_CP009273/3175860‑3176028

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