Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,175,9600GT54.5% 17.3 / 31.9 22G18C (GGC→TGC) zupTzinc transporter ZupT
Reads supporting (aligned to +/- strand):  ref base G (4/6);  new base T (12/0);  total (16/6)
Fisher's exact test for biased strand distribution p-value = 2.81e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-04
Rejected as consensus: Frequency below/above cutoff threshold.

TCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATC  >  NZ_CP009273/3175874‑3176049
                                                                                       |                                                                                          
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)
                                                                                     ttGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGGGCATGTTGctcc    >  1:249481/1‑89 (MQ=255)
                                                                                     ttGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGctca    >  2:283860/1‑90 (MQ=255)
                                                                                       ggCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGctcatc  <  2:119154/90‑1 (MQ=255)
                                                                                       ggCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGctcatc  <  1:213397/90‑1 (MQ=255)
                                                                                       ggCGCGTTTCCCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGctcatc  <  2:347305/90‑1 (MQ=255)
                                                                                       |                                                                                          
TCTGCCTGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATC  >  NZ_CP009273/3175874‑3176049

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