Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,544,9230CA52.9% 9.2 / 25.2 17N177K (AAC→AAAuxuRUxu operon transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (1/7);  new base A (9/0);  total (10/7)
Fisher's exact test for biased strand distribution p-value = 4.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGCG  >  NZ_CP009273/4544842‑4544990
                                                                                 |                                                                   
gTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGGAAAAAAACCAAt                                                             >  2:78235/1‑90 (MQ=255)
 ttCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAAtg                                                            <  1:106802/90‑1 (MQ=255)
  tCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGGAAAAAAACCAAagt                                                           >  2:106802/1‑90 (MQ=255)
     aTCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAAACCAATTTggg                                                        >  2:470006/1‑89 (MQ=255)
                     aGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGGGAAAAAACCAATTAGGGGCCAGGTTCCCAGCCa                                        >  1:458655/1‑90 (MQ=255)
                       agaCGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATc                                      <  2:198042/88‑1 (MQ=255)
                                   gCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGGAAAAAAAACAAAGTGGGGCCCGTTTTACAGCCCTCTTGGTGACaact                          >  1:198042/1‑87 (MQ=255)
                                      tGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAAAAAACCAAAGTGGATCCCGTTGCACAACCCTCTGGATGACAGCCCGt                       >  2:40892/1‑90 (MQ=255)
                                             ggAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCaa                <  1:377382/90‑1 (MQ=255)
                                               aGCTGTTCCGTCAGTCCTGGCAGTGGCGGGGAAAAAAACCAATGTGGGTCCAGTTGCACAGCCCTCCGGGTGGCAGCCTGTATCGCaaaa              >  1:470276/1‑89 (MQ=255)
                                                    ttCCGTCAGTCCTGGCAGTGGCGGGGAAAAAAACCAAAGTGGGTCCCGTTGCACAGCCCTCTGGGTGACAGCCCGTATCGCCAAGAgggg         >  1:185444/1‑90 (MQ=255)
                                                        gTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGcc                          >  2:221554/1‑69 (MQ=255)
                                                        gTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGcc                          <  1:221554/69‑1 (MQ=255)
                                                          cAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGc   <  2:185444/90‑1 (MQ=255)
                                                          cAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGc   <  1:470006/90‑1 (MQ=255)
                                                          cAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGc   <  1:40892/90‑1 (MQ=255)
                                                           aGTCCTGGCAGTGGCGGGAAAAAAAAACAAAGTGGATCCCGTTGCCCAGCCATCTGGGTGACAACCCGTATTGCAAAGAGGGGGTggggg  >  1:321/1‑88 (MQ=255)
                                                                                 |                                                                   
GTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGAAAACAATCCAATGTGGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGCG  >  NZ_CP009273/4544842‑4544990

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