Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092736,9940AC58.3% 13.6 / 24.9 24G322G (GGT→GGGyaaJsodium:alanine symporter family protein
Reads supporting (aligned to +/- strand):  ref base A (9/1);  new base C (0/14);  total (9/15)
Fisher's exact test for biased strand distribution p-value = 7.65e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.04e-03
Rejected as consensus: Frequency below/above cutoff threshold.

CTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCC  >  NZ_CP009273/6910‑7074
                                                                                    |                                                                                
cTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGc                                                                             >  2:334922/1‑90 (MQ=255)
 tCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGcc                                                                            >  1:627848/1‑90 (MQ=255)
          aaGGAACCCATTCCCCCCCCCATCCCCTTCTGAAAAAGCTGAATCCCTCCCAGCGGCTTGTGGTTTGTGCCTTTCCCCGCCAGTAATATc                                                                   <  1:59687/89‑1 (MQ=255)
                 ccATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATgg                                                            >  1:38416/1‑90 (MQ=255)
                   aTTAGCCCCCGCACCCCCTTCTGGATAAGCTAAATCCCTCCCAGCGGCATGTGGGTTGTGCCGTTCCCCGCCAGTAATATCAGCATGgcg                                                          <  2:167754/90‑1 (MQ=255)
                      aGCCCCCGCCTCCCCTTGTGAAAAACCTAAATCCCTCCCGGCGGCATGTGGGTTGTGCCGTTCCCCGCCAGTAATATCAGCATGGCGCtt                                                       <  1:437962/90‑1 (MQ=255)
                          cccGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCg                                                   >  1:173723/1‑90 (MQ=255)
                            cGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTa                                                 >  2:112327/1‑90 (MQ=255)
                             gCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTAc                                                >  1:854644/1‑90 (MQ=255)
                               ccccccTTGTGGAAAACCAAAATCCCTCCCGGGGGCATGGGGGTTGTGCCTTCCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAg                                              <  1:850644/86‑1 (MQ=255)
                                  cccTTCTGAAAAAGCTAAATCCCTCCCGGCGGCTTGTGGGTTGTGCCTTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATg                                           <  1:313952/89‑1 (MQ=255)
                                     ttCTGAAAAAGCTAAATCCCTCCCGGCGGCATGTAGTTTGTGCCGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGAcc                                        <  2:367637/90‑1 (MQ=255)
                                                      cccTCCCAGCGGCATGTAGGTTGTCCCTTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAAt                       <  2:166450/89‑1 (MQ=255)
                                                           ccggCGGCTTGTGGGTTGTGCCTTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCaa                  <  2:854644/87‑1 (MQ=255)
                                                           ccAGGGGCGTGTGGTTTGTGCCTTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCaa                  <  2:626242/90‑1 (MQ=255)
                                                           ccAGCGGAATGTAGGTTGTGCCTTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCaa                  <  1:250989/90‑1 (MQ=255)
                                                                ggCTTGTGGTTTGTGCCTTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAatcatc             <  2:828219/90‑1 (MQ=255)
                                                                ggCCTGTAGGTTGTGCCTTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAatcatc             <  2:698772/90‑1 (MQ=255)
                                                                ggCATGTGGGTTGTGCCTTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAatcatc             <  1:501728/90‑1 (MQ=255)
                                                                ggCAGGTAGTTTGTCCCGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAatcatc             <  2:318953/90‑1 (MQ=255)
                                                                 gCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCt            >  1:828219/1‑90 (MQ=255)
                                                                           tgtgCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATccc  >  1:534227/1‑90 (MQ=255)
                                                                            gtgCCTTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAggg                                    <  1:493169/55‑1 (MQ=255)
                                                                            gtgCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAggg                                    >  2:493169/1‑55 (MQ=255)
                                                                                    |                                                                                
CTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCC  >  NZ_CP009273/6910‑7074

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: