Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273324,5140TG63.2% ‑2.6 / 17.6 19T178P (ACC→CCC) betItranscriptional regulator BetI
Reads supporting (aligned to +/- strand):  ref base T (1/6);  new base G (12/0);  total (13/6)
Fisher's exact test for biased strand distribution p-value = 2.58e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.96e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAG  >  NZ_CP009273/324392‑324650
                                                                                                                           |                                                                                                                                          
gacagggTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTg                                                                                                                            <  1:457276/134‑1 (MQ=255)
gCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTg                                                                                                                            <  1:227408/139‑1 (MQ=255)
                caTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCCGGGGAGTTGGGCGGGGGTTTTTCAGCGGGTTGGCGCTCcgaggccggc                                                                                                            >  2:457276/1‑133 (MQ=255)
                 aTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAATAACGCGGGGTTTTTTCCCCGGGTTTCCGCCCACAGCcgcgcg                                                                                                           >  1:332833/1‑139 (MQ=255)
                       ataAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAATAGCGCGGGGGTTTTTTAGCGGGTTTCCGCCCCGAGCCGCGCgccgcc                                                                                                      >  1:417790/1‑138 (MQ=255)
                          atgctgctGTTCGGCCATGCGGGACTTCTCGTCCCCTCGGTTAATCGGTGTGTAGATGCTGAGTGTTAAAGTGGCGGGTCGGGAAATTAGGGCGGGTTTTATCCAGCGGTTTGCGGCTCGGAGCCGCGCGCAGCCATa                                                                                                  <  2:324824/136‑1 (MQ=255)
                              tgctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGGATTAGCGCGGGGGTTTTTCAGCGGGTTGCCGCCCCCAGACGCGCGCGGCCCTAAACCa                                                                                             >  2:88628/1‑139 (MQ=255)
                                ctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAATAGCGCGGGGGTTTTTCAGCGGTTTGCCGCTCCGAGCCGCGCGCAGCCATAATCCat                                                                                            >  1:86231/1‑138 (MQ=255)
                                             cGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCCGGGGAATATCGCGGGGTTTTTTCACCGGGTTTTCGCCCCCAGCCGCGCGCAGCCCTAATCCACCAAACAACGCGGcc                                                                             >  1:615661/1‑139 (MQ=255)
                                                     cACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGATTAACGCGGGGGTTTTTCCACGGGTTTGCGCCCCGAGCCGCGCGCCGCCATAAATCCTCCAATAACGCGGCCAgggcc                                                                       >  1:394496/1‑134 (MQ=255)
                                                     cACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAGTAGCGGGGGGTTTTTTCAGCCGGTTGGCGCCCCGAGCCGCGCGGCGCCATAAAACATCCATCAAAGCGGCCCGGCCggc                                                                     >  2:688758/1‑137 (MQ=255)
                                                           ccctcGGTTAATCGGTGGTGAGATGCTGAGGGATAAAGTGGCGGGTGAGGGAATTAGCGCGGGTTTTATCAAGCGTTTTGCCGCTCAGGGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCg                                                                <  2:624450/137‑1 (MQ=255)
                                                                     aTCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGGAATAGCGCGGGGTTTTTTCAGCGGGTTGTCGCTCCGAGCCGCGCGCCGCCCTAATCCATCAATCAGCGCGGCCCGGCCGTAGCCCGCCTCCTGTGcc                                                     >  1:114041/1‑139 (MQ=255)
                                                                                 aTGCTGAGTGATAAAGTGGCGGGGGAGGGGATTAAGGGGGGGTTTTTTCAGCGGGTTTGCGCCCACAGGCGCGGCCCGCCACAAACCATCAATAAACCCGGGCCCGGCCGAGCCCCCCTCCTGTTTCCGTTTGCGCggg                                         >  2:443990/1‑138 (MQ=255)
                                                                                    cTGAGTGATAAAGTGGCGGGGCAGGGGATTAACGCGGGGGTTTTCCAACGGGTTTCCGCCCAGAGCCGCGCGCCGCCCTAAACCCTCCATCAGCGCGGCCCGGCCGCAGCCCGCCTCCTGTGCCCGTTTGCGCgggcat                                      >  1:608364/1‑135 (MQ=255)
                                                                                        gTGATAAAGTGGCGGGTCCGGGGATTAGCGCGGGGTTTTTTCCGCGGGTTGCCGCCCAGAGGCGCGCGCAGCCCTAATCCCTCAATCAGCGCGGCCAGGCCGGAGCCCCCCTCCTGTTCCTGTTTCCCCGGCAATTcgc                                  >  1:348351/1‑139 (MQ=255)
                                                                                                gTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAAc                          <  1:411108/139‑1 (MQ=255)
                                                                                                       gTCAGGGAATTAGCGCGGGTTTTATCCAGCGTTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCa                   <  1:772065/139‑1 (MQ=255)
                                                                                                                         ggTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAg  >  2:242808/1‑139 (MQ=255)
                                                                                                                           |                                                                                                                                          
GCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAG  >  NZ_CP009273/324392‑324650

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