Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927338,9430TG53.1% 11.8 / 28.6 32Y58S (TAC→TCC) caiBL‑carnitine CoA‑transferase
Reads supporting (aligned to +/- strand):  ref base T (2/13);  new base G (17/0);  total (19/13)
Fisher's exact test for biased strand distribution p-value = 3.02e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.72e-05
Rejected as consensus: Frequency below/above cutoff threshold.

TTTCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGGCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCC  >  NZ_CP009273/38854‑39008
                                                                                         |                                                                  
tttCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGGGGGt                                                                    >  2:267101/1‑90 (MQ=255)
             cAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGGTGCGGGGAGTTCGGTTGAAc                                                       >  1:1779268/1‑90 (MQ=255)
                   cGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAAt                                                 <  1:902767/90‑1 (MQ=255)
                   cGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAAt                                                 <  1:303323/90‑1 (MQ=255)
                          cgGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCg                                          <  1:1170643/90‑1 (MQ=255)
                               ttCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGGGGGGTGCGGGGAGTTTGGTTGGACGCGAAAGGGGTCGGgccc                                    >  1:1494889/1‑89 (MQ=255)
                                tCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGGCCCAg                                   <  2:1714403/90‑1 (MQ=255)
                                      ttGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGGGAGTTGCGGGGAGTTTGGGTGAACGCGAATGGGGTCGGGCCAGgcggcg                             >  1:521290/1‑90 (MQ=255)
                                      ttGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGGCCCAGGCGACg                             <  1:22524/90‑1 (MQ=255)
                                      ttGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGGCCCAGGCGACg                             <  1:1597048/90‑1 (MQ=255)
                                      ttGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGGAGGTCGGTTGAAAGCGAATGGGGTCGTCCCAGgcggcg                             >  1:1192187/1‑90 (MQ=255)
                                      ttGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGGTGCGGGGAGTTCGGTTGAACGCGAATGGGGTTGGCCCCGGCgggg                             >  2:1705145/1‑87 (MQ=255)
                                        gAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGGGAGGTTCGGGGAGTTTGGGTGAAAGCGGAAGGGGGTGGCCCCGGGGACggt                           >  2:1597048/1‑88 (MQ=255)
                                           aaTATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGGCCCAGGCGACGTtctc                        <  2:1145494/90‑1 (MQ=255)
                                            atatTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGGCCCAGGCGACGTTCTCg                       <  2:535397/90‑1 (MQ=255)
                                              atTTAACGACAGCGCGTGCAAATTGCGGCGGGGGGGTTGCGGGGAGGTCGGTTGAAACGGAAAGGGGTGGGGCCCGGCGGCGGTCTCgct                     >  2:32955/1‑88 (MQ=255)
                                                ttAACGACAGCGCGTGCAAATTGCGGCGGGGGAGGTGGGGGGAGTTTGGGTGGAAGGGGAAGGGGGCGGGCCCGGCGGCGGTCTCGGTcc                   >  2:1112229/1‑90 (MQ=255)
                                                ttAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGGCCCAGGCGACGTTCTCGATcc                   <  1:1082349/90‑1 (MQ=255)
                                                ttAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGGCCCAGGCGACGTTCTCGATcc                   <  1:32955/90‑1 (MQ=255)
                                                     gACAGCGCGTGCAAATTGCGGCGGGAGAGGTGCGGGGAGTTTGGTTGAAAGCGGAAGGGGTCGGGCCCGGGGGCGGTCTCGGTTCCGATa              >  2:563073/1‑90 (MQ=255)
                                                      aCAGCGCGTGCAAATTGCGGCGGGAGAGGTGCGGGGAGGTCGGGTGAACGCGAAAGGGGTCGGCCCAGGGGGCGGTCTCGGTCCAGATaa             >  1:1077876/1‑90 (MQ=255)
                                                        aGCGCGTGCAAATTGCGGCGGGGGAGGTGGGGGTAGTTTGGTTTGACGCGGAAGGGGTTGGGCCCGGGGGCGGTCTCGGTCCCGATaaat           >  2:230993/1‑88 (MQ=255)
                                                        aGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGGAGGTTGGTTTAACGCGAAAGGGGGCGGCCCAGGCGACGGTCTCGATCCAGATaaat           >  1:513057/1‑88 (MQ=255)
                                                        aGCGCGTGCAAATTGCGGCGGGAGAGGTGGGGGGAGTTTGGGTGAAAGCGAAAGGGGTCGGGCCCGGCGGCGGTCTCCGTCCAGAAAACt           >  2:446261/1‑90 (MQ=255)
                                                         gcgcGTGCAAATTGCGGCGGGAGGGGTGCGGGGAGTTTCGGTGAAAGCGAATGGTGTGGGCCCCAGGCACGCTCTCCATCCAGATAACtt          >  1:512724/1‑90 (MQ=255)
                                                         gcgcGTGCAAATTGCGGCGGGAGAGTTGCGGGGAGTTCGGTTGAACGCGAATGGGGGCGGCCCAGGCGGCGTTCTCGATCCAGATAACtt          >  1:1276085/1‑90 (MQ=255)
                                                         gcgcGTGCAAATTGCGGCGGGAGAGGTGGGGGGAGGTTGGTTGAAAGCGGATGGGGTGGGGCCCGGGGGCGGTCTCGATCCAGATAAAtt          >  2:1479271/1‑90 (MQ=255)
                                                         gcgcGTGCAAATTGCGGCGGGAGAGGTGCGGGGGGGTTGGGTGAAAGCGGAAGGGGTCGGGCCCGGCGACGTTCTCGATCCCGATAAAtt          >  2:430765/1‑90 (MQ=255)
                                                           gcgTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGGCCCAGGCGACGTTCTCGATCCAGATAACTTcc        <  1:1695219/90‑1 (MQ=255)
                                                              tGCAAATTGCGGCGGGAGAGGTGCGGGGAGTTCGGGTGAACGCGAAAGGGGTCGGGCCCGGCGGCGGTCTCGATCCCGATAAATTCcgcc     >  2:1793228/1‑89 (MQ=255)
                                                                 aaaTTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGGCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGccc  <  1:430765/90‑1 (MQ=255)
                                                                                agTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTc                                           <  1:392549/35‑1 (MQ=255)
                                                                                agTTGCGGGTAGTTCGGGTGAACGCGAATGGTGTc                                           >  2:392549/1‑35 (MQ=255)
                                                                                         |                                                                  
TTTCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGGCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCC  >  NZ_CP009273/38854‑39008

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