Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,684,7490AC57.7% 26.7 / 39.2 26D164A (GAC→GCC) tusinhibitor of replication at Ter, DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (11/0);  new base C (0/15);  total (11/15)
Fisher's exact test for biased strand distribution p-value = 1.29e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-06

AAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTG  >  NC_000913/1684613‑1684881
                                                                                                                                        |                                                                                                                                    
aaaCTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTAc                                                                                                                        >  2:341211/1‑151 (MQ=255)
       aGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTttgg                                                                                                                 >  2:66241/1‑151 (MQ=255)
                 cGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCtaat                                                                                                       >  1:659066/1‑151 (MQ=255)
                       tatCGCCACGGTTAAATCAGACCCCCCCCCCGCGGCACGTTTTGAAGGGGTGCATCGTCTTTGGCGGGGGCTGATCCCCCTTAATGCTTCCCGCACGCTCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCtaata                                                                                                      <  1:400413/146‑1 (MQ=255)
                             cgcgtttaaaTCAGACCCCCCCCCCGGGGCAGGTTTTGAAGGGGTCCTTCTTCATTTGCCGGGGCGGACCGCCCTTAATGCTTCGGGCGCGCTAACCGTTCTGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACat                                                                                                 <  2:594536/137‑1 (MQ=255)
                             cACGGTTGAATCAAACCCCCCCCCCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCCCCCTTAATGCTTACCGCCCGCCCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATc                                                                                              <  1:614244/148‑1 (MQ=255)
                                  tttaaTAAGAACTCCCCCCCGCGCCCCGTTTTGAAGGGGTGCATCGTAATTTGCCGGGGCTGATCCCCCTTAATTTTTACCGCCCGCCCACCGTTCTGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGa                                                                                       <  1:502147/147‑1 (MQ=255)
                                   taaaccaaaaCCCCCCCCCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCCCCCTTAATGCTCCCCCCCCCCTCACCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGaa                                                                                      <  2:781753/142‑1 (MQ=255)
                                     aaTCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAttt                                                                                   >  1:237274/1‑151 (MQ=255)
                                         agacccccccccCGCGTCCCGTTTTGAATGGGTGCACCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTCCCGCCCGCTCACCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACAt                                                                               <  1:263884/145‑1 (MQ=255)
                                            ccgccccccCGCGCCCTGTTTTGATTGGGTGCATCTTCTTTTGCCGGGGCTGACCGCCCTTATTGTTTACCGCGCGCCCACCGTTCTTCACACCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGt                                                                            <  2:177314/148‑1 (MQ=255)
                                             ccccccccCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCCCCTTTAATGCTTACCGCACGCTCACCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGt                                                                            <  1:169868/148‑1 (MQ=255)
                                                          ccgtttttAAAGGGGTGCATGGTCTTTTGCCGGGGCTAACCCCCCTTATTGCTTCCCGCGCGCTCCCCTTCCGGCACACCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGc                                                              <  1:750618/147‑1 (MQ=255)
                                                             gTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGcac                                                            >  1:306098/1‑150 (MQ=255)
                                                              ttttGAATGGGCGAATCGTCTTTGGCCGGGGCTGATCCCCTTTAATGTTTACCGCCCGCCCACCGTTCTCCACGCCACCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGca                                                             <  2:237329/148‑1 (MQ=255)
                                                                   aaTGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTgg                                                      >  2:179120/1‑150 (MQ=255)
                                                                                    ttGCCGGGGCTGACCCCCCTTAATGCTTCCCGCCCGCTCACCGTTCTGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATca                                     <  1:788637/150‑1 (MQ=255)
                                                                                         ggggCTGCCCCCCTTTAATGCTTCCCCCCCCCTCCCCGTTCTCCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGc                               <  2:328387/151‑1 (MQ=255)
                                                                                         ggggCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGc                               >  1:666928/1‑151 (MQ=255)
                                                                                          gggCTGATCCCCTTTAATGCTTACCGCACGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCa                              <  1:792057/151‑1 (MQ=255)
                                                                                              tGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAgtgt                          >  1:217541/1‑151 (MQ=255)
                                                                                                             cTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcgcg           >  1:517737/1‑151 (MQ=255)
                                                                                                                aCCGCACGCTCACCGTTCTGCACGACCCCGCCCCTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcgcgcg         >  1:775099/1‑150 (MQ=255)
                                                                                                                 ccGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcgcgcg         >  1:244971/1‑149 (MQ=255)
                                                                                                                          ccccGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCgagg      <  2:196224/141‑1 (MQ=255)
                                                                                                                          cacgtttttccACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTg  <  2:59360/137‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                    
AAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTG  >  NC_000913/1684613‑1684881

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