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

GTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCT  >  NZ_CP009273/2621897‑2622062
                                                                                       |                                                                              
gTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGAcc                                                                              >  2:488302/1‑90 (MQ=255)
             caATCTGCAAGAATTTGCGCGGGAAAATTAAGCGGGTACCCCAATGTAAAATGTTGTTCCGCCGGGGAAAAAAGCCCCTGTAACCGcccc                                                                 <  2:566160/90‑1 (MQ=255)
                tCTGCCAGAATATACGCTGGAAAATTAACCGGGTATCCCCATGTGAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAAc                                                              <  1:115143/90‑1 (MQ=255)
                    ccAAAATATCCGGGGGAAAATTAAACGGGTATCACCATGTAAAACGTTTTTCCGCGGGGGAGAAAAGCCCCTGAAACCGCCCCAACagag                                                          <  2:99115/90‑1 (MQ=255)
                    ccAAAATATACGCGGGAAAATAAACCGGTTACCCCCATGTAAAACGTTGTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACagag                                                          <  2:573163/90‑1 (MQ=255)
                           ataCGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaatgata                                                   >  1:35819/1‑90 (MQ=255)
                                  gggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAAtgataatgataa                                            >  2:252995/1‑90 (MQ=255)
                                            aCGGGTACCACCATGTAAAACGTTGTTACGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGc                                  <  1:449178/90‑1 (MQ=255)
                                              gggTACCCCCTTGTAAAACTTTTTTACGCGGGGGAAAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCg                                <  2:2378/90‑1 (MQ=255)
                                                         tgtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGc                     >  1:206249/1‑90 (MQ=255)
                                                                   tgttCCGCGGGGGAGAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaa           <  1:186227/90‑1 (MQ=255)
                                                                           cGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGc   <  1:90574/90‑1 (MQ=255)
                                                                            ggggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct  >  1:142203/1‑90 (MQ=255)
                                                                                       |                                                                              
GTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCT  >  NZ_CP009273/2621897‑2622062

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