Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273328,7020GT35.3% 68.0 / 12.2 34L207F (TTG→TTTpdeLDNA‑binding transcriptional activator/c‑di‑GMP phosphodiesterase PdeL
Reads supporting (aligned to +/- strand):  ref base G (3/19);  new base T (12/0);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 8.30e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGT  >  NZ_CP009273/328617‑328789
                                                                                     |                                                                                       
aGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTGGCAg                                                                                     >  1:166428/1‑90 (MQ=255)
 gACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGGAGc                                                                                    >  2:326486/1‑90 (MQ=255)
  aCTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCg                                                                                   >  1:6472/1‑90 (MQ=255)
   cTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGGAGCgg                                                                                  >  1:250701/1‑90 (MQ=255)
       ggATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGAtt                                                                              <  1:125533/90‑1 (MQ=255)
               tGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATAATGGCATCAACGTCTCGGCGGGGTGTTTTTTTGGAGCGGGGTTTTtaaaaa                                                                      >  1:445747/1‑89 (MQ=255)
                gATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGGTTTTaaaaaaa                                                                     >  2:454814/1‑88 (MQ=255)
                 aTGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAagag                                                                    <  2:234683/90‑1 (MQ=255)
                 aTGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAagag                                                                    <  2:124536/90‑1 (MQ=255)
                     cGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGATTTGGAAAaaaggggc                                                                >  2:461688/1‑86 (MQ=255)
                      ggTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGGGGGATTTGGAAAAGAGGGTCt                                                               >  2:470827/1‑90 (MQ=255)
                           aaCATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAAcc                                                          <  2:59971/90‑1 (MQ=255)
                           aaCATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAAcc                                                          <  1:203549/90‑1 (MQ=255)
                            aCATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCt                                                         <  2:460871/90‑1 (MQ=255)
                             cATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTg                                                        <  2:168773/90‑1 (MQ=255)
                                 tgctgcCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTta                                                    <  2:382600/90‑1 (MQ=255)
                                      ccGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGGTTTTAAAAAAAGGGTCTGAAACTGGGTAATaaa                                               >  2:38651/1‑90 (MQ=255)
                                              tttCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTaa                                       <  2:119293/90‑1 (MQ=255)
                                                     atTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGGTTTTAAAAAAAGGGTCTGAAACTGGGTAATAAATTAGGGAAAGATaaa                                >  2:46813/1‑90 (MQ=255)
                                                        ggCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGGTTTGAAAAAGAGGGGCTGAACCTGGTTAATAAATTTGGGAACGATaaaaac                             >  1:303452/1‑88 (MQ=255)
                                                         gCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCa                            <  2:290463/90‑1 (MQ=255)
                                                             cAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCt                        <  2:363407/90‑1 (MQ=255)
                                                                  tctcGGCGGGTTGTTTTTTTGCAGCGGGGTTTGAAAAAGAGTGTCTGAACCTGGGTAATAAATTAGGGAACGATAAAATCAAGCTGGTtc                   >  1:341493/1‑90 (MQ=255)
                                                                  tctcGGCGGGGTGTTTTTTTGCAGCGGGGTTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGGAAAGATAAAATCAAGCTgggtt                   >  2:207379/1‑87 (MQ=255)
                                                                    tcGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTtctc                 <  1:207379/90‑1 (MQ=255)
                                                                    tcGGCGGGGTGTTTTTTTGCAGCGGGGTTTGAAAAAAAGGGTCTGAAACTGGTTAATAAATTAGGGAAAGATAAAATCAAAGTGGTtttc                 >  2:340364/1‑90 (MQ=255)
                                                                      ggcggGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGa               <  2:6472/90‑1 (MQ=255)
                                                                        cggGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGc             <  2:195477/90‑1 (MQ=255)
                                                                         gggTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTTGGGAACGATAAAATCAAGCTGGTTCTCGAGCt            >  1:351018/1‑90 (MQ=255)
                                                                                   ttGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGt  <  2:407843/90‑1 (MQ=255)
                                                                                   ttGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGt  <  2:18004/90‑1 (MQ=255)
                                                                                   ttGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGt  <  2:175051/90‑1 (MQ=255)
                                                                                   ttGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGt  <  1:46813/90‑1 (MQ=255)
                                                                                   ttGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGt  <  1:110023/90‑1 (MQ=255)
                                                                                     |                                                                                       
AGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGT  >  NZ_CP009273/328617‑328789

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