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

TGGCGCTGGGTACTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGG  >  NZ_CP009273/172675‑172833
                                                                                |                                                                               
tGGCGCTGGGTACTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTa                                                                        <  2:211702/90‑1 (MQ=255)
                gCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGc                                                       <  2:339981/90‑1 (MQ=255)
                 cTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTTCTATTTCCGGGATgggggga                                                      >  2:58375/1‑88 (MQ=255)
                 cTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTTCTATTGCCGGGATGGGGGCa                                                      >  1:215590/1‑90 (MQ=255)
                   gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGGGGTTTCTTTTGCCGGGATGGGGGCAtt                                                    >  1:563072/1‑90 (MQ=255)
                   gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTTCTATTGCCGGGATGGGGGGAtt                                                    >  2:206922/1‑90 (MQ=255)
                   gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGCTATTTCCGGGATGGGGGGAtt                                                    >  2:124256/1‑90 (MQ=255)
                   gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGCTATTGCCgg                                                                >  1:505793/1‑78 (MQ=255)
                   gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTGCTATTGCCGGGATGGGGGCAtt                                                    >  1:546671/1‑90 (MQ=255)
                   gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCgg                                                                <  2:505793/78‑1 (MQ=255)
                   gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGGTTCCGCAATATCACCTTGAGGCGGGGGGGTTTTCTTTTGTCCGGATTGGGGGctt                                                     >  1:329851/1‑88 (MQ=255)
                                           cGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGGTTGCTATTGCCGGAATGGGGGCATTACTGGCGGGATCTATTTGCGCGcc                            >  1:72236/1‑90 (MQ=255)
                                             ttGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCACCTATTCGCGCGCCGt                          <  2:425499/90‑1 (MQ=255)
                                                      ttCCGCAATATCACCTTGAGGCGGGGGGGGTTGCTATTGCCGGGATGGGGGGATTTCTTGGGGGCTCTATTTGCGCGCCCTTTACggggg                 >  1:282492/1‑89 (MQ=255)
                                                      ttCCGCAATATCACCTTGAGGCGGGGGCGGTTGCTATTGCCGGGATGGGGGGATTACTTGCGGGATCTATTTTCGCGCCGGTTACggggg                 >  1:74228/1‑89 (MQ=255)
                                                       tCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGAt                <  1:434912/90‑1 (MQ=255)
                                                             atatCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCAttct          <  2:563072/90‑1 (MQ=255)
                                                                  aCCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCAttctggttc     <  2:175799/90‑1 (MQ=255)
                                                                     ttGAGGCGGGGGCGTTTGCTATTGCCGGGATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGGTCAttttggttttgg  >  1:335467/1‑90 (MQ=255)
                                                                                |                                                                               
TGGCGCTGGGTACTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGG  >  NZ_CP009273/172675‑172833

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