Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,141,2140TG57.9% 10.9 / 18.4 19G20G (GGT→GGGyegDmolecular chaperone
Reads supporting (aligned to +/- strand):  ref base T (2/6);  new base G (11/0);  total (13/6)
Fisher's exact test for biased strand distribution p-value = 1.03e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-03
Rejected as consensus: Frequency below/above cutoff threshold.

CCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATG  >  NZ_CP009273/2141126‑2141300
                                                                                        |                                                                                      
cccTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTa                                                                                       <  1:946071/90‑1 (MQ=255)
   ttcttcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGGGACGGGaaaa                                                                                    >  2:617570/1‑89 (MQ=255)
                      ttattTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGGCCTTGGTGCCGGGAAAACACCCATGCTaaaaaatga                                                                 >  1:790159/1‑86 (MQ=255)
                        attTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGaaa                                                               <  1:617570/90‑1 (MQ=255)
                                 ttATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCa                                                      <  2:603086/90‑1 (MQ=255)
                                         tttGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGGCGGGAAAACGCCTTTTCTAAAAATGGAAAAAGACAGCACGctgctg                                              >  1:63846/1‑90 (MQ=255)
                                                      cAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAAACGCATTTTCTAAAAAAAGAAAAAGACAACACGCCGCTGGCTTCAATGCttt                                 >  1:557548/1‑90 (MQ=255)
                                                             cTGTTCAGTGGCGGTCATGCGTGACGGGAAAACGCATTTGCTAAAAATGGGAAAAGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGcc                          >  2:29546/1‑90 (MQ=255)
                                                               gTTCAGTGGCGGTCATTCGGGGCGGGAAAACGCCTTTTTTAAAAAAAGAAAAAAACAACACGCCGCTGCCTTCCATTCTTTTTGCGCCaa                        >  1:419497/1‑90 (MQ=255)
                                                               gTTCAGTGGCGGTCATGCGTGGCGGGAAAACGCATTTGCTAAAAATGGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGCCaa                        >  2:726172/1‑90 (MQ=255)
                                                               gTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAAAGGAAAAAGACAGCACGCCGCTGCCTTCAATGCTTTGTGCGCCaa                        >  2:951758/1‑90 (MQ=255)
                                                                ttCAGTGGCGGTCATGCGTGGCGGGAAAACCCCTTTTTTAAAAAAGGAAAAAAAAAACACGCCGCTGCCCTCCAATCTTTTTGCGCCCAc                       >  1:925532/1‑90 (MQ=255)
                                                                ttCAGTGGCGGTCATGCGTGACGGGAAACCGCATTTGCTAAAAAAGGAAAAAGACAGCACCCTGCTGCCTTCAATGCTTTTCGCGCCAAc                       >  2:93829/1‑90 (MQ=255)
                                                                 tCAGTGGCGGGCATGCGTGACGGGAAACCGCATTTTCTAAAAAAGGAAAAAGAAAGCACGCGCCCGCCTTCAATGCTTTTTGCGCCAAcg                      >  1:484787/1‑90 (MQ=255)
                                                                    gTGGCGGTCATGCGTGGCGGGAAAACGCCTTTTCTAAAAATGGGAAAAAACAGCACGCGGCTGCCCTCCAAGTTTTTTGCGCCAAcgcct                   >  2:455465/1‑88 (MQ=255)
                                                                     tGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTg                  <  2:925532/90‑1 (MQ=255)
                                                                      ggCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGa                 <  2:1050722/90‑1 (MQ=255)
                                                                         ggTCATGCGTGACGGTAAACCGCATTTTCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGCCAACGCGTGaaac              >  1:769394/1‑88 (MQ=255)
                                                                                  tGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAATCGGTAAGCGa     <  1:951758/90‑1 (MQ=255)
                                                                                        tAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATg  <  2:673728/87‑1 (MQ=255)
                                                                                        tAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATg  >  1:673728/1‑87 (MQ=255)
                                                                                        |                                                                                      
CCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATG  >  NZ_CP009273/2141126‑2141300

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