Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273353,6010CA58.8% ‑2.5 / 16.5 17L267F (TTG→TTTcynRtranscriptional regulator CynR
Reads supporting (aligned to +/- strand):  ref base C (5/2);  new base A (0/10);  total (5/12)
Fisher's exact test for biased strand distribution p-value = 3.39e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.50e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTA  >  NZ_CP009273/353515‑353683
                                                                                      |                                                                                   
tttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                  >  2:412438/1‑90 (MQ=255)
      ccAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                            >  1:118940/1‑90 (MQ=255)
        aacaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                          >  1:458306/1‑90 (MQ=255)
          cacccGCGCTTTTACCCACCCCCTTGCAAAAAAATCCCTCCGGGGCGGCTGCCCCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                        <  2:398654/87‑1 (MQ=255)
          cacccGCGCTTTTACCCAACGCCATGTGCAAAAATGCCTCCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTtc                                                                        <  1:114954/87‑1 (MQ=21)
                  cATTTACCCAACCCCATGTGAAAAAATGCCTTCGGGGGGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagt                                                                <  1:124226/90‑1 (MQ=25)
                                                 tCGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                <  1:103479/90‑1 (MQ=255)
                                                 ccggggcggtTGTCTGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                <  2:456715/86‑1 (MQ=255)
                                                 ccggggcggcTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                <  1:412438/86‑1 (MQ=255)
                                                    cggcggCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAAGAAAAAACTTAAACCCCg                              >  1:314562/1‑90 (MQ=255)
                                                         gcTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                        <  2:348082/90‑1 (MQ=255)
                                                                       tttttcCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc          <  2:314562/89‑1 (MQ=255)
                                                                       tttttcCCCCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc          <  1:216838/89‑1 (MQ=255)
                                                                       tttttcCCCCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGAAAGAGAAATAGCTCTAAGCCCGTCATGTTGTGTGGCAATCgc          <  1:454840/89‑1 (MQ=255)
                                                                         tttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg        <  1:254569/90‑1 (MQ=255)
                                                                               gccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAACCCCCCCTGTTGTTTGGGGAACGCTGCTgggg  >  1:185511/1‑88 (MQ=255)
                                                                               gccgaAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa  <  1:387067/90‑1 (MQ=255)
                                                                                      |                                                                                   
TTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTA  >  NZ_CP009273/353515‑353683

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