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

GCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTT  >  NZ_CP009273/2621895‑2622072
                                                                                          |                                                                                        
gCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGa                                                                                          >  2:639237/1‑90 (MQ=255)
gCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGa                                                                                          >  2:488516/1‑90 (MQ=255)
          aaaCACAATCTGCCAGAATATCCGCGGGAAAATTAAACGGGTACCCCCATGTGAAACTTTTTTACGCCGGGGAGTAAAGCCCCTTTAAcc                                                                                <  1:488516/90‑1 (MQ=255)
                    tGCCAAAATATACGCGGGGAAATTAAACGGTTTCCCCCTTGAAAACTTTTGTTCCCCCGGGGAGTAAAGCCCCTTTAACCGCCCCAACag                                                                      <  1:217877/90‑1 (MQ=255)
                    cgCCAGAATAACCGCGGGAAATTTAACGGGTTATCACCTTGGGAACCGTTGTTACGCGGGGGAGAAAAGCCCTTTTAACCGCCCCAACag                                                                      <  2:189881/89‑1 (MQ=255)
                    agacAGAATATACGGTGGGAAATTAAACGGGTATCCCCATGTGAAACGTTTTTCCGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAACag                                                                      <  2:693929/87‑1 (MQ=255)
                      ccAGAATATCCGCGGGGAAATTAAACGGGTCTCCCCATGTGAAACGTTGTTGCGCCGGGGGGAAAAGACCCTGTAACCGCCCCAACagag                                                                    <  2:636518/90‑1 (MQ=255)
                       cAGAATATACGCTGGAAAATTAAACGGGTATCCCCATGTGAAACGTTGTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAgaga                                                                   <  2:661285/90‑1 (MQ=255)
                                    gggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAAtgataatgataa                                                      >  1:349830/1‑90 (MQ=255)
                                              aggggTACCACCATGAAAAAATTTGTTAGGGCGGGGAAAAAAGCCCCTATAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGc                                            <  1:310631/88‑1 (MQ=255)
                                              aCGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGc                                            >  1:714300/1‑90 (MQ=255)
                                                            tgtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGc                               >  1:728647/1‑90 (MQ=255)
                                                                     ttgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaa                      >  1:19482/1‑90 (MQ=255)
                                                                     ttgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaa                      >  2:253784/1‑90 (MQ=255)
                                                                      ttttACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaa                     <  2:523055/88‑1 (MQ=255)
                                                                           cgcgGGGGAAAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAccc                <  1:203162/90‑1 (MQ=255)
                                                                           acgggggaAGTAAAGCCTCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAccc                <  2:19147/86‑1 (MQ=255)
                                                                              gggggAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGc             <  2:328975/89‑1 (MQ=255)
                                                                              gggggAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGc             <  2:504427/89‑1 (MQ=255)
                                                                              ggggaAGAAAAGCCCTTGAACCCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGc             <  1:610261/89‑1 (MQ=255)
                                                                                gggAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGctc           <  2:685766/90‑1 (MQ=255)
                                                                                gggAGAAAAGCCCTTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGctc           <  2:472468/90‑1 (MQ=255)
                                                                                         gACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCatctgt  >  1:19147/1‑86 (MQ=255)
                                                                                          |                                                                                        
GCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTT  >  NZ_CP009273/2621895‑2622072

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