Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,213,8650GT57.1% 9.0 / 30.2 21S361I (AGT→ATT) pdeGcyclic diguanylate phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (8/1);  new base T (0/12);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 4.42e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.95e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTAGCGCTTCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAATTTCCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTT  >  NZ_CP009273/1213805‑1213949
                                                            |                                                                                    
ttAGCGCTTCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAATTTCCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTc                                                         >  1:651888/1‑90 (MQ=255)
     gcTTCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAATTTCCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGaa                                                    >  2:337149/1‑90 (MQ=255)
      cTTCGTATATTATTTCGCCGACTTTTGTCGACGGGTTTTTAAATTTCCGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAg                                                   <  1:759770/90‑1 (MQ=255)
       ttCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAATTTCCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGt                                                  >  1:492571/1‑90 (MQ=255)
       ttCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAATTTCCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGt                                                  >  1:92701/1‑90 (MQ=255)
            atatTTTTTCGCCGACTTTTGTCGAGGAGTTTTTAATTTTCCGTGAAGGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCg                                             <  2:418921/90‑1 (MQ=255)
                tattTCGCCGACGTTTGTCGACGAGTGTTTAAATTTCCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCg                                         >  2:1132085/1‑90 (MQ=255)
                       ccGACGTTTGTCGACGAGTGTTTAAATTTCCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCa                                  >  2:757109/1‑90 (MQ=255)
                                          gtTTAAATTTCCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTc               >  1:400309/1‑90 (MQ=255)
                                          gtTTAAATTTCCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTc               >  1:716057/1‑90 (MQ=255)
                                           tttAATTTTCCGTGACATTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCt              <  1:122504/90‑1 (MQ=255)
                                                 tttCCGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCa        <  1:989637/90‑1 (MQ=255)
                                                 ttcccGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCa        <  1:300376/87‑1 (MQ=255)
                                                 ttcccGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCa        <  2:465789/87‑1 (MQ=255)
                                                   tCCGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGc      <  2:940445/90‑1 (MQ=255)
                                                   tCCGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGc      <  2:712874/90‑1 (MQ=255)
                                                     cGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCgg    <  2:445744/90‑1 (MQ=255)
                                                     cGTGACATTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCgg    <  2:796034/90‑1 (MQ=255)
                                                       tcaCATTTTCACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGtt  <  2:327343/88‑1 (MQ=255)
                                                       tGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGtt  <  2:354889/90‑1 (MQ=255)
                                                       tGACATTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGtt  <  2:244449/90‑1 (MQ=255)
                                                            |                                                                                    
TTAGCGCTTCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAATTTCCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTT  >  NZ_CP009273/1213805‑1213949

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