Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,175,9460CG55.0% 5.6 / 17.6 20A13G (GCA→GGA) zupTzinc transporter ZupT
Reads supporting (aligned to +/- strand):  ref base C (5/4);  new base G (11/0);  total (16/4)
Fisher's exact test for biased strand distribution p-value = 2.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.33e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCG  >  NZ_CP009273/3175866‑3176031
                                                                                |                                                                                        
tttGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGAGCCACggt                                                                                 >  2:371220/1‑88 (MQ=255)
   gACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGCCACGttttt                                                                              >  1:112007/1‑88 (MQ=255)
             cTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGttt                                                                    <  2:45373/90‑1 (MQ=255)
                   gTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGAGCCCCGCTTTTTTGGGCGGTTCTCTGc                                                              >  1:708790/1‑90 (MQ=255)
                     gAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGAGCCCCGTTTTTTTGCGCGGTTTTCGCCgg                                                           >  1:485475/1‑89 (MQ=255)
                           ttATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCgg                                                   <  1:474619/90‑1 (MQ=255)
                                  gAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGCGCCCCGTTTTTTTGCGCGTTTTTCCGCGGTCTCCGGGCaaaa                                             >  2:378442/1‑90 (MQ=255)
                                    gTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGCCCCGTTTTTTTGCGCGGTTTTCTGCGGTTTCTGGGGAAAAAc                                           >  1:739163/1‑90 (MQ=255)
                                    gTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGAGCCACGTTTTTTTGCGCGTTTTTCGGCGTTCTCCGGGCAAAAAc                                           >  1:2002/1‑90 (MQ=255)
                                            tGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGCCCCGGTTTTTTGTGCGTTTTTTTGCGTTCTTTGGCAAAAAACCCCCAAcc                                  >  1:695125/1‑90 (MQ=255)
                                                  tACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTAc                            <  1:668297/90‑1 (MQ=255)
                                                      tctcATTCTGACCATACTGGCGGGGGGAGCCACGGTTTTTTGCGCGGTTTTCGGCGGTCTCCGGCAAAAAACCTCCCACCCCTTTCTTGc                        >  2:406684/1‑90 (MQ=255)
                                                      tctcATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAAACCTCGAACCGCTTACTggg                        >  2:216772/1‑89 (MQ=255)
                                                              tGACCATACTGGCGGGGGGAGCCACGGTTTTTTGCGCGGTTTTCCGCGGTTTCCGGGCAAAAACCCCCAACCCCCTTCTGGCGTTTTTGc                >  1:719552/1‑90 (MQ=255)
                                                               gACCATACTGGCGGGGGGCGCCACCTTTTTTTGCGCGTTTTTCGGCCTTCTCGGGGCAAAAACCCCCAACCGCCTACTTGGGGTTTTGCt               >  1:498894/1‑90 (MQ=255)
                                                               gACCATACTGGCGGGGGCAGCCCCGGTTGTTTGCGGGTTTTTCGTCGTTTTTGGGGGAAAAACCCCGAACCCCGTACTGTCGGTTTTGCt               >  1:559156/1‑90 (MQ=255)
                                                                   aTACTGGCGGGGGCAGCCACGGTTTTTTGGGCGTTTTTCGGCGTTCTCGGGCAAAAAACCTCGAACCGCTTACTGGCGGTTTTGCTAGGt           >  2:869848/1‑90 (MQ=255)
                                                                    tACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTTTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTTGCTAGGtt          >  2:253466/1‑90 (MQ=255)
                                                                         gCGGGGGCAGCCACGGTTTTTGGCGCGTTTTTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTTCTGGCGGTTTTGCTAGGGTTTtcg     >  2:203756/1‑90 (MQ=255)
                                                                            ggggCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcggcg  <  2:48535/90‑1 (MQ=255)
                                                                                |                                                                                        
TTTGACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCG  >  NZ_CP009273/3175866‑3176031

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