Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273353,6010CA57.1% ‑0.6 / 22.8 21L267F (TTG→TTTcynRtranscriptional regulator CynR
Reads supporting (aligned to +/- strand):  ref base C (8/1);  new base A (0/12);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 4.42e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTC  >  NZ_CP009273/353464‑353717
                                                                                                                                           |                                                                                                                     
cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                       >  2:600053/1‑138 (MQ=255)
     cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                 >  1:539821/1‑139 (MQ=255)
            ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                                                                                                          >  1:123775/1‑139 (MQ=255)
              ccACCCGCAGCCCTGTGTCAGCGGCTCCCGTGATTCTTTCCCGCCAACACCCGCGCATTTACCCAACGCCATGTGAAAAAATGCCTCCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTtctc                                                                                                        <  1:600053/139‑1 (MQ=255)
                  ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctctc                                                                                                      >  1:319008/1‑137 (MQ=255)
                         ccggtgtGAGCGCTACCGGTAATTCATTCCCCCCACCACCCGCCATTTTACCAAACCCCATGTGAAAAAATCCCTCCGGGGGGGTTGTCGGCAGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                             <  2:847416/136‑1 (MQ=255)
                                  gCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                    >  2:157204/1‑139 (MQ=255)
                                  gCGGCTACCGTGATTCATTTCCGCCAACAACCGAGCATTTATCAAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggagac                                                                                     >  1:136667/1‑134 (MQ=255)
                                       tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGaa                                                                               >  1:240432/1‑138 (MQ=255)
                                        cccGTGATTTATTCCCCCCAACACCCGCGCTTTTCCCCAACGCCAGGTGCAAAAAGCCCTCCGCGGCGGTTGTCGGCCGTCTTTTTGTCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                                             <  2:418636/138‑1 (MQ=255)
                                                   aTTCCCCCCAACACCCGCGCATTTCTCCACCGCCATGTGAAAAAATGCCTCCGCGGCGCCTGCCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg                                                                    <  1:650448/138‑1 (MQ=255)
                                                     tccccccAACACCCGCCAATTTATCCACCCCCTTGTGCAAAAATGCCTCCGCGGCGGCTTTCTGCAAGCTATTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc                                                                   <  1:351866/132‑1 (MQ=255)
                                                                         ttATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGtt                                                                                                           >  2:713034/1‑78 (MQ=255)
                                                                         ttATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGAAACAAAACCGCCGtt                                                                                                           <  1:713034/78‑1 (MQ=255)
                                                                          tATCCAACGCGATGTAAAAAATGGCCTTGGGGGCGGCGTTCTCCCAGTTATTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                             <  1:226522/138‑1 (MQ=255)
                                                                            cccAACGCCTGGTGAAAAAATGCCTCCGCGGCGGTTGCCGGCCAGTTATTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct                                           <  2:264683/137‑1 (MQ=255)
                                                                              cAACCCAATGTGAAAAATTGCCTCGGCGGGGGTTTTCTGCCACTTATTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                        <  1:238661/139‑1 (MQ=255)
                                                                                 cGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGGAAAAACATTTAACCCCCCCCTTTGTGTGGGAATAGCCGCTggg                                     >  2:593107/1‑138 (MQ=25)
                                                                                       gtgAAAAATTGCCTCCGGGGCGCCTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                               <  1:365598/139‑1 (MQ=255)
                                                                                       gaaaaaaaaCCCCTCCGGGGCGTCTTCCTGCAAGCTTTTTCCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                               <  2:7730/135‑1 (MQ=255)
                                                                                                                ttCTGCCGCTTTTTTTCCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc       <  1:763241/137‑1 (MQ=255)
                                                                                                                    gCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc  <  2:173686/139‑1 (MQ=255)
                                                                                                                                           |                                                                                                                     
CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTC  >  NZ_CP009273/353464‑353717

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: