Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,175,9570AT46.8% 72.4 / 31.1 62I17F (ATT→TTT) zupTzinc transporter ZupT
Reads supporting (aligned to +/- strand):  ref base A (8/25);  new base T (29/0);  total (37/25)
Fisher's exact test for biased strand distribution p-value = 1.21e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGC  >  NZ_CP009273/3175881‑3176044
                                                                             |                                                                                        
gAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGCCACGTTTTTTTGCGCGttttt                                                                             >  2:1865910/1‑88 (MQ=255)
  aGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGCCCCCCTTTTTTGCGCGTTtttcc                                                                           >  1:1255344/1‑89 (MQ=255)
      gAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGCCCCGTTTTTTTGGGCGGTTTTCCGCGt                                                                       >  2:859994/1‑90 (MQ=255)
      gAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGCGCCCCGTTTTTTTGCGCGTTTTTCcgcgg                                                                       >  1:1166848/1‑89 (MQ=255)
      gAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCCCGTTTTTTTGCGCGTTTTTCCGCGt                                                                       >  1:673825/1‑90 (MQ=255)
                aaGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGCCCCCTTTTTTTGGGGGGTTTTTCGCGGTCTCggggca                                                            >  2:737050/1‑87 (MQ=255)
                  ggAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCaaaa                                                          <  1:687201/90‑1 (MQ=255)
                  ggAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCaaaa                                                          <  2:382184/90‑1 (MQ=255)
                  ggAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCaaaa                                                          <  1:191568/90‑1 (MQ=255)
                   gAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGCCACGTTTTTTTGCGCGGTTTTTCGCGTTCTCCGGCaaaaa                                                         >  2:646338/1‑90 (MQ=255)
                     gTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGCCCCGGTTTTTTGGGCGGTTTTCTGCGTTCTTCGGGAAAAAcc                                                       >  2:1054190/1‑90 (MQ=255)
                     gTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGAGCCCCCGTTTTTTGGGCGTTTTTTGGCGTTTTTGGGGaaaaaac                                                       >  2:1829186/1‑88 (MQ=255)
                     gTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGAGCCACGGTTTTTGGCGCGTTTTTCGGCGCTCTCCGGCCAAAAcc                                                       >  1:1921444/1‑90 (MQ=255)
                     gTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCGGCCGCGTTTTTTGGCGCGGTTTTTTGCGGTCTCTGGCaaaaaac                                                       >  2:1787607/1‑88 (MQ=255)
                     gTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCGGCCACGTTTTTTGGCGCGTTTTTCTGCGGTCTCGGGCAAAAAcc                                                       >  1:1371874/1‑90 (MQ=255)
                       tGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGGGGGCCCCTTTTTTTGGGCGTTTTTTTGCGCTGTCGGGGCAAAAAccc                                                      >  2:1660691/1‑90 (MQ=255)
                            aTGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAAc                                                <  2:420381/90‑1 (MQ=255)
                               tCAGTACCTCTCATTCTGACCATACTGGCGGGGGGAGCCACGTTTTTTTGGGCGTTTTTCGGCGTTCTCCGGCAAAAAACCCCCAAcccc                                             >  2:53932/1‑88 (MQ=255)
                               tCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGc                                             >  1:863447/1‑90 (MQ=255)
                                  gTACCTCTCATTCTGACCATACTGGCGGGGGGAGCCACGGTTTTTTGGGCGGTTTTCCGCGGTCTCCGGGAAAAAACCCCCAACCCCTTa                                          >  2:575181/1‑90 (MQ=255)
                                     ccTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTg                                       <  1:575181/90‑1 (MQ=255)
                                       tctcATTCTGACCATACTGGCGGGGGGGGCCACCCTTTTTTGGGCGGTTCTCGGCGGTCTCGGGCAAAAAACCCTGAACCCCTTAATGGc                                     >  1:276648/1‑90 (MQ=255)
                                       tctcATTCTGACCATACTGGCGGGGGGAGCCACGTTTTTTGGGGCGTTTTTCGGCGTTTTCGGGCCAAAAACCCCGAACCGCTTTCTGGc                                     >  1:396029/1‑90 (MQ=255)
                                       tctcATTCTGACCATACTGGCGGGGGGAGCCACGGTTATTTGGGGGGTTTTCCGCGTTTTCGGGCAAAAAACCCCGAACCGCTTCCTggg                                     >  1:1868904/1‑89 (MQ=255)
                                       tctcATTCTGACCATACTGGCGGGGGCAGCCACGCTTTTTTGGGCGTTTCTCTGCGGTCTCGGGCCAAAAACCCCCAACCCCCTACTGGc                                     >  1:2199321/1‑90 (MQ=255)
                                        ctcATTCTGACCATACTGGCGGGGGGAGCCACGTTTTTTTGCGCGTTTTTCGGCGGTCTCCGGCCAAAAACCCCCAACCCCCTACTGGCg                                    >  1:2036046/1‑90 (MQ=255)
                                          cATTCTGACCATACTGGCGGGGGGAGCCCCGGTTTTTTGCGCGGTTTTCGGCGGTCTCCGGGAAAAAACCCCGCACCCCGTTCTGGCggt                                  >  2:442150/1‑88 (MQ=255)
                                            ttCTGACCATACTGGCGGGGGCAGCCACGTTTTTTTGGGCGGTTTTCGGCGGTTTCGGGCAAAAAACCCCGAACCGCTTTCTTGCGgttt                                >  1:2001291/1‑90 (MQ=255)
                                            ttCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGtttt                                <  2:1265402/90‑1 (MQ=255)
                                                gACCATACTGGCGGGGGGGGCCACGTTTTTTTGGGCGTTTTTCTGCGTTTTCGGGCGAAAAACCCCCAACCGCTTACTGGGGGTTTTGCt                            >  2:1784363/1‑90 (MQ=255)
                                                gACCATACTGGCGGGGGGAGCCACGTTTTTTTGGGCGTTTCTCCGCGTTCTCCGGCCAAAAACCCCGCACCCCCTTCTGGCGgtttttct                            >  2:924760/1‑86 (MQ=255)
                                                gACCATACTGGCGGGGGCAGCCCCGTTTTTTTGCGCGTTTTTCGGCGTTTTCGGGCAAAAAACCTCGAACCGCTTTCTGGCGGTTTTGCt                            >  1:935800/1‑90 (MQ=255)
                                                gACCATACTGGCGGGGGCAGCCCCGGTTTTTTGCGCGGTTTTCGGCGTTCTCCGGCAAAAAACCCCGAACCCCTTTCTTGCGGTTTTGCt                            >  2:86053/1‑90 (MQ=255)
                                                 aCCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTa                           <  1:53932/90‑1 (MQ=255)
                                                 aCCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTa                           <  2:2199321/90‑1 (MQ=255)
                                                    aTACTGGCGGGGGGAGCCACGTTTTTTTGCGCGTTTTTCCGCGTTTTCCGGCAAAAAACCCCGAACCCCTTTCTGGGGGTTTTGCTAggg                        >  1:2216297/1‑89 (MQ=255)
                                                    aTACTGGCGGGGGCAGCCACGTTTTTTGGCGCGTTTTTCGGCGTTCTCGGGCAAAAACCCCCGAACCCCTTACTGGCGGTTTCGCTAGGt                        >  1:267636/1‑90 (MQ=255)
                                                     tACTGGCGGGGGCAGCCCCGGTTTTTTGCGCGGTTTTCCGCGTTCTCCGGGAAAAAACCCCGAACCGCCTTCTTGCGGTTTTGCTAgggt                       >  2:412844/1‑88 (MQ=255)
                                                        tGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTg                    <  1:787142/90‑1 (MQ=255)
                                                          gCGGGGGGAGCCACGTTTTTTTGCGGGTTTTTCCGCGGTCTCGGGCAAAAACCCCCCAACCCCCTTCTGGCGGTTTTGCTAGGGTTTgtg                  >  2:1030261/1‑88 (MQ=255)
                                                          gCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCCGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcg                  >  1:71478/1‑90 (MQ=255)
                                                           cGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcgg                 <  1:833434/90‑1 (MQ=255)
                                                            gggggCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcggc                <  2:1255344/90‑1 (MQ=255)
                                                             ggggCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcggcg               <  1:2215582/90‑1 (MQ=255)
                                                             ggggCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcggcg               <  2:863447/90‑1 (MQ=255)
                                                             ggggCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcggcg               <  2:771527/90‑1 (MQ=255)
                                                             ggggCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcggcg               <  1:1660691/90‑1 (MQ=255)
                                                             ggggCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcggcg               <  1:465781/90‑1 (MQ=255)
                                                             ggggCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcggcg               <  1:1030261/90‑1 (MQ=255)
                                                              gggCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGcggcgg              <  2:2064612/90‑1 (MQ=255)
                                                                  aGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGAt          <  1:824036/90‑1 (MQ=255)
                                                                  aGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGAt          <  1:312179/90‑1 (MQ=255)
                                                                  aGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGAt          <  1:1509617/90‑1 (MQ=255)
                                                                  aGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGAt          <  2:1166848/90‑1 (MQ=255)
                                                                  aGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGAt          <  1:1163164/90‑1 (MQ=255)
                                                                   gCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCggggggc         >  1:2214869/1‑87 (MQ=255)
                                                                   gCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCgggggac         >  2:1507178/1‑87 (MQ=255)
                                                                   gCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGGTc         >  1:852472/1‑90 (MQ=255)
                                                                        cGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGt     >  1:153756/1‑90 (MQ=255)
                                                                         gTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCaaaa                                                          <  1:1935084/36‑1 (MQ=255)
                                                                         gTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCaaaa                                                          >  2:1935084/1‑36 (MQ=255)
                                                                           ttATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGc  <  1:1954483/90‑1 (MQ=255)
                                                                             |                                                                                        
GAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGC  >  NZ_CP009273/3175881‑3176044

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