Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,621,9840AC54.5% 4.7 / 11.0 11V105G (GTC→GGC) pdeFcyclic‑guanylate‑specific phosphodiesterase PdeF
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base C (0/6);  total (3/8)
Fisher's exact test for biased strand distribution p-value = 6.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGC  >  NZ_CP009273/2621898‑2622061
                                                                                      |                                                                             
ttGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGAccc                                                                            <  2:22171/90‑1 (MQ=255)
               tCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCCCCATGTGAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAAc                                                             <  2:7999/90‑1 (MQ=255)
                 tGCCAGAATATACGCTGGGAAATTAAACGGGTATCCCCATGTGAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACag                                                           <  1:452102/90‑1 (MQ=255)
                 tGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCAAACag                                                           <  1:437091/90‑1 (MQ=255)
                                tGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAatgataatgata                                            >  1:241476/1‑90 (MQ=255)
                                 gggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAAtgataatgataa                                           >  1:520506/1‑90 (MQ=255)
                                 gggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAAtgataatgataa                                           >  2:204676/1‑90 (MQ=255)
                                           aggggTATCCCCATGTAAAACGTTTTTACGCGGGGGAGAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGc                                 <  2:475208/88‑1 (MQ=255)
                                                      cttgtaAAACGTTTTTCCGCGGGGGAGAAAAGCCCCTTTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATa                      <  2:242606/88‑1 (MQ=255)
                                                                     tccGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAc       <  2:235494/88‑1 (MQ=255)
                                                                          cGGGGAAAAAAGCCCCTGTACCCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGc  <  1:469436/90‑1 (MQ=255)
                                                                                      |                                                                             
TTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGC  >  NZ_CP009273/2621898‑2622061

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