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

CGATAACCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGT  >  NC_000913/1684564‑1684857
                                                                                                                                                     |                                                                                                                                                 
cGATAACCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCAc                                                                                                                                                  <  2:241716/150‑1 (MQ=255)
    aaCCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTa                                                                                                                                             <  1:235273/151‑1 (MQ=255)
             cAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAGTGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTAcc                                                                                                                                    <  1:649971/151‑1 (MQ=255)
                  agcGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACg                                                                                                                               >  2:284073/1‑151 (MQ=255)
                    cGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCt                                                                                                                             >  1:692478/1‑151 (MQ=255)
                     gTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTc                                                                                                                            >  1:65619/1‑151 (MQ=255)
                              tcaTATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTg                                                                                                                   >  1:693125/1‑151 (MQ=255)
                               caTATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGc                                                                                                                  >  1:714629/1‑151 (MQ=255)
                                atatTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGc                                                                                                                  >  1:275118/1‑150 (MQ=255)
                                     cAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGAcc                                                                                                            >  1:270249/1‑151 (MQ=255)
                                      aGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGAccc                                                                                                           <  1:809865/151‑1 (MQ=255)
                                              aaTAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCCCCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCCCCGTTCTGCACGCCCCCGCCACtt                                                                                                   <  2:603552/151‑1 (MQ=255)
                                              aaTAAACGCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGACCTCCCCACCGCGGCACGTGTTAAAGGGGTGCATCGTCATTTGCCGGGGCTGACCACCCTTAATGCTTCCCGCACCCGCCCCGTTCTGCACGTCCCCGCCACtt                                                                                                   <  2:112318/151‑1 (MQ=255)
                                                 aaaCTCAAGCCCGCGTTCGACCATAGCGACGCGGTTAAATCAGAACCCCCCCCCGCGGCACGTTTTGAAGGGGTGCATCGTCATTTGCCGGGCCTCATCCCCCTTAATCCTTCCCGCACGCCCACCGTTCTCCGCGCCCCCGCCACTTTAc                                                                                                <  2:102912/151‑1 (MQ=255)
                                                       aaacccgcgtcagaacATCTCGTCCCGGTTAAATCAAACCCCCCCCCCGCGGCCTTTTTTGAAGGGGTGCATCGTCATTTCCCGGGGTTCACCCCCTTTAATTCTTCCCCCACCCTCCCCTTTCTGCCCGACCCCGCCACTTTACGCttt                                                                                           <  2:711221/135‑1 (MQ=255)
                                                                ttCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCt                                                                                  >  2:232832/1‑150 (MQ=255)
                                                                 tCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCta                                                                                 >  2:572319/1‑150 (MQ=255)
                                                                       atatCGTCACGGTTGAATCAGACCTCCCCCCCGCGGCACGTTTTGAATGGGTGCATCGTCTTTTGCCGGGGCTGATCCCCCTTAATGCTTCCCGCACGCTCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCtaat                                                                               <  1:614416/146‑1 (MQ=255)
                                                                        tatCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACat                                                                         >  1:209571/1‑151 (MQ=255)
                                                                             tCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCAt                                                                    >  1:823781/1‑151 (MQ=255)
                                                                                 gtttaaaTCAAACCCCCCCCCCGCGGCACTTTTTGAATGGGTCCATCGTCATTTGCCGGGGCTGACCCCCCTTAATGCTTCCCGCACGCTCACCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTa                                                                  <  2:851310/144‑1 (MQ=255)
                                                                                 ggTTGAATCAGAACCCCCCCCCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCCCCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTa                                                                  <  1:138032/149‑1 (MQ=255)
                                                                                 ggTTAAATCAAACCCCCCCCCCGCGGCACGTTTTGAAGGGGGGCATCGTCATTTGCCGGGGCTGACCCCCCTTAAGGCTTCCCGCACGCCCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCtaat                                                                               <  2:692401/136‑1 (MQ=255)
                                                                                    taactcaaacccccccccCGCGGCCCTTTTGGAATGGGGCCATCTTCTTTTCCCGGGGCGGATCCCCCTTAATGTTTCCCGCGCGCTCCCCTTTCTGAACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGaa                                                              <  2:714398/138‑1 (MQ=255)
                                                                                    taaaccAGACCCCCCCCCCGCGGCACGTTTTGATTGGGTGCTTCGTCATTTGCCGGGGCTGTTCCCCCTTAATGTTTCCCGCCCGCTCACCTTTCTCCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGaa                                                              <  2:805979/145‑1 (MQ=255)
                                                                                          agacccccccccCGCGGCACTTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCCCCCTTAATGCTTACCGCACGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTa                                                          <  1:118648/142‑1 (MQ=255)
                                                                                          aGACCTCCCCCCCCCGGCAGGTTTTGAATGGGCGCATCGTCTTTTGCCGGGGCTAATCCCCCTTAATGTTTCCCGCACGCCCCCCTTTCCGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACAt                                                       <  2:448524/151‑1 (MQ=255)
                                                                                            acccccccACCGCCGCACGTTGGAAAGGGGGCCGTCGCCTTTGCCGGGGGTTTATCCCCCTTAATGTTTTCCGTACGATCCTCTTTCGTCAGGCCGCCGCCACTTTGCGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTAc                                                         <  2:182202/143‑1 (MQ=255)
                                                                                               tCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtg                                                   >  2:203969/1‑150 (MQ=255)
                                                                                                 ccccccGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCCCCCTTAATGCTTCCCGCCCGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtgatg                                                <  1:371130/151‑1 (MQ=255)
                                                                                                 cccACCGCGGCTTGTTTTGAATGGTGCCATCGTCATTTGCCGGGTTTACCCCCCCTAATTGCTCACGCCGCCCCCCCCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtgatg                                                <  1:562580/151‑1 (MQ=255)
                                                                                                  cccccGCGCCAGTTTTTGAATGGGTGCATCTTCTTTTGCCGGGGCTGACCCCCTTTAATGCTTACCGCCCGCTCACCGTTCTGCCCGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtgatga                                               <  2:97042/148‑1 (MQ=255)
                                                                                                   ccccGCGGCACGTTTTGAATGGGTGCATCGTCTTTTGCCGGGGCTGATCCCCTTTAATGCTTACCGCCCGCTCACCTTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGaa                                              <  2:97308/149‑1 (MQ=255)
                                                                                                     ccGCGGAAGTTTTTGAAGGGGTGCATCTTAATTCGCCGGGGCGGATCCCCCTTAATTCTTACCGCCCGCTCGCCTTCCTCCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAg                                             <  1:648145/150‑1 (MQ=255)
                                                                                                             cGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGcac                                    >  1:540836/1‑151 (MQ=255)
                                                                                                               ttttGATGGGGTTCATCTTCATTTGCCGGGGCTGATCCCCCTTAATGCTTCCCGCACGATCACCTTTCTGCCCAACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCt                                         <  2:227985/144‑1 (MQ=255)
                                                                                                                 ttGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCt                                >  1:528928/1‑151 (MQ=255)
                                                                                                                    aaTGGTTCCATCTTCTTTTCCCGGTCCTAATCCCCCTTAATCCTTCCCCCCCCCTCCCCTTTCTCCCCCACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTg                               <  1:80669/149‑1 (MQ=255)
                                                                                                                    aaTGGGTGCATCTTCATTGGCCGGGGCTGATCCCCTTTAATGCTTCCCGCCCGCTCACCTTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTg                               <  1:718843/149‑1 (MQ=255)
                                                                                                                               cGTCATTTGCGGGGCCTGACCCCCCTTATTGCTTACCGCACCCCCACCGTTCGGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGa                  <  2:275088/151‑1 (MQ=255)
                                                                                                                                gTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGa                  >  1:276149/1‑150 (MQ=255)
                                                                                                                                    tttGCCGGGGCTGATCCCCCTTAATGCTTCCCGCCCGCTCCCCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATc              <  1:494756/150‑1 (MQ=255)
                                                                                                                                    tttGCCGGGGCTGATCCCCCTTAATGCTTCCCGCCCGCTCACCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAAt               <  1:357748/149‑1 (MQ=255)
                                                                                                                                    tttGCCGGGGCTGACCCCCCTTAATGCTTCCCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTg                   <  1:412645/145‑1 (MQ=255)
                                                                                                                                     ttGCCGGGGCTGACCCCCTTTAATGCTTACCGCACGCTCCCCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATc              <  2:65614/149‑1 (MQ=255)
                                                                                                                                      tGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATcacc           >  1:439106/1‑151 (MQ=255)
                                                                                                                                         cGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCGTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACg        >  2:404643/1‑151 (MQ=255)
                                                                                                                                                tGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAgtgt  >  2:823133/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CGATAACCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGT  >  NC_000913/1684564‑1684857

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