Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273353,5910TC52.2% 9.8 / 32.7 23K271E (AAA→GAA) cynRtranscriptional regulator CynR
Reads supporting (aligned to +/- strand):  ref base T (9/2);  new base C (0/12);  total (9/14)
Fisher's exact test for biased strand distribution p-value = 6.73e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.99e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATC  >  NZ_CP009273/353454‑353712
                                                                                                                                          |                                                                                                                           
aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc                                                                                                                            >  2:372277/1‑139 (MQ=255)
   aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCGGTGTCAGCGACAACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAGGTGCAAAAATGCCGTCGCGGCGGCTGTCTGCCAGCTATTTTTccgc                                                                                                                         <  2:149940/139‑1 (MQ=17)
                                       tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACCAAACCGCCGTTCTCTCCAGCAGTGGCggg                                                                                     >  2:734010/1‑139 (MQ=255)
                                       tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                     >  2:61358/1‑139 (MQ=255)
                                            gCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggca                                                                                 >  1:544244/1‑135 (MQ=255)
                                                   ccGTGTTTCATTCCCGCAAACACCCGCGCTTTTACCCAACGCCAGGTGAAAAAATGCCTCCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                                        <  1:401944/138‑1 (MQ=25)
                                                            ctttCCGCCAACCACCCCGTCTTTTACCCACCGCCTGTTGCAAAAATGCTTTGGGGGCGGTGGTCGGCCAGCTTTTTACCCGCCGCAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa                                                                <  2:398924/137‑1 (MQ=255)
                                                              tcccGCCAACACCCGCCCTTTTACCCACCGCCATGTGAAAAAATCCCTGCGGGGGGGCTGTCTGCCAGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc                                                              <  1:372277/135‑1 (MQ=255)
                                                                ccccaaacaccCGCGTATTTATCCAACCCAATGTGAAAAAATGCTTTCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCg                                                           <  1:436079/133‑1 (MQ=255)
                                                                 ccccAACACCCGCCCATTTACCCACCGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTc                                                         <  1:333120/137‑1 (MQ=38)
                                                                   ccAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAAATTTTAAACCCgcct                                                        >  1:637294/1‑136 (MQ=25)
                                                                            gcgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAACATTTAACCCCCCCTGTTGTGTgg                                              >  1:286958/1‑139 (MQ=25)
                                                                            gcgcATTTACCCACCGCCATGTCAAAAAATGCCTTCGGGGCGGTTGTCTGCCAGCTATTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTgg                                              <  2:651951/139‑1 (MQ=255)
                                                                                    tacccACCGCTATGTGAAAAAATGCCTCCGCGGCGGCTGTCTGCAATTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                        <  2:637294/135‑1 (MQ=255)
                                                                                       ccAACCCCTTGTAAAAAAATCCCTCCGGGGGGGCTGTCTGCCAGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgc                                    <  2:223113/139‑1 (MQ=255)
                                                                                           cGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGGGAAAAAACCTTTAACCCCCCCTGTTTTGTGGGAAACGCTGCTgg                                 >  1:545796/1‑138 (MQ=25)
                                                                                            cccATGTGCAAAAATCCTTTCGGGGCGCCTGCTTGCCAGCTATTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                               <  1:315916/138‑1 (MQ=255)
                                                                                             ccttgtgCAAAAATCCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa                              <  1:703754/136‑1 (MQ=255)
                                                                                                tgtgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCGGAGAAAAAAACTTAATGCCCCTCTGTTGTGTGGGAAGCGCAGCTGGTTACa                            >  2:722682/1‑138 (MQ=25)
                                                                                                tgtgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAACTTTTAACCCGCCATGTTTTGTTGGGATAGCTGCCggggacc                            >  1:505468/1‑133 (MQ=25)
                                                                                                  tgCAAAAACGCTTCCGCGGCTGCTGCCCGCCCTCTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                          <  1:154155/138‑1 (MQ=255)
                                                                                                        atgcccTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa                   <  2:505468/135‑1 (MQ=255)
                                                                                                                           tCTGCCAGCTATTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc  <  1:584260/137‑1 (MQ=255)
                                                                                                                                          |                                                                                                                           
ATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATC  >  NZ_CP009273/353454‑353712

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