Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730567,3890AT78.1% 62.1 / 15.1 32Q95L (CAA→CTA) ctsRCtsR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (4/3);  new base T (12/13);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GATTCACTAATGAACATGGTTATGAAATCGAAAGTAAACGTGGTGGTGGTGGTTACATCCGAATCACTAAAATTGAAAATAAAGATGCAACAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCAACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACGTGAAGCTAAAATGATTCAAGCAGTTATTGATAGAGAAACGCTATCAATGGATATGGTTTCTAGAGATATTATTAGAGCAAAT  >  CP000730/567251‑567528
                                                                                                                                          |                                                                                                                                           
gATTCACTAATGAACATGGTTATGAAATCGAAAGTAAACGTGGTGGTGGTGGTTACATCCGAATCACTAAAATTGAAAATAAAGATGCAACAGGGTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTctac                                                                                                                                           >  2:9988/1‑141 (MQ=255)
      cTAATGAACATGGTTATGAAATCGAAAGTAAACGTGGTGGTGGTGGTTACATCCGAATCACTAAAATTGAAAATAAAGATGCAACAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAg                                                                                                                                     <  1:135338/141‑1 (MQ=255)
       tAATGAACATGGTTATGAAATCGAAAGTAAACGTGGTGGTGGTGGTTACATCCGAATCACTGAAATTGAAAATAAAGATGCAACAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGc                                                                                                                                    <  1:308440/141‑1 (MQ=255)
                tGGTTATGAAATCGAAAGTAAACGTGGTGGTGGTGGTTACATCCGAATCACTAAAATTGAAAATAAAGATGCAACAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGCTtattatat                                                                                                                           >  2:268831/1‑141 (MQ=255)
                  gTTATGAAATCGAAAGTAAACGTGGTGGTGGTGGTTACATCCGAATCACTAAAATTGAAAATAAAGATGCAACAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGCTtattatatta                                                                                                                         >  2:393956/1‑141 (MQ=255)
                            cGAAAGTAAACGTGGTGGTGGTGGTTACATCCGAATCACTAAAATTGAAAATAAAGATGCAACAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGCTTATTATATTATTGATGGGCt                                                                                                               <  1:249219/141‑1 (MQ=255)
                              aaaGTAAACGTGGTGGTGGTGGTTACATCCGAATCACTAAAATTGAAAATAAAGATGCAACAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGCTTATTATATTATTGATGGGCttt                                                                                                             >  2:363206/1‑141 (MQ=255)
                                 gTAAACGTGGTGGTGGTGGTTACATCCGAATCACTAAAATTGAAAATAAAGATGCAACAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCAACAACAAGCTTATTATATTATTGATGGGCTTTTAg                                                                                                          >  1:377836/1‑141 (MQ=255)
                                   aaaCGTGGTGGTGGTGGTTACATCCGAATCACTAAAATTGAAAATAAAGATGCAACAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGCTTATTATATTATTGATGGGCTTTTAGAt                                                                                                        >  2:344255/1‑141 (MQ=255)
                                       gtggtggtggtgGTTACGTCCGAATCACTAAAATTGGAAATAAAGATGAAACAGGTTATAATAATCATTTGCTTCGGCTGATTGGACGTTCTATTTCTCTACAACAAGCTTATTATATTATTGATGGGCTTTTAGATaaaa                                                                                                    >  2:340359/1‑141 (MQ=255)
                                         ggtggtggtggtTACATCCGAATCACTAAAATTGAAAATAAAGATGCAACAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATg                                                                                                  <  2:377882/141‑1 (MQ=255)
                                           tggtggtggtTACATCCGAATCACTAAAATTGAAAATAAAGATGCAACAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGtt                                                                                                <  2:151640/141‑1 (MQ=255)
                                                       cATCCGAATCACTAAAATTGAAAATAAAGATGCAACAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACg                                                                                    >  2:181918/1‑141 (MQ=255)
                                                             aaTCACTAAAATTGAAAATAAAGATGCAACAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACGTGAAGc                                                                              <  1:9988/141‑1 (MQ=255)
                                                                 aCTAAAATTGAAAATAAAGATGCAACAGGTTATTTTAATCATTTGCTTCAGCTGATCGGACCTTCTATTTCTCTACAACAAGCTAATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACGTGAAGCTaaa                                                                          >  2:34834/1‑141 (MQ=255)
                                                                  cTAAAATTGAAAATAAAGATGCAACAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACGTGAAGCTaaaa                                                                         <  2:182983/141‑1 (MQ=255)
                                                                            acaTAAAGATGCAACAGGTTATATTAATCATTTGCTTAAGCTGATTGGACCTTCTATTTCTCAACAACGAGCTTATTATATGGCTGGTAGGCTGTTAGGTGAAATGTTGCTAGATGATAGTGAAGCTAAAATGATTCAGGc                                                               <  2:224689/139‑1 (MQ=255)
                                                                            aaaTAAAGATGCAACAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACGTGAAGCTAAAATGATTCAAGc                                                               >  1:412843/1‑141 (MQ=255)
                                                                            aaaTAAAGATGCAACAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCAACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACGTGAAGCTAAAATGATTCAAGc                                                               <  2:224721/141‑1 (MQ=255)
                                                                            aaaTAAAGATGCAACAGGTCATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCAACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACGTGAAGCTAAAATGATTCAAGc                                                               >  2:105421/1‑141 (MQ=255)
                                                                              aTAAAGATGCAACAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACGTGAAGCTAAAATGATTCAAGCAg                                                             <  2:229698/141‑1 (MQ=255)
                                                                                 aaGATGCAACAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACGTGAAGCTAAAATGATTCAAGCAGTTa                                                          >  2:120189/1‑141 (MQ=255)
                                                                                          cAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTtaa                                                                                              >  1:343404/1‑96 (MQ=255)
                                                                                          cAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTtaa                                                                                              <  2:343404/96‑1 (MQ=255)
                                                                                          cAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACGTGAAGCTAAAATGATTCAAGCAGTTATTGATagag                                                 <  2:122009/141‑1 (MQ=255)
                                                                                           aGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCAACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACGTGAAGCTAAAATGATTCAAGCAGTTATTGATAgaga                                                >  1:210150/1‑141 (MQ=255)
                                                                                              ttATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACGTGAAGCTAAAATGATTCAAGCAGTTATTGATAGAGAAAc                                             <  1:408200/141‑1 (MQ=255)
                                                                                               tataTTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGATTATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACGTGAAGCTAAAATGATTCAAGCAGTTATTGATAGTGAAACg                                            <  2:236/141‑1 (MQ=255)
                                                                                                 tatTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCAACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACGTAAAGCTAAAATGATTCAAGCAGTTATTGATAGAGAAACGCt                                          <  2:124865/141‑1 (MQ=255)
                                                                                                      aTCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACGTGAAGCTAAAATGATTCAAGCAGTTATTGATAGAGAAACGCCATCaa                                     >  2:127942/1‑141 (MQ=255)
                                                                                                            tGCTTCAGCTGATTGGACCTTCTATTTCTCTACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACGTGAAGCTAAAATGATTCAAGCAGTTATTGATAGAGAAACGCTATCAATGGata                               >  2:331199/1‑141 (MQ=255)
                                                                                                                   gCTGATTGGACCTTCTATTTCTCAACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACGTGAAGCTAAAATGATTCAAGCAGTTATTGATAGAGAAACGCTATCAATGGATATGGTTTc                        >  1:434152/1‑141 (MQ=255)
                                                                                                                          ggACCTTCTATTTCTCTACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACGTGAAGCTAAAATGATTCAAGCAGTTATTGATAGAGAAACGCTATCAATGGATATGGTTTCTAGAGat                 <  1:344255/141‑1 (MQ=255)
                                                                                                                                         caacaacaaGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACGTGAAGCTAAAATGATTCAAGCAGTTATTGATAGAGAAACGCTATCAATGGATATGGTTTCTAGAGATATTATTAGAGCAAat  >  2:284769/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                           
GATTCACTAATGAACATGGTTATGAAATCGAAAGTAAACGTGGTGGTGGTGGTTACATCCGAATCACTAAAATTGAAAATAAAGATGCAACAGGTTATATTAATCATTTGCTTCAGCTGATTGGACCTTCTATTTCTCAACAACAAGCTTATTATATTATTGATGGGCTTTTAGATAAAATGTTAATAAATGAACGTGAAGCTAAAATGATTCAAGCAGTTATTGATAGAGAAACGCTATCAATGGATATGGTTTCTAGAGATATTATTAGAGCAAAT  >  CP000730/567251‑567528

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑