Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273172,7550AG57.1% 3.2 / 19.7 21T388A (ACG→GCG) clcAH(+)/Cl(‑) exchange transporter ClcA
Reads supporting (aligned to +/- strand):  ref base A (2/7);  new base G (12/0);  total (14/7)
Fisher's exact test for biased strand distribution p-value = 3.10e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GATGCTGGCGCTGGGTACTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGACC  >  NZ_CP009273/172670‑172841
                                                                                     |                                                                                       
gATGCTGGCGCTGGGTACTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGGGCTGTTTCCGCAATATCACCTTGAGGCGGGGAGGtt                                                                                     >  1:402085/1‑90 (MQ=255)
              gTACTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGaa                                                                       <  2:279465/90‑1 (MQ=255)
                aCTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATg                                                                     <  2:224397/90‑1 (MQ=255)
                 cTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTTGTTTTGCCGGGATgg                                                                    >  1:90856/1‑90 (MQ=255)
                   gtgCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTGCTTTTGCCGGGATtggg                                                                  >  1:246415/1‑90 (MQ=255)
                        gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGGGTTTTGTTTTGCCCGGATGGGGGGAtt                                                             >  1:82330/1‑90 (MQ=255)
                             ccGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTGCTATTGCCGGAATAGGGGTATTACTgg                                                        >  1:176367/1‑90 (MQ=255)
                             ccGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGAGTTTGCTTTTGTCGGGAAGGGGGGATTaaccg                                                        >  2:242406/1‑86 (MQ=255)
                               gCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTggcg                                                      >  2:102789/1‑90 (MQ=255)
                                         tGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTGCTATTTCCGGGATGGGGGGATTACTGGCGGCCTCTATTc                                            >  1:499914/1‑90 (MQ=255)
                                              gCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGGGTTTTCTTTTTCCGGGATGGGGGGATTACTGGCGGCGTCTTTTtccgcg                                       >  1:360481/1‑90 (MQ=255)
                                                cGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTGCTATTTCCGGGATGGGGGGATTAATGGCGGCATCTATTTCCGCGcc                                    >  1:388737/1‑90 (MQ=255)
                                                    gAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGCTTTTTCCGGGATGGGGGGATTACTTGCGGCATCTATTCGCGCGCCGttt                                >  1:436760/1‑89 (MQ=255)
                                                      gCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGCTATTGCCGGGATGGGGGCATTACTTGCGGCATCTATTTGCGCGCCGTTTAc                              >  1:171073/1‑90 (MQ=255)
                                                           ttCCGCAATATCACCTTGAGGCGGGGGCGTTTGTTTTTTCCGGGGTGGGGGGGTTTCTTGCGGGCTCTATTTGCGCGCCCGTTAAGgggg                          >  1:235934/1‑90 (MQ=255)
                                                                  atatCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCAttct                  <  2:176367/90‑1 (MQ=255)
                                                                  atatCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCAttct                  <  1:224786/90‑1 (MQ=255)
                                                                  atatCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCAttct                  <  2:499914/90‑1 (MQ=255)
                                                                  atatCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCAttct                  <  2:86685/90‑1 (MQ=255)
                                                                            gAGGCGGGGGCGTTTGCTATTGCCGGGATGGGGGGATTTTTTGCGGCATCTATTTGCGCGCCGGTAACGGGGATCATTCTTGGTCTGgag        >  2:341094/1‑90 (MQ=255)
                                                                                  gggACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGAcc  <  1:421324/90‑1 (MQ=255)
                                                                                     |                                                                                       
GATGCTGGCGCTGGGTACTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGACC  >  NZ_CP009273/172670‑172841

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