Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,684,7490AC48.9% 48.6 / 13.2 48D164A (GAC→GCC) tusinhibitor of replication at Ter, DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (15/9);  new base C (0/23);  total (15/33)
Fisher's exact test for biased strand distribution p-value = 2.39e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGA  >  NC_000913/1684601‑1684896
                                                                                                                                                     |                                                                                                                                                   
cAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGAcc                                                                                                                                                   <  2:712101/151‑1 (MQ=255)
          aTAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACttt                                                                                                                                         >  2:135716/1‑151 (MQ=255)
                            tCGGGCATACCGTCACGGTTAAACCAAACCCCCCCCCCGCGGCACGTTTTGAATGGGGGCATCGTCATTTCCCGGGGCTGATCCCCCTTAATTCTTACCGCACGCCCCCCGTTCTCCACGACCCCGCCACTTTACGCTTTGGTTggg                                                                                                                           <  2:594827/147‑1 (MQ=255)
                             cGAGCATATCGTCAGGGTTGAATCAGACCCCCCCCCCGCGGCACTTTTTGATTGGGTGCATCGTCTTTTGCGGGGGCTGATCCCCCTTAATGTTTCCCGCCCGCTCCCCTTTCTGCCCGACCCCGCCACTTTACGCTTTGGTTGGGCtaat                                                                                                                      <  2:889822/151‑1 (MQ=255)
                                    atCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACat                                                                                                                >  2:45890/1‑150 (MQ=255)
                                      cgtcggggTTGAAACAAACCCCCCCCGGGCGGCGCGTTTTGAAGGGGTGCATCGTCTTTGGCGGGGTCTAACCCCCCTTATTGCTTCCGGCGGGCCCACCGTCCTCAGGGTCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACat                                                                                                                <  1:702337/142‑1 (MQ=255)
                                      cGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACatat                                                                                                              >  1:626503/1‑150 (MQ=255)
                                      cGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATc                                                                                                             >  2:12/1‑151 (MQ=255)
                                       gtcggggTGGAACCAGACCCCCCCCCCGCGGCACGTTTTGAAGGGGTGCATCGTCATTTGCGGGGGCTGACCCCCCTAAATGCTTCCCGCCCGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACatat                                                                                                              <  2:918752/144‑1 (MQ=255)
                                       gcccgggTTAAATCAAACCCCCCCCCCGCGGCCCTTTTGAAATGGGTGCATCGTCTTTTGCGGGGGCTGACCCCCCTTATTGCTTCCCGCACGCTCACCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCa                                                                                                            <  1:305078/146‑1 (MQ=255)
                                       gTCAGGGTTGAACCAGACCCCCCCCCCGCGCCACTTTTTAATTGGGTGCATCGTCTTTTGCCGGGGCTGACCCCCTTTAATGCTTCCCGCCCCCCCACCGTTCTGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATc                                                                                                             <  1:871081/150‑1 (MQ=255)
                                         cagggTTAAACCAGACCCCCCCCCCGCGCCCCTTTTTGAATGGGTGCATCGTCTTTGGCCGGGGCTGACCCCCCTTAATGCTTCCCGCCCGCTCCCCTTTCTCCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACa                                                                                                                 <  1:715798/141‑1 (MQ=255)
                                         cACGGTTGAACCAGACCCCCCCCCCGCGCCCCTTTTTGAATGGGTGCATCGTCTTTTGCCGGGGCTGATCCCCCTTAATGCTTCCCGCCCGCTCCCCTTTCTCCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCAt                                                                                                           <  2:859670/150‑1 (MQ=255)
                                          aCGGTTAAATCAGACCCCCCCCCCGCGCCTCTTTTTGAATGGGTGCCTCTTCTTTTGCGGGGGCTGACCCCCCTTATTTTTTCCCGCACGCTCACCGTTCCCCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCa                                                                                                            <  1:231574/148‑1 (MQ=255)
                                            ggTTGAACAGGACCTCCCCCCCGGGGCCGTTTTTGAGGGGGTGCATCTTCTTTGGCCGGGGCGAATCCCCCTAAAGGCTACCCCCACGCTCCCCTTCCGCCGCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAg                                                                                                       <  2:34530/151‑1 (MQ=255)
                                               tGAAACAAAACCCCCCCCCGGCGGCATTTTTTAAGTGGGTGCCTCTTCTTTTGCCGGGGCGGCCCCCCTTTAATTCTTCCCGCCCGCCCCCCTTTCTCCACACCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCAtt                                                                                                          <  1:346908/146‑1 (MQ=255)
                                                 aataaaacccccccccCGCGGCATGTTTTGAATGGGTGCATCGTCTTTTGCCGGGGCTGATCCCCCTTAATGCTTCCCGCACGCTCCCCTTTCTCCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAtt                                                                                                   <  2:896228/140‑1 (MQ=255)
                                                 aaccaaacccccccccCGCGGCACTTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGACCCCCCTTATTGCTTACCGCCCGCCCACCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAttt                                                                                                  <  1:166834/141‑1 (MQ=255)
                                                 aaTCAGACCCCCCCGCCGCGGCACTTTTTGAAGGGGTGCATCGTCTTTTGCCGGGGCTGACCACCCTTATTGCTTCCCGCCCCCTCACCTTCCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAtt                                                                                                   <  1:375910/150‑1 (MQ=255)
                                                   tCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTAc                                                                                                >  1:140383/1‑151 (MQ=255)
                                                            cccACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtgatg                                                                                       >  1:123431/1‑151 (MQ=255)
                                                             ccccccgggccggtttttGATTGGGTGTACCGTCTTTTCCGGGGGCTGACCCCCCTAATGGTTTCCCGCACGCTCACCGTCCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtgatg                                                                                       <  2:12509/137‑1 (MQ=255)
                                                               cccGCGCCCCGTTTTGAATGGGTCCTTCGTCATTTGCGGGGGCTGTTCCCCCTTAATGCTCCCCGCCCACCCACTGTTCGCCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAg                                                                                    <  2:197074/150‑1 (MQ=255)
                                                                 cgcgGCCCGTTTTGAATGGGTGCATCGTCTTTGGCCGGGGCTGACCCCCCTTATTGCTTCCCGCGCGCTCCCCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTc                                                                                  <  1:233679/151‑1 (MQ=255)
                                                                     ggCACTTTTTGAATGGGTGCATCGTCTTTTGCCGGGGCTGACCCCCCTTAATGCTTCCCGCCCGCTCCCCGTTCTGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTg                                                                               <  2:21961/151‑1 (MQ=255)
                                                                                aagggggGCTTCGTCTTTGGCGGGGGCTGTTCCCCCTTATTTCTTCCCGCCCGCTCCCCTTTCTCCCCACCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTgg                                                                     <  1:935206/147‑1 (MQ=255)
                                                                                aaggggTGCCCCTTCTTTTCCCGGGGCGATTCACCCTTATTGTTTGCCGGACGTCCCCTGTTCTCCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCt                                                                       <  1:941331/145‑1 (MQ=255)
                                                                                   gggTGCATCTTCTTTTGCCGGGGCTGATCCCCCTTATTCCTTCCCGCACGCTCCCCTTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGaaa                                                                  <  1:370993/150‑1 (MQ=255)
                                                                                   gggTGCATCGTCTTTTGCCGGGGCTGACCCCCTTTAATGCTTCCCGCCCGCTCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGaaa                                                                  <  1:303582/150‑1 (MQ=255)
                                                                                                tttGCCGGGGCTGATCCCCCTTAATGCTTCCCGCCCCCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGaaa                                                       <  2:464447/148‑1 (MQ=255)
                                                                                                      ggggCTGATCCCCTTTAATGCTTCCCGCCCCCTCACCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATcacc                                                  <  2:772252/147‑1 (MQ=255)
                                                                                                      ggggCTGATCCCCCTTAATGCTTACCGCCCGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATcacca                                                 <  2:253773/148‑1 (MQ=255)
                                                                                                       ggccTGATCACCTTAAATCCCTCCCCCACGCTCCCCTTTCTCCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCa                                             <  2:952455/148‑1 (MQ=255)
                                                                                                          cTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAgtg                                          >  1:55220/1‑151 (MQ=255)
                                                                                                          cTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAgtg                                          >  1:875590/1‑151 (MQ=255)
                                                                                                            gATCCCCCTTAATGCTTACCGCACGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAgtgt                                         <  2:62121/150‑1 (MQ=255)
                                                                                                             aTCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCg                                       >  2:563618/1‑151 (MQ=255)
                                                                                                               cACCCTTAATCCTTCCCGCACGCTCACCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTc                                        <  1:58188/148‑1 (MQ=255)
                                                                                                                       aTGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAc                             >  2:644960/1‑151 (MQ=255)
                                                                                                                           ttACCGCCCGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcg                            <  2:197518/148‑1 (MQ=255)
                                                                                                                                     cTCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTgg                <  2:221759/150‑1 (MQ=255)
                                                                                                                                     cTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGc               >  1:478067/1‑151 (MQ=255)
                                                                                                                                     cTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGc               >  1:881402/1‑151 (MQ=255)
                                                                                                                                      tCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGc               >  2:868761/1‑150 (MQ=255)
                                                                                                                                       ccccGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGt                  <  2:806977/144‑1 (MQ=255)
                                                                                                                                       cACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACAACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGc               >  1:141767/1‑149 (MQ=255)
                                                                                                                                           gTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGa          >  1:831883/1‑150 (MQ=255)
                                                                                                                                                 cACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTg    <  1:458859/150‑1 (MQ=255)
                                                                                                                                                     aCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGa  <  2:443733/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGA  >  NC_000913/1684601‑1684896

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