Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,504,6520CA55.6% 3.3 / 11.1 9S131R (AGC→AGA) ydcRPLP‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base C (4/0);  new base A (0/5);  total (4/5)
Fisher's exact test for biased strand distribution p-value = 7.94e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.17e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGC  >  NZ_CP009273/1504568‑1504735
                                                                                     |                                                                                   
gtcggtcgTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAAc                                                                                >  1:506940/1‑90 (MQ=255)
       gTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCcaca                                                                         >  1:437050/1‑90 (MQ=255)
                   cccgCCTTTCCCCCCCCCCCTCTTTCCCCCCCCCAAACAACAAACCCGCCCGCGGGCGCAGGAAAAAAAACCCGCCACAGCGATGAGCGt                                                              <  2:242187/87‑1 (MQ=255)
                          ttCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGaaa                                                      >  2:15705/1‑90 (MQ=255)
                                                         aacAAACCCGCCCGCTGGCGCAGGAAAGAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTc                        <  1:132341/90‑1 (MQ=255)
                                                                ccgcccgcGGGCGAAGGAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTAt                 <  1:546039/90‑1 (MQ=255)
                                                                          gcgcAGGAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGc       <  1:149424/90‑1 (MQ=255)
                                                                            gcAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTa     >  1:308608/1‑90 (MQ=255)
                                                                               ggaaaGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGc  <  1:115113/87‑1 (MQ=255)
                                                                                     |                                                                                   
GTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGC  >  NZ_CP009273/1504568‑1504735

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑