Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273172,7720AG52.9% 8.4 / 16.4 17G393G (GGA→GGGclcAH(+)/Cl(‑) exchange transporter ClcA
Reads supporting (aligned to +/- strand):  ref base A (2/6);  new base G (9/0);  total (11/6)
Fisher's exact test for biased strand distribution p-value = 2.26e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGACCGATAACTACCAGC  >  NZ_CP009273/172691‑172854
                                                                                  |                                                                                  
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)
   gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTGCTATTGCCGGGATGGGGGCAtt                                                                         >  1:546671/1‑90 (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)
                                                                          ttGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGACCGATAACTACCAGc  <  2:329851/90‑1 (MQ=255)
                                                                                  |                                                                                  
GCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGACCGATAACTACCAGC  >  NZ_CP009273/172691‑172854

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: