Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273353,5910TC56.2% 10.4 / 19.6 16K271E (AAA→GAA) cynRtranscriptional regulator CynR
Reads supporting (aligned to +/- strand):  ref base T (6/1);  new base C (0/9);  total (6/10)
Fisher's exact test for biased strand distribution p-value = 8.74e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-04
Rejected as consensus: Frequency below/above cutoff threshold.

CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATC  >  NZ_CP009273/353502‑353673
                                                                                         |                                                                                  
cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                    >  2:293897/1‑90 (MQ=255)
     gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                               <  1:509304/90‑1 (MQ=255)
       gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                                                             >  1:436368/1‑90 (MQ=255)
         ttCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaa                                                                           >  2:444618/1‑90 (MQ=255)
                       cacccGCGTTTTTACCCACCGCCATGCAAAAAAACCCCTCCGGGGGGGCTGTCTGCCATTTTTTTTCCCCCCGAAAAAAAACCGCCGTtc                                                             <  2:513055/87‑1 (MQ=255)
                       cacccGCGCATTTACCCAACCCCATGTGAAAAAATGCTTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTtc                                                             <  1:231286/87‑1 (MQ=21)
                       aacccGCGCTTTTACCCAACCCCTTGTGAAAAAATGCCTCCGGGGGGGCTGTCGGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTtc                                                             <  1:351095/87‑1 (MQ=255)
                            gcgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                        >  1:11634/1‑90 (MQ=255)
                                         aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                           >  1:278085/1‑90 (MQ=255)
                                         aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                           >  1:405481/1‑90 (MQ=255)
                                         aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                           >  2:345405/1‑90 (MQ=255)
                                                   aaaaaaGCCCTCCGGGGGGGCTGTCCCCCATCTATTTTCCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                 <  1:418023/89‑1 (MQ=255)
                                                         gcccTCCGCGGGGGCTGCCGCCCTGCTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt                           <  2:84363/88‑1 (MQ=255)
                                                            cTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg                        <  1:90769/90‑1 (MQ=255)
                                                              tCGGGGGGGCTTTCTCCCAGTTTTTTTCCCGCCCCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                      <  1:549470/90‑1 (MQ=255)
                                                              tCGGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                      <  1:30744/90‑1 (MQ=255)
                                                              gcgcgGGGCCTGTCTGCGATTTTTTTCCCCGCCACAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                      <  2:622408/89‑1 (MQ=255)
                                                                                  tttttttcCCGCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc  <  2:499793/87‑1 (MQ=255)
                                                                                         |                                                                                  
CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATC  >  NZ_CP009273/353502‑353673

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