Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,141,2140TG72.2% 7.4 / 12.1 18G20G (GGT→GGGyegDmolecular chaperone
Reads supporting (aligned to +/- strand):  ref base T (0/5);  new base G (13/0);  total (13/5)
Fisher's exact test for biased strand distribution p-value = 1.17e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGA  >  NZ_CP009273/2141133‑2141297
                                                                                 |                                                                                   
tcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAAACGCa                                                                             >  1:286530/1‑90 (MQ=255)
   ccGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGGCATTCGGGACGGGAAAACGCAttt                                                                          >  2:593064/1‑90 (MQ=255)
               ttattTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCTTGTGTGTCGGGGAAACGCATTTTCTTAAAATAGa                                                              >  2:667191/1‑90 (MQ=255)
               ttattTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAACCGCATTTTCTAAAAAAGGa                                                              >  1:123532/1‑90 (MQ=255)
               ttattTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAACCGCATTTGCTAAAAAAGGa                                                              >  1:262800/1‑90 (MQ=255)
               ttattTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGGGACGGGAAAACGCATTTGGTAAAAAAGGa                                                              >  1:160794/1‑90 (MQ=255)
                      gtgtTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGGGGCGGGAAACCGCATTTGCTAAAAATGGAAAAAGAc                                                       >  1:534805/1‑90 (MQ=255)
                      gtgtTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGGGGCGGGAAAACGCATTTGCTAAAAAAGGAAAAAGAc                                                       >  1:529963/1‑90 (MQ=255)
                           tATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAACCGCATTTTCTAAAAAAGGAAAACGACAGCAc                                                  >  1:223160/1‑90 (MQ=255)
                                    tGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGcc                                         <  2:529963/90‑1 (MQ=255)
                                                      cTGTTCAGTGGCGGTCATGCGTGACGGGAAAACCCATTTTTTAAAAAAAGAAAAAAACAACACCCCGCTGCCTTCAATTCTTTTCGCGcc                       >  2:560951/1‑90 (MQ=255)
                                                        gTTCAGTGGCGGTCATGCGGGGCGGGAAAACGCCTTTGCTAAAAAAAGAAAAAGACAACACGCCGCTGCCCTCAATGCTTTTTGCGccca                     >  2:87983/1‑88 (MQ=255)
                                                          tCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCAAGCTGCTGCCTTCAATGCTTTGCGCGCCAAcg                   <  1:560951/90‑1 (MQ=255)
                                                             gTGGCGGTCATGCGTGGCGGGAAAACCCATTTTCTAAAAAAGGAAAAAGAAAGCACGCCGCTGCCTTCAATGCTTTTTGCGCCCAAGCGt                >  2:84234/1‑90 (MQ=255)
                                                                gCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGaa             <  1:84234/90‑1 (MQ=255)
                                                                  ggTCATGCGGGGCGGGAAACCGCATTTTCTAAAAATGGAAAAAGACAGCACGCCGCTGCCTTCAATGCTTTGCGCGCCAACGCCGGAAGc           >  2:73335/1‑90 (MQ=255)
                                                                      aTGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTa       <  1:415479/90‑1 (MQ=255)
                                                                           tGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGa  <  2:534805/90‑1 (MQ=255)
                                                                                 |                                                                                   
TCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGA  >  NZ_CP009273/2141133‑2141297

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