Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273353,6010CA64.5% ‑0.5 / 20.5 31L267F (TTG→TTTcynRtranscriptional regulator CynR
Reads supporting (aligned to +/- strand):  ref base C (9/2);  new base A (0/20);  total (9/22)
Fisher's exact test for biased strand distribution p-value = 2.73e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.18e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTA  >  NZ_CP009273/353520‑353683
                                                                                 |                                                                                   
gCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                            >  1:637342/1‑90 (MQ=255)
     cacccGCTCTTTTCTCCACCCCCATGAAAAAAAAGGCCTCCGGGGGGGTTGTCGGCAAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTtc                                                                       <  2:951807/87‑1 (MQ=255)
     cacccGCGCTTTTACCCACCCCCTTGTAAAAAAAGGCTTTCGGGGGGCTTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTtc                                                                       <  1:118683/87‑1 (MQ=255)
     cacccGCGCTTTTACCCACCCCCATGTCAAAAAAGCCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                       <  1:102477/87‑1 (MQ=19)
     cacccGCGCTTTTACCCAACCCCATGTGAAAAAATCCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTtc                                                                       <  1:151883/87‑1 (MQ=21)
     cacccGCCCTTTTACCCACCCCCTTGTAAAAAAATCCCTCCGGGGGGGTTTTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                       <  2:280635/87‑1 (MQ=255)
        ccGCGCTTTTACCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTtctct                                                                    <  1:90356/90‑1 (MQ=21)
          gcgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                  >  1:794764/1‑90 (MQ=255)
          ccctaTTTACCCAACCCCATGTAAAAAAATCCTTTCGCGGCGGTTGTCGGCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTcc                                                                  <  1:697181/86‑1 (MQ=255)
           agaaTTTACCCACCCCCATGTAAAAAAATCCCTCCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCa                                                                 <  1:1102463/87‑1 (MQ=21)
                       accccATGTGCAAAAATGCCTTCGGGGGGGCTTTCTGCCAGTTTTTTTTCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                     <  1:41075/87‑1 (MQ=37)
                       aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggt                                                     >  2:1102463/1‑89 (MQ=255)
                       aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                     >  2:933182/1‑90 (MQ=255)
                       aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                     >  2:765090/1‑90 (MQ=255)
                            aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCCGTGGTGGCGGGGccgga                                                >  1:676079/1‑86 (MQ=37)
                            aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGagag                                               >  2:480905/1‑90 (MQ=37)
                            aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGagag                                               >  2:379426/1‑90 (MQ=37)
                                  aaaatTCCCCTCGCGGGGGTTGTCGCCCATTTTTTTTCCCCCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                                          <  2:1131791/90‑1 (MQ=255)
                                  aaaaaTCCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                                          <  1:933182/90‑1 (MQ=38)
                                  aaaaaCCCCTCCGCGGGGGTTCTCGGCCGGTTTTTTTCCCCCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                                          <  1:648260/90‑1 (MQ=255)
                                            tCGGGGCGGCTGTCTGCCAGCTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                <  1:299548/90‑1 (MQ=255)
                                            tCGCGGGGGCTGTCTGCCGTTTTTTTTCCCCCGCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                <  1:1130163/90‑1 (MQ=255)
                                            ccggggcggcTTTCTCCCGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                <  2:155858/86‑1 (MQ=255)
                                               gggcggtTGTCGGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt                             <  1:226742/89‑1 (MQ=255)
                                                 gcggcTGTCTCCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCa                           <  1:1101121/90‑1 (MQ=255)
                                                     ttGTCGGCCACTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgtt                       <  1:527248/89‑1 (MQ=255)
                                                                  tttttcCCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc          <  1:379426/89‑1 (MQ=255)
                                                                  tttttcCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc          <  1:485822/89‑1 (MQ=255)
                                                                  tttttcCCCCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc          <  1:183814/89‑1 (MQ=255)
                                                                   tttttCCCCCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc          <  2:472067/90‑1 (MQ=255)
                                                                          gccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAAACCCCCCTGTTGTGTGGGGAACGCTGCTgggg  >  1:1130400/1‑88 (MQ=255)
                                                                          cccacaaaaaACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAACCCGTCATGTTGTGTGGCAATCGCTGCTGGTa  <  1:765090/86‑1 (MQ=255)
                                                                                 |                                                                                   
GCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTA  >  NZ_CP009273/353520‑353683

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