Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273353,5860AT41.2% 36.8 / 18.6 34N272K (AAT→AAAcynRtranscriptional regulator CynR
Reads supporting (aligned to +/- strand):  ref base A (15/5);  new base T (0/14);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 9.80e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAT  >  NZ_CP009273/353497‑353659
                                                                                         |                                                                         
aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa                                                                           >  2:213152/1‑90 (MQ=255)
aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa                                                                           >  1:259599/1‑90 (MQ=255)
aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa                                                                           >  2:339651/1‑90 (MQ=255)
aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa                                                                           >  2:36606/1‑90 (MQ=255)
aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa                                                                           >  2:195241/1‑90 (MQ=255)
 gCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAt                                                                          >  2:332501/1‑90 (MQ=255)
     cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                      >  2:102818/1‑90 (MQ=255)
     cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                      >  1:287325/1‑90 (MQ=255)
     cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                      >  1:388346/1‑90 (MQ=255)
     cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                      >  1:88118/1‑90 (MQ=255)
      tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc                                                                     >  1:193228/1‑90 (MQ=255)
          gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                 <  1:132722/90‑1 (MQ=255)
          gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                 >  2:10200/1‑90 (MQ=255)
          gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                 <  1:449017/90‑1 (MQ=255)
          gTGATTCATTTCCGCCAACAACCCCCCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                 <  1:95447/90‑1 (MQ=255)
            gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                                               >  2:327135/1‑90 (MQ=255)
             attcattcCCGCCAACAACCGCCCTTTTACCCAACGCCATGGGCAAAAAGGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGca                                                              <  2:341104/90‑1 (MQ=14)
               tCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                            >  1:199927/1‑90 (MQ=255)
                 aTTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                          >  2:21503/1‑90 (MQ=255)
                          aacaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                 >  2:57034/1‑90 (MQ=255)
                          aacaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                 >  2:35415/1‑90 (MQ=255)
                            cacccGCGCTTTTACCCACCCCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                               <  1:130638/87‑1 (MQ=21)
                            cacccGCGCTTTTACCCAACGCCTTGTAAAAAAATGCCTTCGCGGGGGCTGTCTGCCAGTTTTTTTTCCGCCCAAAAAAAACCGCCGTtc                                               <  1:439003/87‑1 (MQ=16)
                            cacccGCGCTTTTACCCAACGCCATGTACAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                               <  1:327135/87‑1 (MQ=21)
                            cacccGCGCTTTTACCCAACGCCAGGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTtc                                               <  1:52638/87‑1 (MQ=25)
                            cacccGCGATTTTACCCAACCCCTTGTGAAAAAATCCTTTCGGGGCGGTTGTCTGCCACTTTTTTTTCCGCCGCAAAAAAACCGCCGTtc                                               <  1:36606/87‑1 (MQ=255)
                              cccGCGCTTTTACCAACCCCCTTGTAAAAAAATCCTTTCGGGGCGGCTGTCGGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTtctc                                             <  1:102818/89‑1 (MQ=255)
                               ccGCGCATTTCCCCAACCCAAGGTAAAAAAATGCCTTCGGGGCGGCTGTCTCCCATCTTTTTTCCCGCCGAAAAAAAACCGCCGTtctct                                            <  1:458444/90‑1 (MQ=16)
                                 gcccATTTACCCAACGCCATGTGAAAAAATGCCTCCGCGGCGGCTGCCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTcc                                          <  1:7075/87‑1 (MQ=255)
                                              aCGCCATGTGCTAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                             >  2:329467/1‑90 (MQ=255)
                                              aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                             >  1:163562/1‑90 (MQ=255)
                                                cccATGTGAAAAAATCCCTCCGCGGCGGCTGTCCGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                           <  1:213152/89‑1 (MQ=255)
                                                   aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                                              <  2:308270/52‑1 (MQ=255)
                                                   aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                                              >  1:308270/1‑52 (MQ=255)
                                                               gCCTCCGGGGCGGCTGTCCCCCAGTTTTTTTCCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTa            <  1:324677/90‑1 (MQ=255)
                                                                   tCGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc        <  1:157136/90‑1 (MQ=255)
                                                                      ggggggcTTTCGGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt     <  1:216078/86‑1 (MQ=255)
                                                                        ggggcTGTCTGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCa   <  1:112372/88‑1 (MQ=255)
                                                                         gggcTTTCTGCCATTTTTTTTCCCCCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAt  <  2:43288/89‑1 (MQ=255)
                                                                                         |                                                                         
AGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAT  >  NZ_CP009273/353497‑353659

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: