Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273353,5910TC56.2% 11.3 / 18.6 16K271E (AAA→GAA) cynRtranscriptional regulator CynR
Reads supporting (aligned to +/- strand):  ref base T (7/0);  new base C (0/9);  total (7/9)
Fisher's exact test for biased strand distribution p-value = 8.74e-05
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:163481/1‑90 (MQ=255)
cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                    >  2:542693/1‑90 (MQ=255)
cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                    >  2:433046/1‑90 (MQ=255)
 tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc                                                                                   >  1:208315/1‑90 (MQ=255)
          tCATTCCCGCCAACACCCGCGCTTTTACCCAACCCCATGTGAAAAAATGCCTCCGCGGCGGCTGTCGGCCAGCTTTTTTCCCGCCGcaaa                                                                          <  1:147389/90‑2 (MQ=17)
             tttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                       >  1:336128/1‑90 (MQ=255)
             tttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                       >  2:167140/1‑90 (MQ=255)
                     aacaacCGCGCATTTATCCAACGCCATGTTCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                               >  1:72552/1‑90 (MQ=255)
                       cacccGCGCTTTTACCCACCCCCTTGTAAAAAAATCCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                             <  1:139512/87‑1 (MQ=255)
                       cacccGCGCTTTTACCCACCCCCATGTGAAAAAATCCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTtc                                                             <  2:496974/87‑1 (MQ=21)
                            gcgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                        >  1:388251/1‑90 (MQ=255)
                                        aaCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                            >  1:349976/1‑90 (MQ=255)
                                         aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                           >  1:12551/1‑90 (MQ=255)
                                              ttgtgAAAAAATGCCTTCGGGGGGGCTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                      <  2:242092/89‑1 (MQ=255)
                                                   aaaaaaTCCTTTCGGGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                 <  2:208315/89‑1 (MQ=255)
                                                   aaaaaaTCCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                 <  2:15629/89‑1 (MQ=38)
                                                    aaaaaTGCTTTCGCGGCGGCTGTCTCCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                                <  2:23718/90‑1 (MQ=255)
                                                            ctcggcggggGCTGTCTGCCATTTTTTTTCCCGCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg                        <  1:505460/86‑1 (MQ=255)
                                                                          tctcCCAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtg          <  2:284651/90‑1 (MQ=255)
                                                                                  tttttttcCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc  <  2:252155/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 < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: