Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,825,061	0	G	A	42.3%	28.1 / 13.8	26	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 (12/3);  new base A (0/11);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 4.87e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTG‑ATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAA  >  NC_000913/1824920‑1825202                                                                                                                                                |                                                                                                                                              cATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTG‑ATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgcc                                                                                                                                      <  1:347043/151‑1 (MQ=255)       aaaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTG‑ATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCt                                                                                                                                >  1:253240/1‑151 (MQ=255)         aGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTG‑ATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGt                                                                                                                              >  2:500400/1‑151 (MQ=255)                    cccgcccccgccccTTTTTTTTCTTTTTTTTCTTTCAGTTAACGGTGGTTGTTTCTCGGAAACGTGGTTT‑TTTTACGTTA‑TTGTTACCCGGTTTTTTTTTACTCTTTTTGGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCAt                                                                                                                     <  2:653833/147‑1 (MQ=255)                            cGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTG‑ATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                                                           >  1:670310/1‑151 (MQ=255)                             ggCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTG‑ATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGt                                                                                                          >  1:696299/1‑151 (MQ=255)                                   ttCTTTTAATTTTTCTTTTTAAGGTACCGGTGACTGTTAATGGGCATCTGGGGAG‑ATTAACGTTA‑TGGGATATAAGGTTTTTTTTAAAACTGTTTAAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg                                                                                                     <  1:435066/150‑1 (MQ=255)                                    tCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTG‑ATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAc                                                                                                   >  1:594800/1‑151 (MQ=255)                                               ttATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTG‑ATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcgg                                                                                        >  1:694121/1‑151 (MQ=255)                                               ttATCGTTAAGGTAACGTGGGATGGTTTCCGGCACCTGGGGTG‑ATTTACGTTA‑GTGTTAATCAGTTTTTTTTTACACCTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                                                                                          <  2:822364/149‑1 (MQ=255)                                                         ggTACCGGTGGAGGTTTATGGGAAACTCGGTTA‑AGTTA‑GTTAGGTGTTGATGGGGTTTTTTTTAAAACTGTTGGAAAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTa                                                                                   <  2:270493/146‑1 (MQ=255)                                                                  ggttgtttATCGGCAACTGGGGTT‑TTTTACGTTA‑GTGGTGATCGGTTTTTTTTTACAATTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTgg                                                                         <  1:337321/141‑1 (MQ=255)                                                                   gATGGTTATCGGCAACTCGGGTG‑ATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCt                                                                       >  2:82480/1‑148 (MQ=255)                                                                    aTGGTTATCGGCAACTCGGGTG‑ATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaa                                                                    >  1:350421/1‑150 (MQ=255)                                                                        ttATCGGCAACTGGGGTG‑ATTTACGTTG‑GTGGTGATCAGGTTTTTTTTACAATTTTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat                                                                <  1:576788/150‑1 (MQ=255)                                                                             ggCAACTCGGGTG‑ATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTc                                                           >  1:416373/1‑150 (MQ=255)                                                                                caactgggTGAATTTCCATTA‑GTGCTGCCCCGGCTTTTTTACAAATTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagc                                                        <  2:383413/146‑1 (MQ=255)                                                                                 aCTCGGGTG‑ATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagc                                                        >  1:738824/1‑149 (MQ=255)                                                                                    ggggTG‑TTTTACTTTA‑GTGGTGATCGGTTTTTTTTTACAGTTTTTGGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTg                                                   <  1:596820/150‑1 (MQ=255)                                                                                     gggTG‑TTTTACTTTA‑GTGGTAATCGGGTTTTTTTTACAATTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGt                                                     <  1:230695/148‑1 (MQ=255)                                                                                     gggTG‑ATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTCAAAGATTATTcag                                                         >  2:504090/1‑144 (MQ=255)                                                                                            tttACGTTG‑GTGGTAACCAGGTTTTTTTTACACTTTTTTGAAAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCa                                             <  2:100466/150‑1 (MQ=255)                                                                                                        ggtgATCGGTTTTTTTTTACACCTTTTTAAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAttt                                          <  2:58880/142‑1 (MQ=255)                                                                                                          tattcagttttttttaACAACTTTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTg                                 <  1:500509/140‑1 (MQ=255)                                                                                                             tCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGATTTTTGGCTGAAAGATTATTCAGCAGGTGCAGGCATTTTTCCTTTTGcccc                             >  1:583724/1‑148 (MQ=255)                                                                                                                                       aagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaa  <  2:385213/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                              CATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTG‑ATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAA  >  NC_000913/1824920‑1825202

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