Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,383,7220GA34.2% 56.9 / 17.4 38F80F (TTC→TTT) ycjWLacI family putative transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (24/1);  new base A (0/13);  total (24/14)
Fisher's exact test for biased strand distribution p-value = 2.59e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTT‑CGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATAAA  >  NC_000913/1383488‑1383948
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                  
ggCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTT‑CGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTT‑GGCGGTGTTGGCATGCAAAACTTAAGGGaaaa                                                                                                                                                                                                                                >  1:404648/1‑239 (MQ=255)
                                       aaTATCATATTTACCGTGTTGGAAAAAAAAGCACTCATCCGGGGGAGAACTGAACATAATAATCCCTTTAATCTTTTTCGGTTTAATTTTGTTTT‑CGATTTTTTGAAGTTCTTCTGCCGGATTGGGCGACGTTTGAAATAACACGTCAAAACCTTCTTCTTCAGCTTG‑GGGGGTGATGGCATGAAAAACTTAAGAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGaa                                                                                                                                                                                          <  2:107384/238‑1 (MQ=255)
                                                   aTCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTT‑CGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATg                                                                                                                                                                              >  1:781894/1‑238 (MQ=255)
                                                     cGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTT‑CGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGc                                                                                                                                                                            >  1:809781/1‑238 (MQ=255)
                                                                                                aataatGCCTTTAATCATTTTCTGCTTAATTTTGCTTT‑CGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTg                                                                                                                                   >  1:419275/1‑236 (MQ=255)
                                                                                                       ccTTTAATCATTTTCTGCTTAATTTTGCTTT‑CGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGAt                                                                                                                             >  1:283566/1‑235 (MQ=255)
                                                                                                        cTTTAATCATTTTCTGCTTAATTTTGCTTT‑CGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAAAGGATTCCCCGCCTTAGTTTTGGTCGAACGGGTAGAAATCACCAAAAACGCCCACCCACCCGAAGAGGGTAGCGCGAGGGCCAGCTTGGGTGGCTGAt                                                                                                                             >  2:142002/1‑234 (MQ=255)
                                                                                                          ttAATCTTTTTCGTTTTATTTTTGCTTT‑CGCTTTGTGGCAGTTCTTCTGCCGTATTGTGCGTCTTTTGTAATATCACGTCAAAACCTTCTTCTTCAGTTTG‑GGCGGTGATGGCATGAAAAACTTCAGAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTaa                                                                                                                      <  2:208197/239‑1 (MQ=255)
                                                                                                           tAATCATTTTCTGCTTAATTTTGCTTT‑CGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAAt                                                                                                                     >  1:208197/1‑239 (MQ=255)
                                                                                                               tttttttgtgttAATTTTTCTTTCCGTTTTTTTCCAATTCTTCGGCCGTTTTGTGCAACGTTTGAAATATCCCGTCAAACCCTTCTTCTCCAGCTGT‑GGCGGTAATGGCATGAAAACCTAAAAAAAAAAACGAATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCGTCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAAtt                                                                                                                    <  2:868954/227‑1 (MQ=255)
                                                                                                                 ttttCTGCTTAATTTTGCTTT‑CGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGtt                                                                                                                                     >  1:848491/1‑217 (MQ=255)
                                                                                                                  tttCTGCTTAATTTTGCTTT‑CGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCaa                                                                                                              >  1:913621/1‑239 (MQ=255)
                                                                                                                        cTTAATTTTGCTTT‑CGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGc                                                                                                          >  1:868954/1‑237 (MQ=255)
                                                                                                                                    tt‑CTCATTTTGGCATATCTTCTGCCGGTTGGTGCGATGTTGGAAATATCACGTCGAAACCTTCTTCTCCAGCTTG‑GGGGGTGATGGCAGAAAAACCTTCAGAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGt                                                                                                                                      <  2:764315/197‑1 (MQ=255)
                                                                                                                                               tGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATa                                                                                                                            >  1:764315/1‑197 (MQ=255)
                                                                                                                                                      cttctGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGgag                                                                           >  1:891392/1‑239 (MQ=255)
                                                                                                                                                                           tttGTAATATCCCTTGCAACCCTTCTTCTTCAGCTTG‑GGCGGTGAGGGAATGAAAAACTTAAGAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCa                                                                                                                                            <  2:773841/153‑1 (MQ=255)
                                                                                                                                                                                  tatCAGCTCAAACCCTTCTTCTTCAGTTTT‑GGCGGTGAGGCCATCAAAACCTTCAGAAAAAACCGTATCACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTcgcg                                                                                        <  2:607517/198‑1 (MQ=255)
                                                                                                                                                                                     cACGTCGAAACCTTCTTCTTCAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTgg                                                                                                                                  >  1:773841/1‑153 (MQ=255)
                                                                                                                                                                                       cGTCGAAACCTTCTTCTTCAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACg                                          >  1:762408/1‑239 (MQ=255)
                                                                                                                                                                                        gtaaaacccTTTTTCTTCAGTTTG‑GGCGGTGAGGGCATAAAAAACTTAAGAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGt                                                                                                                                      <  2:538193/139‑1 (MQ=255)
                                                                                                                                                                                         tCGAAACCTTCTTCTTCAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCCATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGt                                        >  2:756694/1‑239 (MQ=255)
                                                                                                                                                                                          agaagctttcttctt‑aGCTTTGGGCGTTGATGGCATGAAAAACTTCAAAAAAAAACGATTTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCa                                                        <  2:283566/215‑1 (MQ=255)
                                                                                                                                                                                            aaaCCTTCTTCTTCAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCggg                                                                              >  1:607517/1‑198 (MQ=255)
                                                                                                                                                                                            aaaCCTTCTTCTTCAGCTTT‑GGCGGTGATGGCATGAAAACCTTCAGAAAAAACCGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATc                                                                                                                                                                         <  2:465466/107‑1 (MQ=255)
                                                                                                                                                                                                  tcttcttcAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATa                                                                                                                            >  1:538193/1‑146 (MQ=255)
                                                                                                                                                                                                  tcttcttcAGCTTT‑GGCGGTGATGGCATGAAAATCTACAGAAAAAACCGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGa                               <  2:762408/239‑1 (MQ=255)
                                                                                                                                                                                                   cttcttcAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGa                                                                                                                                                                                                    <  1:561354/73‑1 (MQ=255)
                                                                                                                                                                                                   cttcttcAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGa                                                                                                                                                                                                    >  2:561354/1‑73 (MQ=255)
                                                                                                                                                                                                   cttcttcAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGAt                              >  1:821590/1‑239 (MQ=255)
                                                                                                                                                                                                     tcttcAGTTTG‑GGCGGTGATGGCATGAAAAACTTCAGAAAAAACCGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAgcgc                                                                                                                                                    <  2:531543/119‑1 (MQ=255)
                                                                                                                                                                                                      cttcAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGaa                                                                                                                                                               >  1:465466/1‑107 (MQ=255)
                                                                                                                                                                                                             tttgggCGGTGATGGCATGCAAACCTTCAGAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATAc                                                                                                                           <  2:444280/133‑1 (MQ=255)
                                                                                                                                                                                                               t‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGc                                                                                                                                          >  1:531543/1‑119 (MQ=255)
                                                                                                                                                                                                                       gtggacATGCAACTTTTCAGAAAACTGCGAATCACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCa            <  2:627015/233‑1 (MQ=255)
                                                                                                                                                                                                                       gATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAAttctt                                                                                                                 >  1:444280/1‑137 (MQ=255)
                                                                                                                                                                                                                       gATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCCTGGTTGGCTGATACTGTAATTCTTCAATGGCCCGCAACACCTTTTCCGTCGCTTCCGGGga                                                                            >  2:170112/1‑174 (MQ=255)
                                                                                                                                                                                                                                cAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATaaa  >  1:627015/1‑239 (MQ=255)
                                                                                                                                                                                                                                          gagaAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGa                                                                                                                              <  2:472165/105‑1 (MQ=255)
                                                                                                                                                                                                                                          gagaAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGa                                                                                                                              >  1:472165/1‑105 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                  
GGCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTT‑CGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTT‑GGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATAAA  >  NC_000913/1383488‑1383948

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑