Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,590,9850GT40.5% 49.4 / 11.3 37intergenic (+94/+30)PP_1394/PP_1395acetolactate synthase large subunit/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (14/8);  new base T (0/15);  total (14/23)
Fisher's exact test for biased strand distribution p-value = 7.00e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGGACCAGAATCAGTGGCAGGCCGCGGTGTTGGGTTGAAGCATTGCGATTGCCTGTTTCGGCCTCTCGCGCCTGCAGGCAGCGCGCAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAAC  >  NC_002947/1590854‑1591112
                                                                                                                                   |                                                                                                                               
gggACCAGAATCAGTGGCAGGCCGCGGTGTTGGGTTGAAGCATTGCCATTGCCTGTTTCGGCCTCTCGCGCCTGCAGGCAGCGCGCAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTTACCCGCGAAAg                                                                                                                             <  2:693326/135‑1 (MQ=255)
  gACCAGAATCAGTGGCAGGCCGCGGTGTTGGGTTGAAGCATTGCGATTGCCTGTTTCGGCCTCTCGCGCCTGCAGGCAGCGCGCAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTTACCCGCTAAAgcg                                                                                                                           <  2:1361616/136‑1 (MQ=255)
     aaGAATCTGTGGCCGGCGGCGGTGTTGGGTTGAATCTTGGCGATTGCCTGTTTCGGCCTCTCGCGCCTGCAGGCAGCGCGAATTTTTCCCATCTGTGAAAAAACCTGTTGGTGCGGGTTTACCCGCTAAAGCGCCg                                                                                                                        <  2:56987/135‑1 (MQ=255)
          tCAGTGGCAGGCCGGGGTGTTGGGTTGAAGCATTGCGATTGCCTGTTTCGGCCTCTCGCGCCTGCAGGCAGCGCGCAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTTACCCGCTAAAGCGCCGGTGCa                                                                                                                   <  1:382117/136‑1 (MQ=255)
                 cAGGCCGCGGTGTTGGGTTTAAGCATTGCGATTGCCTGTTTCGGCCTCTCGCGCCTGCAGGCAGCGCGCAATTTTCGCATCTGTCCAAAAACCTGTAGGAGCGGGTTTACCCGCTAAAGCGCCGGTGCAGGAGAAt                                                                                                            <  1:98841/136‑1 (MQ=255)
                     ccGCGGTGTTGGGTTGAAGCATTGCGATTGCCTGTTTCGGCCTCTCGCGCCTGCAGGCAGCGCGCAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCAGGAGAATGTTg                                                                                                        >  1:124378/1‑136 (MQ=255)
                       gcgGTGTTGGGTTGAAGCATTGCGATTGCCTGTTTCGGCCTCTCGCGCCTGCAGGCAGCGCGCAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCAGGAGAATGTTGAg                                                                                                      >  1:299825/1‑136 (MQ=255)
                                 gTTGAAGCATTGCGATTGCCTGTTGCGGCCGCTCCCGCCTGCAGGCAGCGCGCAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTAACCCGCTAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCg                                                                                            <  1:172499/136‑1 (MQ=255)
                                 gTTCAACCAGTGCGATTCCCGGTTTCGGCCTCTCGCGCCTGCTGGCAGCGCGAAATTTCCCCAGCTGTCAAAAAACCTGTAGCCGGGGGTTTACCCGCTAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCg                                                                                            <  2:822490/136‑1 (MQ=255)
                                   tGAAGCATTGCGATTGCCTGTTGCGGCCTCTCGCGCCTGCTGGCAGCGCGAAATTTTCGCATCTGTGCAAAAACCTGTCGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCa                                                                                          <  2:1054614/136‑1 (MQ=255)
                                      tgCATTGCGATTGCCTGTTGCGGCCTCTCGCGCGTGTAGGCAGCGCGCAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTTACCCGCTAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTa                                                                                       <  1:400832/135‑1 (MQ=255)
                                        cATTGCGATTGCCTGTTTCGGCCTCTCGCGCCTGCAGGCAGCGCGCAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGc                                                                                     >  1:1172962/1‑136 (MQ=255)
                                        cATTGCGATTGCCTGTTTCGGCCTCTCGCGCCTGCAGGCAGCGCGCAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGc                                                                                     >  1:1024154/1‑136 (MQ=255)
                                          ttGCGATTGCCTGTTTCGGCCTCTCGCGCCTGCAGGCAGCGCGCAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGccc                                                                                   >  2:1127983/1‑136 (MQ=255)
                                          ttGCGATTGCCTGTTTCGGCCTCTCGCGCCTGCAGGCAGCGCGCAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGccc                                                                                   >  2:1299636/1‑136 (MQ=255)
                                             cGATTGCCTGTTTCGGCCTCTCGCGCCTGCAGGCAGCGCGCAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCa                                                                                >  1:492439/1‑136 (MQ=255)
                                                     tGTTTCGGCCTCTCGCGCCTGCAGGCAGCGCGCAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCAACCGAGcc                                                                        >  1:93763/1‑136 (MQ=255)
                                                                     cccTGCAGGCAGCGCGCAATTTGCGCATCTGTGCAAAAACCTGTCGGAGCGGGTTAACCCGCTAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCAACCGAGCCACATCGCGAATcggcg                                                        <  1:1248798/135‑1 (MQ=255)
                                                                           aggcagCGCGAAATTTTCGCATCTGAGAAAAACCGTGTCGGCGGGGTTTTACCCGCGAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCAACCGAGCCACATCGCGAATCGGCGgcgcgc                                                  <  2:975692/136‑1 (MQ=255)
                                                                              cagCGCGCAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGt                                               <  2:1371555/136‑1 (MQ=255)
                                                                                   cgcAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGc                                          <  2:492439/136‑1 (MQ=255)
                                                                                   cgaatttttCGCTTCTGTGCAAAAACCTTTCGGTGGGGGTTTACCCGCTAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGc                                          <  2:314693/131‑1 (MQ=255)
                                                                                   cgaaaTTTTCGCATCTGTCAAAAAACCTGTAGCAGCGGGTTTACCCGCTAACGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGc                                          <  2:289442/133‑1 (MQ=255)
                                                                                    gcAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCa                                                                                                                   <  2:1187802/62‑1 (MQ=255)
                                                                                    gcAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCa                                                                                                                   >  1:1187802/1‑62 (MQ=255)
                                                                                    gcAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGcc                                         >  1:783274/1‑136 (MQ=255)
                                                                                     cAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCg                                        <  2:1024154/136‑1 (MQ=255)
                                                                                     aaaTTTTCGAATCTGTCAAAAAACCTGTAGGAGCGGGTTTACCCGCTAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCTCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCg                                        <  2:1172962/135‑1 (MQ=255)
                                                                                                    ttcAAAAACGTGTAGGAGCGTTTTAACCCGCTAACGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCt                         <  1:1030163/134‑1 (MQ=255)
                                                                                                          aaCCTGTAGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTg                   <  2:783274/136‑1 (MQ=255)
                                                                                                           aCCTGTAGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGc                  >  1:583902/1‑136 (MQ=255)
                                                                                                               gttggAGCGGGTTTACCCGCTAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACg              <  1:437899/133‑1 (MQ=255)
                                                                                                               gTAGGAGCGGGTTTACCCGCTAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACg              <  2:1372940/136‑1 (MQ=255)
                                                                                                                 aGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACggg            >  1:943847/1‑136 (MQ=255)
                                                                                                                 aGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACggg            >  1:757723/1‑136 (MQ=255)
                                                                                                                          tttACCCGCGAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATaa   >  2:282483/1‑136 (MQ=255)
                                                                                                                           ttACCCGCTAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAac  <  2:782986/136‑1 (MQ=255)
                                                                                                                                   |                                                                                                                               
TGGACCAGAATCAGTGGCAGGCCGCGGTGTTGGGTTGAAGCATTGCGATTGCCTGTTTCGGCCTCTCGCGCCTGCAGGCAGCGCGCAATTTTCGCATCTGTGCAAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGCGCCGGTGCAGGAGAATGTTGAGCTTCAGACCGCAGTAGCCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAAC  >  NC_002947/1590854‑1591112

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑