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

GACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGT  >  NZ_CP009273/3175869‑3176041
                                                                                         |                                                                                     
gACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCCGGGGCAGCCACGTTTAt                                                                                      >  1:715965/1‑90 (MQ=255)
     tCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGGCCCGGTTATTgtgc                                                                                 >  1:508989/1‑87 (MQ=255)
      cTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGgcgc                                                                                <  2:124002/90‑1 (MQ=255)
            gAAGGTGAAGATTTACAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGCCACGCTTTTTTGGGCGgtttt                                                                          >  1:622498/1‑88 (MQ=255)
                                 gTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGAGCCCCGCTTTTTTGGGCGGTTTTCGGCGTTCTCGGGCAAAAAcc                                                    >  2:377717/1‑90 (MQ=255)
                                 gTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTTTTGGCGCGTTTTTCGGCGTTCTCCGGCAAAAAcc                                                    >  2:151730/1‑90 (MQ=255)
                                         tGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGCCACGTTTTTTTGGGCGGTTTTCCGCGTTCTCTGGGGAAAAACCCCGAAcc                                            >  1:592784/1‑90 (MQ=255)
                                         tGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCCCGGTTTTTTGCGCGTTTTTCCGCGTTTTCGGGCAAAAAACCCCCAAcc                                            >  2:676890/1‑90 (MQ=255)
                                              gTACCTCTCATTCTGACCATACTGGCGGGGGCGGCCCCGGTTTTTTGCGCGGTTTTCCGCGTTCTCGGGCAAAAAACCCCGAACCGCTTa                                       >  1:30385/1‑90 (MQ=255)
                                                   tctcATTCTGACCATACTGGCGGGGGGGGCCCCGTTTTTTTGCGCGGTTCTCGGCGTTCTCGGGCAAAAAACCCCGAACCGCCTTCTggg                                  >  1:256350/1‑89 (MQ=255)
                                                   tctcATTCTGACCATACTGGCGGGGGGAGCCCCGTTTTTTTGCGCGGTTTTCGGCGTTCTCCGGCAAAAAACCCCGAACCCCTTACTGGc                                  >  2:592041/1‑90 (MQ=255)
                                                   tctcATTCTGACCATACTGGCGGGGGCAGCCACGTTTTTTTGCGCGTTTCTCGGCGTTCTCGGGCAAAAAACCCCGAACCGCTTACTGGc                                  >  2:264175/1‑90 (MQ=255)
                                                   tctcATTCTGACCATACTGGCGGGGGCAGCCACCTTTTTTTGGGCGGTTTTCCGCCTTCTCGGGCAAAAACCCCCGAACCCCCTACTggg                                  >  2:730683/1‑89 (MQ=255)
                                                          cTGACCATACTGGCGGGGGCAGCCCCGGTTTTTTGCGCGTTTTTCGGCGTTCTCTGGGCAAAAACCCCCAACCGCTTACTGGCGgttttg                           >  1:626483/1‑88 (MQ=255)
                                                             aCCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTa                        <  2:637371/90‑1 (MQ=255)
                                                                 tACTGGCGGGGGCAGCCACGTTTTTTTGCGCGTTTTTCGGCGTTTTCGGGCAAAAACCCCCCAACCCCCTACTGGCGGTTTCGCTAgggt                    >  2:663822/1‑88 (MQ=255)
                                                                 tACTGGCGGGGGCAGCCACGTTTATTTGCGCGGTTTTCGGCGTTCTCGGGCAAAAAACCCCGAACCGCTTTCTGGGGGTTTTGCCAGGtt                    >  2:745455/1‑90 (MQ=255)
                                                                      gCGGGGGGAGCCACCGTTTTTTGGGCGTTTCTCGGCGGTCTCCGGGAAAAAACCCCCAACCCCCTTCTTGCGGTTTTGCTAGGGTTTtcg               >  1:174751/1‑90 (MQ=255)
                                                                         ggggCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcggcg            <  2:554611/90‑1 (MQ=255)
                                                                            gCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCgggg         <  2:662235/90‑1 (MQ=255)
                                                                              aGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGAt       <  1:377717/90‑1 (MQ=255)
                                                                                    cGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGt  <  1:151730/90‑1 (MQ=255)
                                                                                         |                                                                                     
GACGATCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGT  >  NZ_CP009273/3175869‑3176041

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