Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,684,7120AC42.9% 88.2 / 15.1 56T152P (ACC→CCC) tusinhibitor of replication at Ter, DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (25/7);  new base C (0/24);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 8.52e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAACCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCAC  >  NC_000913/1684567‑1684862
                                                                                                                                                  |                                                                                                                                                       
tAACCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCtt                                                                                                                                                    >  1:165346/1‑151 (MQ=255)
          cAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTAcc                                                                                                                                          >  1:991376/1‑151 (MQ=255)
           aaGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCg                                                                                                                                         <  2:707374/151‑1 (MQ=255)
             gcagcGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCa                                                                                                                                       <  2:811069/151‑1 (MQ=255)
              cagcGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCAc                                                                                                                                      >  1:690691/1‑151 (MQ=255)
                gcGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGc                                                                                                                                    >  1:10035/1‑151 (MQ=255)
                     ggTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCAcc                                                                                                                               >  1:736304/1‑151 (MQ=255)
                       tcagtcaTATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGt                                                                                                                             >  2:183212/1‑151 (MQ=255)
                          gtcaTATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCt                                                                                                                         >  2:1006177/1‑151 (MQ=255)
                             atatTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGc                                                                                                                       >  1:1023323/1‑150 (MQ=255)
                                  cAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGAcc                                                                                                                 <  1:163913/151‑1 (MQ=255)
                                     cacaTCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCg                                                                                                              <  1:133224/151‑1 (MQ=255)
                                      acaTCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGc                                                                                                             <  1:436481/151‑1 (MQ=255)
                                        aTCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCa                                                                                                           <  2:437454/151‑1 (MQ=255)
                                              aaaCTCAAGCCCACGTTCGAGCATATCGTCAGGTTGGAATCAGAACCCCCCCCCGCGGCATGTTTTGAAGGGGTGCATCGTCTTTTGCCGGGGCTGCCCCCCTTTAATGCTCACCGCACGCTCCCCGTTCTGCCCGCCCCCGCCACTTTa                                                                                                      <  1:585622/150‑1 (MQ=255)
                                                    aaGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTttg                                                                                               >  2:333907/1‑151 (MQ=255)
                                                        ccccGTTGGAGAATTTCGTCAGGTTTAAATCAGACCCCCCCCCCGCGCCACGTTTTAAATGGGTGCATCTTCTTTTGCCGGGGCTAACCCCCTTTATTTCTTCCCGCCCGCTCCCCGTTCTGCACGACCCCGCCACTTTACGCTttgg                                                                                              <  1:762229/145‑1 (MQ=255)
                                                        ccccGTTCGAGCATATCGTCACGGTTAAATCAAACCCCCCCCCCGCGGCACTTTTTGAAGGGGTGCTTCGTCTTTTGCCGGGGCTGATCCCCTTTAATGCTCCCCGCACGCCCCCCTTTCTGCACGCCCCCGCCACTTTACGCTttggttg                                                                                           <  2:140287/148‑1 (MQ=255)
                                                          aCGTTCGAGCATAGCGTCAGGGTTGAATAAGACCCCCCCCCCGCGGCGCTTTTTGAAGGGGTGCATCGTCATTTGCGGGGGCTGACCACCCTTAATGCTTCCCGCCCGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTggg                                                                                         <  2:649676/151‑1 (MQ=255)
                                                           cGTTCGAGCATATGGCCCCGGTTAAATCAAACCCCCCCACCGCGGCGCTTTTTGAATGGGTGCGTCGTCATTTGCCGGGGCTGATCCCCCTAAATGCTTCCCGCCCGCTCCCCTTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGc                                                                                        <  2:794855/151‑1 (MQ=255)
                                                             ttCGAGCATATCGTCAGGGTTGAATCAGACCCCCCCCCCGCGCCCCTTTTTGAAGGGGTGCATCGTCTTTTGCGGGGGCTGATCCCCTTTAATGCTTCCCCCACGCCCCCCTTTCTCCACGCCCCCGCCACTTTACGCTTTGGTTggg                                                                                         <  2:589085/148‑1 (MQ=255)
                                                               cgaccatacgtgccgggTTGACTCAGAACCCCCCCCCGCGGCATTTTTAAATGGGTTGGACCTTCATTTCCCGGGGCTCATCCCCCTTAATGCTTCCCGCCCCCTCCCGTTACTCCACGCCCCCGCCACTTTACGCTTTGGTTGGGCtaat                                                                                    <  1:988476/136‑1 (MQ=255)
                                                               cGAACATATCGTCACGGTTAAATCAGAACCCCCCCCCGCGGCCCGTTTTGAAGGGGTGCATCGTCATTTGCGGGGGCTGATCCCCCTTATTGCTTACCGCCCGCTCACCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCtaat                                                                                    <  2:388634/151‑1 (MQ=255)
                                                                gAGCATATCGTCACGGTTAAATCAGACCCCCCCCCCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCCCCCTTAATGCTTACCGCACGCTCCCCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGc                                                                                        <  2:359596/146‑1 (MQ=255)
                                                                   cATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAAc                                                                                >  2:578491/1‑151 (MQ=255)
                                                                             cGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTa                                                                       >  2:205510/1‑150 (MQ=255)
                                                                                ttaattCAGACCCCCCCCCCGCGGCAGGTTTTGAGGGGGTGCATCGTCTTTTCCCGGGCCTGATCCCCATTACTGCTAACCGCTCGCTCCCCTTTCTGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCAt                                                                         <  2:683452/140‑1 (MQ=255)
                                                                                   aaacaaaaCCCCCCCCCGCGCCACTTTTGGAAGGGGTGCATCGTCTTTGGCCGGGGCTGACCCCCTTTATTGCTTCCCGCACGCCCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAtt                                                                 <  1:87957/144‑1 (MQ=255)
                                                                                   aaTCAGACCCCCCCCCCGCGCCCCTTTTTGAATGGGTCCATCGTCTTTTGCCGGGGCTGACCCCCCTTAATGCTTACCGCGCGCTCCCCGTTCTGCACCCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGa                                                                    <  1:21423/147‑1 (MQ=255)
                                                                                   aaTCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAttt                                                                >  1:651352/1‑151 (MQ=255)
                                                                                   aaTCAGAACCCCCCCCCGCGGCACTTTTTGAAGGGGTGCATCGTCATTTGCGGGGGCTGATCCCCCTTAATGCTTCCCGCACGCTCCCCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAttt                                                                <  2:679744/151‑1 (MQ=255)
                                                                                      caaacccccccccCGCGGCATGTTTTGAAGGGGTGCTTCGTCATTTGCCGGGGCTGATCCCCTTTATTGCTTCCCGCACGCTCACCTTTCTCCCCGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAtt                                                                 <  1:233025/140‑1 (MQ=255)
                                                                                       aGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACAt                                                            >  1:659537/1‑151 (MQ=255)
                                                                                       aGAACCCCCCCCCGCGGCATTTTTTGAATGGGTGCATCGTTTTTTGCGGGGGCTGATCCCCTTTAATTCTTCCCGCACGCCCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTAc                                                              <  1:977713/149‑1 (MQ=255)
                                                                                         acccccccccCGCGGCGCTTTTTGAGTGGGTGCATCGTCATTTGCGGGGGCTGTTCCCCCTTAATGCTTCCCGCCCGCTCCCTTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCg                                                          <  1:1095389/147‑1 (MQ=255)
                                                                                          aCTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGt                                                         >  1:608517/1‑151 (MQ=255)
                                                                                           ccccccACCGCGGCACGTTTTGAGTGGGTCCACCCCCATTTGCGGGGCCTGACCCCCTTTATTGCTTCCCGCCCGCCCCCCGTTCTCCCCCACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCg                                                          <  1:738590/147‑1 (MQ=255)
                                                                                            tCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAGTAAACATATCATTAAGAATTTACATCGt                                                         >  1:418462/1‑149 (MQ=255)
                                                                                               ccccggcgGCCCTTTTTGAGTGGGTGCGTCGTCTTTGGCGGGGGCTGACCCCCCTTATTTTTTCCCGCGCGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtgatga                                                    <  1:391557/146‑1 (MQ=255)
                                                                                                       gCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTgg                                            >  1:649556/1‑151 (MQ=255)
                                                                                                        cACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGc                                           >  2:581199/1‑151 (MQ=255)
                                                                                                           gTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGcac                                         >  1:540331/1‑150 (MQ=255)
                                                                                                            ttatGAATGGGTGCACCGCAATTCGCCGGGGCGAATCCCCCTTATTCCTTCCCCCACGCTCTCCTTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGcac                                         <  2:402104/146‑1 (MQ=255)
                                                                                                                 acggggTGCTTCCTCATTTCCCGGGCCTGATCCCCCTTAATGTTTCCCCCCCGCTCCCCTTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGc                                      <  2:503328/144‑1 (MQ=255)
                                                                                                                     ggTGCATCGTCATTGGCCGGGGCTGATCCCCTTTAATGCTTACCGCACTCTCCCCGTTCTGCACCCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGaaaaa                              <  1:678593/151‑1 (MQ=255)
                                                                                                                      gtttcttCGTCTTTTGCGGGGGCTGATCCCCTTTAATGCTTCCCGCCCGCCCCCCGTTCGGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTgg                                   <  2:275210/140‑1 (MQ=255)
                                                                                                                       gcgaatcttaaTTCGCCGGGCCTGTCCCCCTTTATTGCTTCCCCCACGCCCACCTTTCTCCACAACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGaaaaaa                             <  1:68337/141‑1 (MQ=255)
                                                                                                                             cGTCTTTTGCCGGGGCTAATCCCCCTTAATGCTTCCCGCCCGCTCCCCTTTCTCCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAg                            <  1:686385/146‑1 (MQ=255)
                                                                                                                                  tttGCCGGGGCTGATCCCCCTTAATGCTTCCCGCACGCCCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAAt                    <  2:538618/149‑1 (MQ=255)
                                                                                                                                   ttGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATcac                 >  1:713266/1‑151 (MQ=255)
                                                                                                                                         gggCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGc            >  2:282101/1‑150 (MQ=255)
                                                                                                                                            cTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAgtg        >  2:138344/1‑151 (MQ=255)
                                                                                                                                             tGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAgtgt       >  1:211341/1‑151 (MQ=255)
                                                                                                                                              gATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTc      >  1:465602/1‑151 (MQ=255)
                                                                                                                                               aTCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCg     >  1:1127287/1‑151 (MQ=255)
                                                                                                                                                 cACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGc    >  2:852835/1‑150 (MQ=255)
                                                                                                                                                  aCCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCAc  >  1:165374/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                       
TAACCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCAC  >  NC_000913/1684567‑1684862

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