Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,245,9890CT55.0% ‑0.6 / 15.8 20E136K (GAA→AAA) plsBglycerol‑3‑phosphate 1‑O‑acyltransferase PlsB
Reads supporting (aligned to +/- strand):  ref base C (1/8);  new base T (11/0);  total (12/8)
Fisher's exact test for biased strand distribution p-value = 7.14e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTAC  >  NZ_CP009273/4245913‑4246057
                                                                            |                                                                    
cgACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGgcgcg                                                         <  1:196798/90‑1 (MQ=255)
 gACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGgcgcgc                                                        >  2:155271/1‑90 (MQ=255)
 gACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCCTTTTTACCCCCCGgcgcgc                                                        >  2:626421/1‑90 (MQ=255)
      aGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCa                                                   <  1:272555/90‑1 (MQ=255)
           caGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACAt                                              <  1:458139/90‑1 (MQ=255)
             gTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACCTTCCCTTTTTTACCCCCCCGCGCGCGCCCCAACaacc                                            >  1:220717/1‑87 (MQ=255)
                cAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCg                                         <  2:570523/90‑1 (MQ=255)
                cAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCTCCCTTTTTTCCCCCCCGCGCGCGACCAAACaccaccc                                         >  1:618376/1‑89 (MQ=255)
                cAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTTTACGCCCCGGCGCGCGGCCAAACATCACCg                                         >  2:255053/1‑90 (MQ=255)
                  gCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTTCCTTTGCCTTTTTTTCCCCCCCGCGCGCGGCCAAAAaccacccgc                                       >  1:370037/1‑87 (MQ=255)
                  gCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTTTACCCCCCCGCGCGCGACCAAAAAACACCGac                                       >  1:224782/1‑90 (MQ=255)
                  gCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCCCCCTTTTTTCCCCCCCGCGCGCGGCCAAAAATCAccccc                                       >  1:356462/1‑87 (MQ=255)
                    gAAAAATTTCTGTCCGCCGTTAAGCATACGCAGCGGCGGGTTCCCTTCGCCTTTTTCCCGCCCCGGCGCGCGACCAAACATCACCGacac                                     <  1:588873/90‑1 (MQ=255)
                       aaaTTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTTTACGCCCCCGCGCGCGGCCCAACCACACCCACACCgg                                  >  2:375025/1‑90 (MQ=255)
                                      gTTAAGCATACGCAGCGGCGGGTTTCCTTCGCCTTTTTTACGCCCCCGCGCGCGACCCAACAACACCCGCACTGGCCCCCCCTGCACCTc                   >  1:64395/1‑90 (MQ=255)
                                            cATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGAtt             <  1:255053/90‑1 (MQ=255)
                                             aTACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGAttt            <  1:227237/90‑1 (MQ=255)
                                              tACGCAGCGGCGGGTTCACTTCGCCTTTTTTACGCCCCCGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTg           >  2:100464/1‑90 (MQ=255)
                                                    gcggcgGGTTCACTTCGCCTTTTTTACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTTGGGTGc     >  2:148217/1‑90 (MQ=255)
                                                       gcggGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTAc  <  1:155271/90‑1 (MQ=255)
                                                                            |                                                                    
CGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTAC  >  NZ_CP009273/4245913‑4246057

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