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

GCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGC  >  NZ_CP009273/2621895‑2622061
                                                                                         |                                                                             
gCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGa                                                                               <  2:57590/90‑1 (MQ=255)
        gaaaaCACAATCTGCCAGAATATCCGCGGGAAAATAAAACGGGTACCCCCATGTAAAACTTTGTTCCGCGGGGGAGAAAAGCCCCTGaaa                                                                       <  2:673164/88‑1 (MQ=255)
            acacAATCTGCCAGAATATCCGCGGGGAAATTAACGGGGTATCCCCATGTGAAACGTTTTTACGCGGGGGAGTAAAGCCCCTGTAACCGc                                                                   <  1:101563/90‑1 (MQ=255)
               aaaTCTGCCAGAATATCCGCGGGAAAATTAAACGGGTACCCCCATGTAAAACGTTTTTACGCCGGGGAGAAAAGCCCCTGTAACCGcccc                                                                <  1:299018/89‑1 (MQ=255)
                 aTCTGCCTGAATATACGCTGGGAATTTAAACGGGTATTACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCaa                                                              >  2:19890/1‑90 (MQ=255)
                  tCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAAc                                                             <  1:387410/90‑1 (MQ=255)
                  tCTGCAAAAATATCCGCGGGAAAATTAAACGGGTATCACCATGTGAAACGTTGTTCCGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAAc                                                             <  1:532722/90‑1 (MQ=255)
                    tGCCAGAATACCCGCCGGGAAATTACACGGGTCCCCCCATGTTAACCGTTGTTCCCCCGGGGAGTAAAGCCCCTGTAACCGCCCCAACag                                                           <  2:618684/90‑1 (MQ=255)
                      ccAAAATATCCGCTGGAAAATAAAACGGGTATCCCCATGTGAAACGTTGTTCCGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAACagag                                                         <  1:339002/90‑1 (MQ=255)
                       cAGAAAATACGCTGGAAAATTAACCGGGTACCCCCATGTGAAAGGTTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAACAgaga                                                        <  2:206583/90‑1 (MQ=255)
                         gAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaat                                                      >  1:365915/1‑90 (MQ=255)
                                cGCTGGGAAATAAAACGGATAACCCCTTGAGAAAGGTTTTTACGCGGGGGAGAAAAGCCCCTTTAACCGCCCCAACAGAGaatgataatg                                               <  2:502702/90‑1 (MQ=255)
                                   tGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAatgataatgata                                            >  1:762910/1‑90 (MQ=255)
                                   tGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAatgataatgata                                            >  2:119540/1‑90 (MQ=255)
                                     ggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATgataatgataat                                          >  1:561186/1‑90 (MQ=255)
                                              ccGGGTATCACAATGTAAAACTTTTTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGc                                 <  1:559539/89‑1 (MQ=255)
                                                  gTATCACCATGTAAAACGTTGTTACGCGGGGGAGAAAAGCCCCTTTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGtt                             <  2:561186/90‑1 (MQ=255)
                                                              aaaaGGTTGTTACGCGGGGAAAAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGtt                 <  2:609278/89‑1 (MQ=255)
                                                                aaCTTTTTTACGCGGGGGAGAAAAGCCCCTTTACCCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTc               <  2:365915/90‑1 (MQ=255)
                                                                  cGTTGTTACGCCGGGGAGTAAAGGCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGt             >  1:697098/1‑90 (MQ=255)
                                                                    ttgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaa           >  2:552987/1‑90 (MQ=255)
                                                                      gttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaaa         >  2:616224/1‑90 (MQ=255)
                                                                             cGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGc  <  1:616224/90‑1 (MQ=255)
                                                                                         |                                                                             
GCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGC  >  NZ_CP009273/2621895‑2622061

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