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

CGATAACCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGT  >  NC_000913/1684564‑1684857
                                                                                                                                                     |                                                                                                                                                 
cGATAACCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCAc                                                                                                                                                  <  2:241716/150‑1 (MQ=255)
    aaCCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTa                                                                                                                                             <  1:235273/151‑1 (MQ=255)
             cAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAGTGGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAATGCTTAcc                                                                                                                                    <  1:649971/151‑1 (MQ=255)
                  agcGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAATGCTTACCGCACg                                                                                                                               >  2:284073/1‑151 (MQ=255)
                    cGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAATGCTTACCGCACGCt                                                                                                                             >  1:692478/1‑151 (MQ=255)
                     gTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAATGCTTACCGCACGCTc                                                                                                                            >  1:65619/1‑151 (MQ=255)
                              tcaTATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTg                                                                                                                   >  1:693125/1‑151 (MQ=255)
                               caTATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGc                                                                                                                  >  1:714629/1‑151 (MQ=255)
                                atatTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGc                                                                                                                  >  1:275118/1‑150 (MQ=255)
                                     cAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGAcc                                                                                                            >  1:270249/1‑151 (MQ=255)
                                      aGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGAccc                                                                                                           <  1:809865/151‑1 (MQ=255)
                                              aaTAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCCCCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAATGCTTACCGCACGCTCCCCGTTCTGCACGCCCCCGCCACtt                                                                                                   <  2:603552/151‑1 (MQ=255)
                                              aaTAAACGCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGACCTCCCCACCGCGGCACGTGTTAAAGGGGTGCATCGTCATTTGCC‑GGGGCTGACCACCCTTAATGCTTCCCGCACCCGCCCCGTTCTGCACGTCCCCGCCACtt                                                                                                   <  2:112318/151‑1 (MQ=255)
                                                 aaaCTCAAGCCCGCGTTCGACCATAGCGACGCGGTTAAATCAGAACCCCCCCCCGCGGCACGTTTTGAAGGGGTGCATCGTCATTTGCC‑GGGCCTCATCCCCCTTAATCCTTCCCGCACGCCCACCGTTCTCCGCGCCCCCGCCACTTTAc                                                                                                <  2:102912/151‑1 (MQ=255)
                                                       aaacccgcgtcagaacATCTCGTCCCGGTTAAATCAAACCCCCCCCCCGCGGCCTTTTTTGAAGGGGTGCATCGTCATTTCCC‑GGGGTTCACCCCCTTTAATTCTTCCCCCACCCTCCCCTTTCTGCCCGACCCCGCCACTTTACGCttt                                                                                           <  2:711221/135‑1 (MQ=255)
                                                                ttCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCt                                                                                  >  2:232832/1‑150 (MQ=255)
                                                                 tCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCta                                                                                 >  2:572319/1‑150 (MQ=255)
                                                                       atatCGTCACGGTTGAATCAGACCTCCCCCCCGCGGCACGTTTTGAATGGGTGCATCGTCTTTTGCC‑GGGGCTGATCCCCCTTAATGCTTCCCGCACGCTCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCtaat                                                                               <  1:614416/146‑1 (MQ=255)
                                                                        tatCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACat                                                                         >  1:209571/1‑151 (MQ=255)
                                                                             tCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCAt                                                                    >  1:823781/1‑151 (MQ=255)
                                                                                 gtttaaaTCAAACCCCCCCCCCGCGGCACTTTTTGAATGGGTCCATCGTCATTTGCC‑GGGGCTGACCCCCCTTAATGCTTCCCGCACGCTCACCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTa                                                                  <  2:851310/144‑1 (MQ=255)
                                                                                 ggTTGAATCAGAACCCCCCCCCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCCCCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTa                                                                  <  1:138032/149‑1 (MQ=255)
                                                                                 ggTTAAATCAAACCCCCCCCCCGCGGCACGTTTTGAAGGGGGGCATCGTCATTTGCC‑GGGGCTGACCCCCCTTAAGGCTTCCCGCACGCCCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCtaat                                                                               <  2:692401/136‑1 (MQ=255)
                                                                                    taactcaaacccccccccCGCGGCCCTTTTGGAATGGGGCCATCTTCTTTTCCC‑GGGGCGGATCCCCCTTAATGTTTCCCGCGCGCTCCCCTTTCTGAACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGaa                                                              <  2:714398/138‑1 (MQ=255)
                                                                                    taaaccAGACCCCCCCCCCGCGGCACGTTTTGATTGGGTGCTTCGTCATTTGCC‑GGGGCTGTTCCCCCTTAATGTTTCCCGCCCGCTCACCTTTCTCCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGaa                                                              <  2:805979/145‑1 (MQ=255)
                                                                                          agacccccccccCGCGGCACTTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCCCCCTTAATGCTTACCGCACGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTa                                                          <  1:118648/142‑1 (MQ=255)
                                                                                          aGACCTCCCCCCCCCGGCAGGTTTTGAATGGGCGCATCGTCTTTTGCC‑GGGGCTAATCCCCCTTAATGTTTCCCGCACGCCCCCCTTTCCGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACAt                                                       <  2:448524/151‑1 (MQ=255)
                                                                                            acccccccACCGCCGCACGTTGGAAAGGGGGCCGTCG‑CCTTTGCCGGGGGTTTATCCCCCTTAATGTTTTCCGTACGATCCTCTTTCGTCAGGCCGCCGCCACTTTGCGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTAc                                                         <  2:182202/143‑1 (MQ=255)
                                                                                               tCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtg                                                   >  2:203969/1‑150 (MQ=255)
                                                                                                 ccccccGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCCCCCTTAATGCTTCCCGCCCGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtgatg                                                <  1:371130/151‑1 (MQ=255)
                                                                                                 cccACCGCGGCTTGTTTTGAATGGTGCCATCGTCATTTGCC‑GGGTTTACCCCCCCTAATTGCTCACGCCGCCCCCCCCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtgatg                                                <  1:562580/151‑1 (MQ=255)
                                                                                                  cccccGCGCCAGTTTTTGAATGGGTGCATCTTCTTTTGCC‑GGGGCTGACCCCCTTTAATGCTTACCGCCCGCTCACCGTTCTGCCCGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtgatga                                               <  2:97042/148‑1 (MQ=255)
                                                                                                   ccccGCGGCACGTTTTGAATGGGTGCATCGTCTTTTGCC‑GGGGCTGATCCCCTTTAATGCTTACCGCCCGCTCACCTTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGaa                                              <  2:97308/149‑1 (MQ=255)
                                                                                                     ccGCGGAAGTTTTTGAAGGGGTGCATCTTAATTCGCC‑GGGGCGGATCCCCCTTAATTCTTACCGCCCGCTCGCCTTCCTCCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAg                                             <  1:648145/150‑1 (MQ=255)
                                                                                                             cGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGcac                                    >  1:540836/1‑151 (MQ=255)
                                                                                                               ttttGATGGGGTTCATCTTCATTTGCC‑GGGGCTGATCCCCCTTAATGCTTCCCGCACGATCACCTTTCTGCCCAACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCt                                         <  2:227985/144‑1 (MQ=255)
                                                                                                                 ttGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCt                                >  1:528928/1‑151 (MQ=255)
                                                                                                                    aaTGGTTCCATCTTCTTTTCCC‑GGTCCTAATCCCCCTTAATCCTTCCCCCCCCCTCCCCTTTCTCCCCCACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTg                               <  1:80669/149‑1 (MQ=255)
                                                                                                                    aaTGGGTGCATCTTCATTGGCC‑GGGGCTGATCCCCTTTAATGCTTCCCGCCCGCTCACCTTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTg                               <  1:718843/149‑1 (MQ=255)
                                                                                                                               cGTCATTTGCG‑GGGCCTGACCCCCCTTATTGCTTACCGCACCCCCACCGTTCGGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGa                  <  2:275088/151‑1 (MQ=255)
                                                                                                                                gTCATTTGCC‑GGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGa                  >  1:276149/1‑150 (MQ=255)
                                                                                                                                    tttGCC‑GGGGCTGATCCCCCTTAATGCTTCCCGCCCGCTCCCCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATc              <  1:494756/150‑1 (MQ=255)
                                                                                                                                    tttGCC‑GGGGCTGATCCCCCTTAATGCTTCCCGCCCGCTCACCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAAt               <  1:357748/149‑1 (MQ=255)
                                                                                                                                    tttGCC‑GGGGCTGACCCCCCTTAATGCTTCCCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTg                   <  1:412645/145‑1 (MQ=255)
                                                                                                                                     ttGCC‑GGGGCTGACCCCCTTTAATGCTTACCGCACGCTCCCCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATc              <  2:65614/149‑1 (MQ=255)
                                                                                                                                      tGCC‑GGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATcacc           >  1:439106/1‑151 (MQ=255)
                                                                                                                                         c‑GGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCGTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACg        >  2:404643/1‑151 (MQ=255)
                                                                                                                                                tGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAgtgt  >  2:823133/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CGATAACCTTTCGCAAGCAGCGTTGGTCAGTCATATTCAGCACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCC‑GGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGT  >  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: ‑