Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,213,8650GT65.0% ‑3.6 / 21.2 20S361I (AGT→ATT) pdeGcyclic diguanylate phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (6/1);  new base T (0/13);  total (6/14)
Fisher's exact test for biased strand distribution p-value = 1.81e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.25e-03
Rejected as consensus: E-value score below prediction cutoff.
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: