Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,175,9570AT54.5% ‑0.5 / 10.4 11I17F (ATT→TTT) zupTzinc transporter ZupT
Reads supporting (aligned to +/- strand):  ref base A (0/5);  new base T (6/0);  total (6/5)
Fisher's exact test for biased strand distribution p-value = 2.16e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGC  >  NZ_CP009273/3175884‑3176030
                                                                         |                                                                          
ggTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCgg                                                            <  2:270127/90‑1 (MQ=255)
          tataAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGCGCCCCGGTTTTTTGCGCGGTTTTCTGCGTTCTCggg                                                 >  2:198473/1‑90 (MQ=255)
               ggAGTGGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCCCGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCaaaa                                            <  1:209622/90‑1 (MQ=255)
                gAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGGTTTTTTGGGCGGTTTTCGGCGTTCTCCGGCaaaaa                                           >  2:198507/1‑90 (MQ=255)
                  gTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGCCCCGTTTTTTTGCGCGGTTTTCGGCGGTTTCCGGGCAAAAAc                                          >  1:307308/1‑90 (MQ=255)
                  gTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGCCACGTTTTTTTTGGCGTTTTTCGGCGGTCTCGGGGgaaaaac                                         >  1:355446/1‑88 (MQ=255)
                           gTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCg                                <  1:198473/90‑1 (MQ=255)
                            tCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTTTTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAAACCTCGAACCGc                               >  2:160692/1‑90 (MQ=255)
                               gTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCCCGGTTTTTGGCGCGTTTCTCCGCGTTCTCGGGCAAAAACCCCCGAACCGCTTa                            >  2:17805/1‑90 (MQ=255)
                                              aCCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTa             <  1:198507/90‑1 (MQ=255)
                                                         gggggCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcggc  <  1:105564/90‑1 (MQ=255)
                                                                         |                                                                          
GGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGC  >  NZ_CP009273/3175884‑3176030

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