Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,175,9570AT59.1% 8.3 / 25.1 23I17F (ATT→TTT) zupTzinc transporter ZupT
Reads supporting (aligned to +/- strand):  ref base A (1/8);  new base T (13/0);  total (15/8)
Fisher's exact test for biased strand distribution p-value = 2.81e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.16e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCAT  >  NZ_CP009273/3175871‑3176039
                                                                                      |                                                                                    
cGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGAGCCCCGGTttttt                                                                                   >  2:204761/1‑89 (MQ=255)
    cTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGgcgc                                                                               <  2:590587/90‑1 (MQ=255)
     tGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGcgcg                                                                              <  1:553720/90‑1 (MQ=255)
          gAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGTAGCCGCGGTTTTTTGGGCGGTtct                                                                         >  2:711569/1‑90 (MQ=255)
                gAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGAGCCACGGTTTTTGGGGCGGTTTTTGGCGt                                                                   >  1:590535/1‑90 (MQ=255)
                               gTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCGGGCACGTTTGTTTGCGCCTTTTTCTCCGCTCTCCGCCCAAAAcc                                                  >  2:244574/1‑90 (MQ=255)
                               gTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTTTTGGCGCGTTTTTCGGCGTTCTCGGGCaaaaaac                                                  >  1:59105/1‑88 (MQ=255)
                                 tGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGCCACGGTTTTTTGCGCGCTTTTCTGCGTTTTCCGGCGAAAAAccc                                                 >  1:701177/1‑90 (MQ=255)
                                      aTGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAAc                                           <  2:54966/90‑1 (MQ=255)
                                       tGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGAGCCCCGGTTTTTGGGGGGGTTTTCCGCGTTCTCCGGGAAAAAACCCCCAAcc                                          >  2:43242/1‑90 (MQ=255)
                                            gTACCTCTCATTCTGACCATACTGGCGGGGGGAGCCACGTTTTTTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAAACCTCCAACCGCTTa                                     >  2:281237/1‑90 (MQ=255)
                                                ctctCATTCTGACCATACTGGCGGGGGCAGCCACGTTTTTTTGCGCGTTTTTCCGCGTTTTCCGGCAAAAAACCCCCAACCCCCTTCTgg                                 >  1:761210/1‑90 (MQ=255)
                                                 tctcATTCTGACCATACTGGCGGGGGGGGCCCCGGTTTTTTGCGCGGTTCTCTCCGCTCTCGGGGGAAAAACCCCCCAACCCCTTCTGGc                                >  2:586247/1‑90 (MQ=255)
                                                    cATTCTGACCATACTGGCGGGGGGAGCCACGGTTTTTTGCGCGGTTTTCCGCGGTCTCCGGCAAAAAACCCCCAACCCCTTACTGGCggt                             >  2:413402/1‑88 (MQ=255)
                                                         tGACCATACTGGCGGGGGCAGCCACGTTTATTTGGGCGTTTTTTCGCGTTCTCGGGCAAAAAACCCCGAACCCCTTTCTGGCGGTTTTGc                        >  1:596059/1‑90 (MQ=255)
                                                          gACCATACTGGCGGGGGCAGCCCCGTTTTTTTGGGCGGTTTTCGGCGTTTTCGGGCAAAAAACCCCCAACCCCTTTCTGGGGGTTTTGCt                       >  1:955512/1‑90 (MQ=255)
                                                          gACCATACTGGCGGGGGCAGCCACGTTTTTTTGGGCGTTTTTCGGCGTTCTCGGGCAAAAAACCTCGAACCGCTTTCTGGGGTTTTTGCt                       >  1:366031/1‑90 (MQ=255)
                                                          gACCATACTGGCGGGGGCAGCCACGTTTTTTTGCGCGTTTCTCGGCGTTTTCGGGCAAAAAACCCCGAACCGCTTTCTGGCGTTTTTGCt                       >  1:257943/1‑90 (MQ=255)
                                                           aCCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTa                      <  1:586247/90‑1 (MQ=255)
                                                               tACTGGCGGGGGCGGCCACGTTTTTTTGCGCGTTTCTCTGCGTTTTCTGGCAAAAACCCCCGAACCGCCTACTGGCGGTTTTGCTAGGtt                  >  1:191106/1‑90 (MQ=255)
                                                                     cGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcgg            <  2:257943/90‑1 (MQ=255)
                                                                       ggggCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcggcg          <  2:191106/90‑1 (MQ=255)
                                                                       ggggCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcggcg          <  2:988624/90‑1 (MQ=255)
                                                                               cACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCAt  <  2:596059/90‑1 (MQ=255)
                                                                                      |                                                                                    
CGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCAT  >  NZ_CP009273/3175871‑3176039

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