Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,684,749	0	A	C	41.5%	84.4 / 11.0	53	D164A (GAC→GCC)	tus	inhibitor of replication at Ter, DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (23/8);  new base C (0/22);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 1.43e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGG  >  NC_000913/1684605‑1684895                                                                                                                                                  |                                                                                                                                                   acaTCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGc                                                                                                                                              <  1:564/151‑1 (MQ=255)      aaTAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACtt                                                                                                                                         <  2:511357/151‑1 (MQ=255)            cTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCCGGGGCTGATCACACTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCt                                                                                                                                   >  1:494826/1‑151 (MQ=255)                 gACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTttggt                                                                                                                              >  2:839281/1‑151 (MQ=255)                 gACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTttgg                                                                                                                               >  1:361133/1‑150 (MQ=255)                   cccctttCGAGCATATCGTCACGGTTGAATCAGACCCCCCCCCCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCCGGGGCTGATCCCCTTTAATGCTTCCCGCCCGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTttggtt                                                                                                                             <  2:715628/145‑1 (MQ=255)                        ttCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCt                                                                                                                        >  1:54727/1‑150 (MQ=255)                           gaccATACCTTCAGGGTTAAATCAGAACCCCCCCCCGCGGCACGTTTTAAAG‑GGGTGCATCGTATTTGGCGGGGGCTGATCCCCCTTAATGCTTACCCCCCGCCCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCtaat                                                                                                                     <  2:61938/147‑1 (MQ=255)                           gAGCATATCGTCACGGTTGAATCAGAACCCCCCCCCGCGGCACGTTTTGAAT‑GGGTGCATCGTCTTTTGCCGGGGCTGATCCCCCTTAATGCTTACCGCACGCTCCCCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCtaat                                                                                                                     <  2:468490/150‑1 (MQ=255)                               atatCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAAc                                                                                                                 >  2:890869/1‑150 (MQ=255)                                    gTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATc                                                                                                            >  2:226298/1‑150 (MQ=255)                                         ggTTGAATCAGACCCCCCCCCCGCGGCACGTTTTGAAG‑GGGTGCATCGTCATTTGCCGGGGCTGATCCCCCTTAATGCTTCCCGCACGCCCCCCTTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAg                                                                                                      <  2:585823/151‑1 (MQ=255)                                              aaccaaacccccccccCGCGCCCCTTTTTGAGG‑GGGTCCTCCTTCTTTTGCGGGGCCTAACCCCCCTTAATCTTTCCCCCTCGCTCCCCTTTCTCCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTaa                                                                                                       <  1:575612/135‑1 (MQ=255)                                              aaTCAGACCCCCCCCCCGCGCCACTTTTTGAAG‑GGGTGCTTCGTCTTTTGCCGGGGCTGATCCCCCTTAATGCTTCCCGCCCGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGaa                                                                                                    <  2:1156660/148‑1 (MQ=255)                                              aaTCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAttt                                                                                                 >  2:388808/1‑151 (MQ=255)                                                acAGACCCCCCCCCCCCGGCATGTTTTGAAT‑GGGTCCATCGTCATTGCCCGGGGCTAATCCCCCTAAATGTTTCCCGCACGCCCACCATTCGGCTAGCCCCCGCCACTTGACGCTTTGGTTGGGCTAATAAACATATCATTAAGaa                                                                                                    <  1:1047322/145‑1 (MQ=255)                                                      ccccccACCGCGCCACGTTTTGAAG‑GGGTGCTTCTTCTTTTGCCGGGGCTGACCCCCCTAAATGTTTCCCGCACGCTCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtg                                                                                         <  2:67675/149‑1 (MQ=255)                                                        ccccACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtgat                                                                                       >  1:51639/1‑151 (MQ=255)                                                         cccACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGtgatg                                                                                      >  2:437668/1‑151 (MQ=255)                                                            cccGCGCCACGTTTAAAAT‑GGGTGCATCGTCATTTGCTGGGGCTGCCCCCCCTTAATCCTTACGGCTCGCACCCCGTTCTCCACCCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAg                                                                                   <  1:311478/150‑1 (MQ=255)                                                                  gccatttttGAAAGGGGTGCTTCTTCTTTTGCCGGGGCTGATCCCCTTTAATGCTTCCCGCACGCTCCCCGTGCTCCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCt                                                                               <  1:951626/145‑1 (MQ=255)                                                                    aCGTTTTGAAT‑GGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGc                                                                            >  2:1365334/1‑150 (MQ=255)                                                                         ttGAAT‑GGGCGCATTTTCTTTTGCCGGGCCTGATCCCCCTTAATTCTTCCCGCACGCTCCCCTTTCCCCCCACCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGc                                                                       <  1:432788/150‑1 (MQ=255)                                                                            aaT‑GGTTGCATCGTCTTTTGCCGGGGCTGACCCCCCTTAATGCTTCCCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGa                                                                   <  2:239093/151‑1 (MQ=255)                                                                            aaT‑GGGTCCTTCGTCTTTGGCGGGGGCGGACCCCCTTTAATGCTCCCCTCGCGCCCACTTTTCTCCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGc                                                                       <  1:68122/147‑1 (MQ=255)                                                                             aT‑GGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGa                                                                   >  1:1404069/1‑150 (MQ=255)                                                                                gggTGCATCGTCTTTTGCCGGGGCTGATCCCCTTTATTGCTTCCCGCACGCTCCCCTTTCTCCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGaaa                                                                 <  1:117761/150‑1 (MQ=255)                                                                                 ggccaaccGTCATTCCGCGGGCCTGCTCCCCCTTATTGCTACCCGCCCGCCCCCATTTCTGCCCCCCGCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGaaa                                                                 <  1:652018/142‑1 (MQ=255)                                                                                 ggTGCTTCGTCTTTGGCCGGGCCTGTTCCCCCTTAATGTTTCCCGCACTCTCCCCTTTCTGCACGCCTCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGaaaa                                                                <  1:964805/150‑1 (MQ=255)                                                                                    gCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGc                                                            >  1:50648/1‑151 (MQ=255)                                                                                      aTCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCt                                                          >  2:292400/1‑151 (MQ=255)                                                                                              ttGCCGGGGCTGATCCCCTTTAATGCTTCCCGCCCGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATca                                                   <  2:1254236/150‑1 (MQ=255)                                                                                                      gcggATCCCCCTTAATACTTCCGCCCCGCACGCGTTTCGGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATcacca                                                <  1:781796/142‑1 (MQ=255)                                                                                                        tattcccccTTAATGCTCCCCGCCCGCTCACCTTTCGGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAg                                           <  2:1360791/142‑1 (MQ=255)                                                                                                        tGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAgtgt                                        >  2:876917/1‑151 (MQ=255)                                                                                                            cccccTTAATGCTTCCCGCACGCTCACCTTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTc                                       <  1:472416/146‑1 (MQ=255)                                                                                                                 ttAATCCTTCCCGCACTCTCACCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCg                                 <  1:284278/149‑1 (MQ=255)                                                                                                                  tAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTgg                              >  1:163604/1‑151 (MQ=255)                                                                                                                   aaTTCTTCCCGCCCGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTg                               <  2:768121/149‑1 (MQ=255)                                                                                                                     tGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcg                           >  1:403933/1‑151 (MQ=255)                                                                                                                       cTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcgcg                         >  1:1209563/1‑151 (MQ=255)                                                                                                                           ccgcccgcTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCga                      <  1:386222/150‑1 (MQ=255)                                                                                                                            cGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAg                                                             >  1:312426/1‑110 (MQ=255)                                                                                                                                cGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTg                >  1:1144132/1‑151 (MQ=255)                                                                                                                                  cTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTgg               <  1:208755/150‑1 (MQ=255)                                                                                                                                  cTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGc              >  2:488477/1‑151 (MQ=255)                                                                                                                                    ccccGTTCTGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCa             <  2:494181/148‑1 (MQ=255)                                                                                                                                    cACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCa             >  2:963711/1‑150 (MQ=255)                                                                                                                                     ccctttCTGCCCGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCaaa           <  2:754569/147‑1 (MQ=255)                                                                                                                                     aCCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCaa            >  2:1209915/1‑150 (MQ=255)                                                                                                                                           cTGCCCGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGaaa       <  2:1225862/149‑1 (MQ=255)                                                                                                                                              cACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAAc     <  2:54733/148‑1 (MQ=255)                                                                                                                                               aCGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTgg  >  1:100039/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                   ACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAAT‑GGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGG  >  NC_000913/1684605‑1684895

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