Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,825,061	0	G	A	41.2%	61.6 / 22.2	51	intergenic (‑124/+79)	ves/spy	cold‑ and stress‑inducible protein/periplasmic ATP‑independent protein refolding chaperone, stress‑induced
Reads supporting (aligned to +/- strand):  ref base G (20/10);  new base A (0/21);  total (20/31)
Fisher's exact test for biased strand distribution p-value = 7.08e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGC  >  NC_000913/1824917‑1825208                                                                                                                                                  |                                                                                                                                                    aCGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAc                                                                                                                                                <  2:885159/150‑1 (MQ=255) aCGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAc                                                                                                                                                <  2:579460/150‑1 (MQ=255)  cGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccg                                                                                                                                              >  1:755990/1‑151 (MQ=255)         cAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATcc                                                                                                                                       >  2:227498/1‑151 (MQ=255)          aaaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCt                                                                                                                                      >  2:91581/1‑151 (MQ=255)               aTTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTcc                                                                                                                                  >  1:958221/1‑150 (MQ=255)                 tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAcc                                                                                                                               >  2:127600/1‑151 (MQ=255)                    aatctcccccacGGCCATTTCTTTTCATTTTTATCTTTCAGTTACCCTGGT‑TTGTTTATCGCCAACGTTGTTGATTATCATAAGTGAAAACCAGTTTTTTTTCACACCTGTTTAAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCa                                                                                                                            <  2:329575/140‑1 (MQ=255)                    aTACCCGCTCCCGGCCCTTTTTTTTCTTTTTTATCTTTCAGGTACCCGTGG‑TTGTTTATCGGCAACTGGGGTTATTTACTTTAGTGGTGATCAGTTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCg                                                                                                                              <  1:519378/149‑1 (MQ=255)                      aCTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGGGG‑ATGGTTATCGGCAACCCGGGTGATTTACGTTAGTGCTGATCAGGCTTTTTTTACAACTGTTTGAGAAGCGAAAAGACAACCGCCGATCCTGTCCACCGCAtt                                                                                                                          >  2:669472/1‑151 (MQ=255)                          gCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTg                                                                                                                      >  2:954916/1‑151 (MQ=255)                          gCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTg                                                                                                                      >  1:452348/1‑151 (MQ=255)                           cTCCCCGCCCTTTGTTTTCTTTTTTATCGTTCTGCTAATCGTGG‑ATAATTATGAGCATCTGGGGTTATTTATTTTAGTGTTACTCGTTTTTTTTTTAAATCTTTTGAAAAATAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACt                                                                                                                       <  1:513913/149‑1 (MQ=255)                             cccGGCCCTTTTTTTTCATTTTTATCGTTCAGGTAACCGTGG‑TTGGTTACCGGCACCTGGGGTGTTTTACGTTAGTGGTGATCGGTTTTTTTTTACACCTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCa                                                                                                                    <  1:200547/150‑1 (MQ=255)                               cGGCCCGATCTTTTCATTGTTATCGTTCAGGTAACCGTGG‑CTGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGGGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                                                                 >  2:597331/1‑151 (MQ=255)                                ggcctttttttttCATTTTTATCGTTCAGGTAACCGTGG‑TTGTTTTTCGCCACCTCGGGTGATTTACTTTAGTGGTGATCAGGCTTTTTTTACACTTTTTTAAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                                     <  2:679642/137‑1 (MQ=255)                                  cctgTTCTTTTCTTTGTTATCTTTAAGGTAACCGTGG‑ATGGTTTTCGGCAACTGGGGTGATTTACTTTAGTGGTGATCAGGCTTTTTTTACAGTTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTa                                                                                                               <  2:290384/147‑1 (MQ=255)                                   ccGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg                                                                                                              >  2:505717/1‑150 (MQ=255)                                         tttactttGTTATCGTTGAGGTAAGCGTGG‑ATGTTTTTGGGCAACTCGGGTTTTTTCCGTTGGTGTTGATCGGTTTTTTTTTACATCTGTTTGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAc                                                                                                         <  2:717172/143‑1 (MQ=255)                                               tttttACCGTCCAGGTACCCGTTG‑ATGGTTCTGGCCAACGGGGGTGATTTACGTTAGTGGTGCTCAGTTTTTTTTTACACTTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcg                                                                                                  <  2:39381/147‑1 (MQ=255)                                               ttGTTATCTTTCGGGT‑ACCGTGGTGTGTTTATCGGCACCTGGGTTAATTTAGGTTGTTGGTATGCAGGTTTTTTTAACAATTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAcc                                                                                                   <  2:50744/149‑1 (MQ=255)                                                  ttaccgtcaaGGTAACGGTGG‑TGGGTAATCGGAACCGGGGGGTATTTACTTTATTGGTGATCAGGTTTTTTTTACAACTGTTTGGAAAGAGAAAGGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                                                                                                <  2:8644/140‑1 (MQ=255)                                                  ttATCTTTCAGGTACCCGGGG‑ATGTTTATCGGCAACGGGGGTAATTTGCGTTAGTGGTGACCGGGCTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg                                                                                               <  1:927561/150‑1 (MQ=255)                                                    aTCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTc                                                                                            >  2:617431/1‑151 (MQ=255)                                                                acgtgcggtttgtgaTCGGCAACCGGTGTGATTCCCGTTAGTGGTTCTCCGTTTTTTTTCACAACTGTTCGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttttt                                                                                 <  2:470354/137‑1 (MQ=255)                                                                aCCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttttt                                                                                 >  2:121352/1‑150 (MQ=255)                                                                     gg‑ATGTTTATCGGCAACTCGGGTGATTTACGTTGGTGGTGATCAGGTTTTTTTTACAACTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTg                                                                            <  2:962974/150‑1 (MQ=255)                                                                     gg‑ATGTTTACCGCCAACTGGGGGAATTTCCTTTAGTGGTGATCAGGTTTTTTTTACAACTGTTGGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTgg                                                                               <  2:674217/147‑1 (MQ=255)                                                                        aTGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaa                                                                          >  1:16465/1‑150 (MQ=255)                                                                          gtttACCGGAACCTGGGGTGATTTAGGTTAGGGGTGATCAGGTTTTTTTAACACCTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaaa                                                                         <  1:43588/147‑1 (MQ=255)                                                                            ttATCGGCACCTGGGGTTTTTTACGTTAGTGGTGATCGGGTTTTTTTTACAACTGTTGGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaaa                                                                         <  1:816873/147‑1 (MQ=255)                                                                                 gccAACTGGGGTGATTTATGTTAGTGGTGATCAGGTTTTTTTTACACCTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattat                                                                   <  1:456539/146‑1 (MQ=255)                                                                                   aatccggggTGTTTTATGTAAGGGGTGATCGGTTTTTTTTTACACTTTTTTAAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag                                                               <  1:753861/144‑1 (MQ=255)                                                                                    aaCTGGGGTGATTTAGTTTAGTGGTAATCGGGTTTTTTTTCCAACTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagca                                                             <  1:522213/151‑1 (MQ=255)                                                                                         gggTGTTTTACGTTAGTGGTGAGCAGGTTTTTTTTACAGCTTTTTAAGAAAAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCAGAAAGATTATTCAGCAGtt                                                          <  1:907543/149‑1 (MQ=255)                                                                                              aTTTACGTTAGTGGTGAGCGGGTTTTTTTTACAGCTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGc                                                    <  2:104613/150‑1 (MQ=255)                                                                                               tttACGTGGGTGGTGATCAGGTTTTTTTTACAACTTTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAg                                                      <  2:997519/147‑1 (MQ=255)                                                                                                    gTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa                                              >  1:720258/1‑150 (MQ=255)                                                                                                          ggtgTTCCGGTTTTTTTTTCAACTGTTGGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGGGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtt                                          <  1:535497/148‑1 (MQ=255)                                                                                                          ggtgATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTg                                       <  1:101165/151‑1 (MQ=255)                                                                                                          ggtgATAAGTTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt                                         <  1:926375/149‑1 (MQ=255)                                                                                                           gtgATCAGGCTTTTTTTACAACTGTTTGAGAATAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtttt                                        >  1:764363/1‑149 (MQ=255)                                                                                                              aTCAGGCTTTTTTTACAACTGTTTGAGACGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGc                                      >  2:332395/1‑148 (MQ=255)                                                                                                               tCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGGTGGAGGCATTTTACCTTTTGcc                                     >  2:633573/1‑148 (MQ=255)                                                                                                                  gtttttttttACACCTTTGTAAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTgg                                <  2:545678/147‑1 (MQ=255)                                                                                                                   gCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCCGGCATTTTACCTTTTGCCGCTGGAc                              >  1:150776/1‑151 (MQ=255)                                                                                                                         ttaaaaaCTTTTAGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTGTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt                         <  2:607652/145‑1 (MQ=255)                                                                                                                          ttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGt                       >  1:704603/1‑151 (MQ=255)                                                                                                                           tACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGATTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGt                       >  2:258163/1‑150 (MQ=255)                                                                                                                                 tGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACg                 <  1:389782/150‑1 (MQ=255)                                                                                                                                     tGAAAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACg                 <  2:385393/146‑1 (MQ=255)                                                                                                                                            agaaaagaaaACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTTCCTTTTTCCGCTGGACGTTCTGTCAGACGCTTCTCAAAAtt     >  2:812762/1‑151 (MQ=255)                                                                                                                                             aaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaaag      <  2:64402/148‑2 (MQ=255)                                                                                                                                              aaaagaaacCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTTAGCAGTTGCAGGCATTTTACCTTTTGACGCTGGACGTTCTGTCAGACGCTTTTCAAAAtt     >  2:756556/1‑149 (MQ=255)                                                                                                                                                aagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGc  >  1:663686/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                    ACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGC  >  NC_000913/1824917‑1825208

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