Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273353,6010CA56.7% 4.7 / 32.1 30L267F (TTG→TTTcynRtranscriptional regulator CynR
Reads supporting (aligned to +/- strand):  ref base C (10/3);  new base A (0/17);  total (10/20)
Fisher's exact test for biased strand distribution p-value = 9.52e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAA  >  NZ_CP009273/353523‑353687
                                                                                |                                                                                      
aacaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                             >  2:514924/1‑90 (MQ=255)
  cccccGCGTTTTTATCCACCCCCTTGTCAAAAACGCCCTCCGCGGCGGCTGTCTGCAAGATTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                           <  1:682420/87‑1 (MQ=255)
  cacccGCGTTTTTACCCAACCCCATGGGAAAAAATCCTTTCGGGGGGGCTTTCGGCCGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                           <  1:23364/87‑1 (MQ=255)
  cacccGCGCTTTTACCCACCGCCTTGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                           <  2:662704/87‑1 (MQ=21)
  cacccGCGCTTTTACCCACCGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                           <  2:32330/87‑1 (MQ=21)
  cacccGCGCATTTCCCCAACCCCATGTGAAAAAATCCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                           <  2:536579/87‑1 (MQ=16)
  cacccGCGCATTTACCCACCCCCATGTGAAAAAATCCCTTCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTtc                                                                           <  2:281065/87‑1 (MQ=21)
      cgcgCTTTTACCCACCCCCATGGAAAAAAATCCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAACCCGCCGTtctctc                                                                       <  1:108208/90‑1 (MQ=255)
       gcgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                      >  2:628305/1‑90 (MQ=255)
       gcgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                      >  1:297003/1‑90 (MQ=255)
              tATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                               >  1:288609/1‑90 (MQ=255)
                    aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                         >  2:684425/1‑90 (MQ=255)
                    aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                         >  1:408244/1‑90 (MQ=255)
                    aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                         >  2:451181/1‑90 (MQ=255)
                    aCGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                         >  1:808956/1‑90 (MQ=255)
                      cccttgtgAAAAAATGCGTCCGCGGCGGCTGTCTGCCAGTTTTTTTTCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                       <  2:297003/86‑1 (MQ=255)
                         aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGcccgg                                                    >  1:330967/1‑86 (MQ=255)
                              aaaaaaCGCCTCGGGGGCGGCTGTCTGCCGGTTTTTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                               <  1:574974/89‑1 (MQ=255)
                               aaaaaTTCCTTCCGGGGGGGCTGTCGGCCAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                               <  2:288609/90‑1 (MQ=255)
                                aaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAGaaaaaa                                             >  2:374540/1‑87 (MQ=39)
                                          tCGGGGCGGCTGTCTGCCCTTTTTTTTCCCGCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                    <  2:666556/90‑1 (MQ=255)
                                           gacgggcggcTTCCTCCCAGTTTTTTTCCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                    <  2:527761/86‑1 (MQ=255)
                                           ccggggcggcTGTCTGCCACTTTTTTTCCCCCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                    <  2:590147/86‑1 (MQ=255)
                                           ccgcgGGGGCTGTCTGCCACTTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                    <  2:196576/89‑1 (MQ=255)
                                                   gcTGTCGGCCAGTTTTTTTCCCCCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                            <  2:67443/90‑1 (MQ=255)
                                                                 tttttcCCGCCCGAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc              <  2:394515/89‑1 (MQ=255)
                                                                 tttttcCCCCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc              <  2:43887/89‑1 (MQ=255)
                                                                 tttttcCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc              <  2:882148/89‑1 (MQ=255)
                                                                  ttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct             <  2:869798/90‑1 (MQ=255)
                                                                   tttcCCCCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg            <  2:502482/90‑1 (MQ=255)
                                                                          ccacaaaaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa     <  2:383530/87‑1 (MQ=255)
                                                                           cacaaaaaACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc    <  1:221962/88‑1 (MQ=255)
                                                                             caacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa  <  1:684425/90‑1 (MQ=255)
                                                                                |                                                                                      
AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAA  >  NZ_CP009273/353523‑353687

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