Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273353,5860AT53.3% 8.4 / 42.0 30N272K (AAT→AAAcynRtranscriptional regulator CynR
Reads supporting (aligned to +/- strand):  ref base A (13/1);  new base T (0/16);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 1.17e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.56e-06
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAT  >  NZ_CP009273/353497‑353659
                                                                                           |                                                                         
aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa                                                                           >  2:622874/1‑90 (MQ=255)
aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa                                                                           >  1:1131367/1‑90 (MQ=255)
aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa                                                                           >  2:1108002/1‑90 (MQ=255)
aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa                                                                           >  1:278552/1‑90 (MQ=255)
  cGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAtt                                                                         >  1:391929/1‑90 (MQ=255)
     cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                      >  2:429335/1‑90 (MQ=255)
            gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                                               >  2:1444965/1‑90 (MQ=255)
            gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                                               >  1:841433/1‑90 (MQ=255)
                          aacaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                 >  1:1124077/1‑90 (MQ=255)
                            cacccGCGCTTTTATCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTTCCGCCGCAAAAAAACCGCCGTtc                                               <  2:1077356/87‑1 (MQ=21)
                            cacccGCGCTTTTACCCACCCCCTTGTGAAAAATGCCCTTCGCGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTtc                                               <  2:715100/87‑1 (MQ=255)
                               ccGCGCTTTTACCCAACCCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTtctct                                            <  2:1043565/90‑1 (MQ=25)
                               ccGCGATTTTATCCACCCCCTTGTGAAAAAATGCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTtctct                                            <  2:1124077/90‑1 (MQ=25)
                                 gcgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                          >  1:557449/1‑90 (MQ=255)
                                 gcgcATTTACCCACCCCAATGTGAAAAAATGCCTTCGGGGCGCCTGTCTCCCAGCTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTcc                                          <  2:165777/90‑1 (MQ=255)
                                  cgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCa                                         >  1:159314/1‑90 (MQ=255)
                                          cccACCGCAATGTGAAAAAAGCCTTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGc                                 <  2:1131367/89‑1 (MQ=255)
                                          ccaaaCGCCTTGTCAAAAATTTCCTTCGGGGGGCCTCCCCGCCAGTTTTTTTTCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGc                                 <  2:667318/87‑1 (MQ=255)
                                              aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                             >  1:285280/1‑90 (MQ=255)
                                              aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                             >  1:514276/1‑90 (MQ=255)
                                              aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                             >  1:1119834/1‑90 (MQ=255)
                                              aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                             >  1:590008/1‑90 (MQ=255)
                                              aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                             >  1:1137999/1‑90 (MQ=255)
                                                         aaaagTCCCTCCGGGGCGGCTTTCTCCCGTCTTTTTTGCCCCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                  <  2:788039/90‑1 (MQ=255)
                                                         aaaaaTCCCTCCGGGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                  <  2:735710/90‑1 (MQ=255)
                                                               tcccTCCGGGGCGGCTGTCTCCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt             <  2:882600/88‑1 (MQ=255)
                                                               gCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                 >  2:697820/1‑37 (MQ=255)
                                                               gCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                 <  1:697820/37‑1 (MQ=255)
                                                                    tCGGGGGGGCTGTCTCCCAGTTTTTTTCCCCCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc        <  2:841433/90‑1 (MQ=255)
                                                                    tCGCGGCGGCTGTCTGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc        <  2:1331504/90‑1 (MQ=255)
                                                                     ccggggcggtGGTCTGCCAGCTTTTTTTCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc        <  2:1099553/86‑1 (MQ=255)
                                                                     ccgcgGGGGCTTTCTCCCGTTTTTTTTCCCCCCACAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc        <  1:1154375/89‑1 (MQ=255)
                                                                          gcggcTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCa   <  2:1290062/90‑1 (MQ=255)
                                                                           gggcTGTCTGCGGGTTTTTTTCCCCCCGCAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAt  <  2:1209035/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 < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: