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

CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGC  >  NZ_CP009273/353502‑353675
                                                                                         |                                                                                     
cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                       >  1:429917/1‑90 (MQ=255)
cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                       >  1:1209045/1‑90 (MQ=255)
cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                       >  2:625582/1‑90 (MQ=255)
     gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGTTGTCTGCCAGCTATTTTTccgcc                                                                                  >  2:329839/1‑90 (MQ=255)
     gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                  <  1:812271/90‑1 (MQ=255)
       gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                                                                >  1:1138327/1‑90 (MQ=255)
                  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:765090/1‑90 (MQ=255)
                                         aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                             >  2:933182/1‑90 (MQ=255)
                                              aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCCGTGGTGGCGGGGccgga                                        >  1:676079/1‑86 (MQ=37)
                                              aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGagag                                       >  2:379426/1‑90 (MQ=37)
                                              aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGagag                                       >  2:480905/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)
                                                                                         |                                                                                     
CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGC  >  NZ_CP009273/353502‑353675

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: