Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,213,8700AC66.7% 11.9 / 19.9 21T363P (ACC→CCC) pdeGcyclic diguanylate phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (7/0);  new base C (0/14);  total (7/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 = 8.84e-05
Rejected as consensus: Frequency below/above cutoff threshold.

TTTCATATTGGAATTAATTTTAGCGCTTCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAATTTCCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACA  >  NZ_CP009273/1213786‑1213954
                                                                                      |                                                                                    
tttCATATTGGAATTAATTTTAGCGCTTCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAATTTCCGTGACAGTTTTACccgc                                                                                 >  2:641499/1‑90 (MQ=255)
        tGGAATTAATTTTAGCGCTTCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAATTTCCGTGACAGTTTTACCCGCCGCGAttt                                                                         >  1:818085/1‑90 (MQ=255)
         gaaattatttTTAGCGCTTCGCATATTATTCCGCCGACTTTTGTCGACGAGTTTTTAAATTTCCGTGACATTTTTCCCCGCCGCGATTTa                                                                        <  2:99659/88‑1 (MQ=255)
         gaaattatttTTAGCGCTTCGATTTTTTTTCCGCCGATTTTTGTCGACGTGTTTTTAATTTTCCGTGAAATTTTTCCCCGCCGCGATTTa                                                                        <  1:866854/88‑1 (MQ=255)
              taattTTAGCGCTTCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAATTTCCGTGACAGTTTTACCCGCCGCGATTTAAACCt                                                                   >  2:3789/1‑90 (MQ=255)
               atttttAGCGCTTCGCATATTTTTTCGCCGCCGTTTGTCGAGGATTTTTTAAATTCCCGTGACATTTTTCCCCGCCGCGATTTAAACCtt                                                                  <  1:676936/88‑1 (MQ=255)
                   ttAGCGCTTCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAATTTCCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTc                                                              >  1:263397/1‑90 (MQ=255)
                           tCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAATTTCCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTc                                                      >  1:337335/1‑90 (MQ=255)
                                   tttttCGCCGACGTTTGTCGACGAGTTTTTAAATTTCCGTGACAGTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCg                                              <  2:241905/88‑1 (MQ=255)
                                       ccGCCGACGTTTGTCGAGGAGTTTTAAATTTCCCGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAg                                          <  2:266137/89‑1 (MQ=255)
                                             atgTTTGTGGAGGGGTTTTTAATTTTCCGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTg                                    <  1:59741/88‑1 (MQ=255)
                                             aCGTTTGTCGACGATTTTTTAATTTTCCGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTg                                    <  1:281588/90‑1 (MQ=255)
                                                 ttGTCGACGAGTTTTTAATTTTCCGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATg                                <  2:797985/90‑1 (MQ=255)
                                                   gTCGACGAGTGTTAAAATTTCGGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGtt                              <  2:677100/90‑1 (MQ=255)
                                                              tttAAATTTTCCGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTc                    <  2:486046/90‑1 (MQ=255)
                                                               ttAAATTTTCGGTGAAATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCt                   <  2:589878/90‑1 (MQ=255)
                                                                      ttCCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAg            >  1:589878/1‑90 (MQ=255)
                                                                       tCCGTAACTTTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGc           <  1:326571/90‑1 (MQ=255)
                                                                        cccGTGACTTTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGc           <  1:343123/89‑1 (MQ=255)
                                                                            tGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGtt       <  2:818085/90‑1 (MQ=255)
                                                                                 gTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACa  >  2:92976/1‑90 (MQ=255)
                                                                                      |                                                                                    
TTTCATATTGGAATTAATTTTAGCGCTTCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAATTTCCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACA  >  NZ_CP009273/1213786‑1213954

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