Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273353,6010CA54.3% 2.0 / 38.2 35L267F (TTG→TTTcynRtranscriptional regulator CynR
Reads supporting (aligned to +/- strand):  ref base C (10/6);  new base A (0/19);  total (10/25)
Fisher's exact test for biased strand distribution p-value = 4.36e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTA  >  NZ_CP009273/353516‑353683
                                                                                       |                                                                                    
ttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                  >  1:26133/1‑90 (MQ=255)
  ccGCAAACACCCGCCCATTTATCCACCCCCATGTGCAAAAACGCTTCCGCGGCGGCTGTCTGCCCGATTTTTTCCCGCCGCAACAAAG‑‑cc                                                                              <  2:221333/90‑1 (MQ=255)
    gCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑ccgc                                                                            >  2:112471/1‑90 (MQ=255)
         cagccGCGCTTTTACCCACCCCCATGTGCAAAAATCCCTTCGGGGGGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAA‑‑CCGCCGTtc                                                                       <  2:789221/87‑1 (MQ=21)
         cacccGCGTTTTTACCCACCCCCATGTGAAAAAATCCCTTCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑‑CCGCCGTtc                                                                       <  2:364560/87‑1 (MQ=16)
         cacccGCGTTTTTACCCACCCCCATGTAAAAAAATGCCTTCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑‑CCGCCGTtc                                                                       <  2:751823/87‑1 (MQ=16)
         cacccGCGTTTTTACCAACCGCCATGTGAAAAATTCCCTTCGGGGGGGCTGTCTCCCAGTTTTTTTCCCGCCGCAACAAAA‑‑CCGCCGTtc                                                                       <  2:361212/87‑1 (MQ=255)
         cacccGCGCTTTTATCCAACGCCATGTGAAAAAATGCCTTCGCGGGGGCTGTCGGCCAGTTTTTTTCCCGCCGAAACAAAA‑‑CCGCCGTtc                                                                       <  2:706271/87‑1 (MQ=21)
         caacCGCGCATTTACCCAACCCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAA‑‑CCGCCGTtc                                                                       <  2:136165/90‑1 (MQ=21)
            ccGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTtctct                                                                    >  1:849809/1‑90 (MQ=255)
            ccGCGATTTTACCCACCCCCATGTGAAAAAAGCCCTTCGGGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑‑CCGCCGTtctct                                                                    <  2:402763/90‑1 (MQ=255)
              gcgcTTTTACCCACCGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAA‑‑CCGCCGTTCTCTcc                                                                  <  2:26133/90‑1 (MQ=25)
              gcgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTcc                                                                  >  2:317479/1‑90 (MQ=255)
                      aTCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTgg                                                          >  1:79586/1‑90 (MQ=255)
                           ccaccATGTCAAAAAATGCCTCCGGGGCGGTTTCCTCCCGGCTTTTTTCCCCCCAAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCgggg                                                     <  2:493122/87‑1 (MQ=255)
                           aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGcctgt                                                     >  1:584720/1‑86 (MQ=255)
                           aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCgggg                                                     >  1:622331/1‑90 (MQ=255)
                           aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCgggg                                                     >  2:271343/1‑90 (MQ=255)
                                aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAaga                                                >  2:880831/1‑90 (MQ=255)
                                aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCCaga                                                >  2:823344/1‑90 (MQ=255)
                                     cAAAAAGCCCTCCGGGGCGGCTGTCGCCCATTTTTTTTCCCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                           <  1:16884/90‑1 (MQ=255)
                                      aaaaaTGCTTTCGGGGGGGTTGTCGGCCGGTTTTTTTCCCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                                          <  2:441221/90‑1 (MQ=255)
                                      aaaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                                          <  2:53291/90‑1 (MQ=38)
                                      aaaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGGCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                                          <  1:437405/90‑1 (MQ=255)
                                      aaaaaTCCCTTCCGCGGCGCCTGTCGGCCATTTTTTTTCCCCCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                           <  2:625053/90‑1 (MQ=255)
                                      aaaaaTCCCCTCGGGGGGGGCTGTCGGCCGTTTTTTTTCCCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                           <  1:14392/90‑1 (MQ=255)
                                                tCGCGGCGGCTGTCTGCCAGTTTTTTTCCCCCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                <  2:79586/90‑1 (MQ=255)
                                                 cggcggcgGCTTTCTCCCTTTTTTTTCCCCGCCAAAAAAAAC‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                <  1:833336/88‑1 (MQ=255)
                                                          gcTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                        <  2:100965/90‑1 (MQ=255)
                                                                        ttttttCCCCCCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc            <  1:872566/89‑1 (MQ=255)
                                                                         tttttCCCCCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc          <  1:391665/90‑1 (MQ=255)
                                                                          tttcCCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg        <  2:300595/90‑1 (MQ=255)
                                                                          tttcCCGCCCCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg        <  2:811090/90‑1 (MQ=255)
                                                                          ttccccgcccAAAAAAAC‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg        <  1:763435/86‑1 (MQ=255)
                                                                                cccgcAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa  <  2:812190/89‑1 (MQ=255)
                                                                                       |                                                                                    
TTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTA  >  NZ_CP009273/353516‑353683

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