Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,222,596	0	G	A	44.2%	50.3 / 37.8	53	F13F (TTC→TTT)	ymgD	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (23/6);  new base A (0/23);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 1.64e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGAC  >  NC_000913/1222449‑1222739                                                                                                                                                       |                                                                                                                                                aaCGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGaaa                                                                                                                                              <  2:119071/151‑1 (MQ=255) aaCGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGaa                                                                                                                                               >  1:596218/1‑150 (MQ=255)    gCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAaaaca                                                                                                                                           <  2:813718/151‑1 (MQ=255)        gCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAaaacaaac                                                                                                                                        >  2:1028949/1‑150 (MQ=255)          attacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAg                                                                                                                                      <  1:696186/150‑1 (MQ=255)             acccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGccc                                                                                                                                  >  1:201110/1‑151 (MQ=255)             acccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGcc                                                                                                                                   <  2:468555/150‑1 (MQ=255)              cccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGccc                                                                                                                                  <  2:979933/150‑1 (MQ=255)               ccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGCCCGc                                                                                                                                <  2:531276/151‑1 (MQ=255)                      cGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTaa                                                                                                                          >  1:45732/1‑150 (MQ=255)                              tttCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAAttt                                                                                                                 >  1:1017594/1‑151 (MQ=255)                                   gCACCCGGTGAGCTGTTAGTTTTTGCCTTTTCGTTGGTTTCTGCAAAATAAATATAGGTTTCAGAGGGAAAC‑CCTTTTTAAATATTGTAGCCCTGCA‑CCAA‑GTTTTTTGGTGAAAAAAAAAACGGGCCCGCCAGTAAAGCAAATTTTTtc                                                                                                             <  1:864594/150‑1 (MQ=255)                                              gcTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCt                                                                                                 >  1:513614/1‑151 (MQ=255)                                                tGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCtt                                                                                                >  1:106396/1‑150 (MQ=255)                                                 gTCAGTTTTTGCCTTTTCGTTGGCTTCTGCGAAAGTAATATAGGTTTCTGACGGAAAC‑CCTTTTTAAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCtt                                                                                                <  2:398439/149‑1 (MQ=255)                                                  ttattttttGCCTTTTCTTGGGCTTCTGCGAAAGTAATAAGGTTTTCTGACGAAAAC‑CCTTTTTTAATATTGTAGCCCTGCC‑CCCG‑GGTAATTGGTGGAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAt                                                                                              <  2:360083/146‑1 (MQ=255)                                                  tCAGTTTTTGCCTTTTCGTTGGTTTCTGCGCAAGTAATATAGGTTTCTGCCGGCAAC‑CCTTTTTTAATATTGTAGTCCGGCG‑CCCA‑GGTCTTTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAt                                                                                              <  2:256387/150‑1 (MQ=255)                                                   cAGTTTTTGCCATTTCGTTGTCTTCTGCGCAATTAATATAGTTTTCGGACGGCAAC‑CCTTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCTTTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAt                                                                                              <  1:668356/149‑1 (MQ=255)                                                    aGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTc                                                                                           >  2:209987/1‑151 (MQ=255)                                                      tttttGCCTTTTCGTTGGTTTCTGCGCAAGAAATATAGGTTTCTGACGGCAAA‑CCTTTTTTAATATTGTAGCCCTGCG‑CCCA‑GGTTTTTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAt                                                                                              <  1:937522/146‑1 (MQ=255)                                                      tttttCCTATTTCGTTGGCTTCCCCGCAAGTAATATCGTTTGCCGACGAAAAC‑CATTTTTTAATATT‑TTGCACCGCGACCCT‑GGTATTTGGTGCAAAAACAAGCAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCa                                                                                          <  2:1029335/150‑1 (MQ=255)                                                        tttGCCATTTCGTTGGTTTCTGCGCAAGTAATATAGTTTTCGGCCGCCAAC‑CTTTTTTAAATTTTTTATCCCGGCC‑CCCA‑GGTTTTTGGTGCAAAAACAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtt                                                                                        <  1:896167/150‑1 (MQ=255)                                                         ttgctttttCGTTGGCTTCTGCGAAAGTAATATAGTTTTCTGACGGCAAC‑CCTTTTTTAATTTTGTAGTCCGGCG‑CCCA‑GGTCATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAttt                                                                                       <  2:530125/144‑1 (MQ=255)                                                            cctttctgTTGTCTTCTGCAAAAAAAATATAGTTTTCTG‑CGGAAAACCCTTTTTTAATTTTGTAGCCCTGCG‑CCCA‑GGTTTTTGGTGCAAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAc                                                                                    <  1:627841/143‑1 (MQ=255)                                                            ccATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAcc                                                                                   >  1:967100/1‑151 (MQ=255)                                                            caaTTTCGTTGTCTTCTGCGAAAGTATTATAGTTTCCCGCCGCCAAC‑CCTTTTTTAAGTTTGAAGTCTGGCG‑CATA‑GGTGATTGGTGCAAAAAAAACCAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAcc                                                                                   <  2:385146/149‑1 (MQ=255)                                                             cATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAccc                                                                                  >  2:107066/1‑151 (MQ=255)                                                               tttCTTTGGTTTCTGCGCAATTAATATGGGTTTCTGACGGCAAT‑CCGTTTTTAATATTGTAGTCCTGCC‑CCCA‑GTTCATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCa                                                                                 <  1:518579/150‑1 (MQ=255)                                                                    ttGGCTTCTGCGCAAGAAATATAGTTTTCGGACGGCAAC‑CCTTTTTAAATATTGTAGTCCTGCC‑CCCG‑GTTCTTTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGa                                                                               <  2:96176/147‑1 (MQ=255)                                                                        cTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGcacc                                                                       >  2:29016/1‑151 (MQ=255)                                                                                 aaGTAATTTGGGTTAAAGCCGGAAAT‑CCTTTTTTATTATTGTCTGCTTGCC‑CCCA‑GGTTTTTGGTGCCAAACAAAATCGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTCCCCAGAATTGcacc                                                                       <  1:94396/142‑1 (MQ=255)                                                                                 aaGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTAcc                                                              >  1:833765/1‑151 (MQ=255)                                                                                 aaGAAATATAGGTTTCGGACGGCAAA‑CTTTTTTTAATATTGTACTCCTGCG‑CCCA‑GGTTTTTGGTCCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCAccc                                                                   <  1:568536/146‑1 (MQ=255)                                                                                           ggTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACg                                                    >  1:63417/1‑151 (MQ=255)                                                                                              ttCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGc                                                  >  1:733260/1‑150 (MQ=255)                                                                                                     ggCAAC‑CCTTTTTTAATATTGTAGTCCGGCG‑CCGA‑GTTTTTTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTaa                                          <  2:181546/151‑1 (MQ=255)                                                                                                          accctttttaAATTTTGTATTCC‑GCC‑CCCATGTTCTTTGGTACAAAAAAAACCAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGc                                      <  1:763297/146‑1 (MQ=255)                                                                                                               ttttaAATATTGTACTCCGGCC‑CCCA‑GGTGATTGGGCAAAAAACAAAAAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGccc                                    <  1:887442/148‑1 (MQ=255)                                                                                                                ttttAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTc                                >  2:329856/1‑151 (MQ=255)                                                                                                                 ttaaatttggTAGCCCGGCC‑CCCG‑GGTCATGGGTGGAAAAAAAACCAGCGCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTc                                <  2:345713/141‑1 (MQ=255)                                                                                                                 ttaaaTTTTGTAGCCCGCCC‑CCCA‑GGTCATGGGTGCAAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTc                                <  2:55939/147‑1 (MQ=255)                                                                                                                  ttAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCg                               >  2:599267/1‑150 (MQ=255)                                                                                                                    aaTATTTTAGCCCCGCC‑CCCA‑GTTCATGGGTGCAAAAACAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAg                             <  2:195646/150‑1 (MQ=255)                                                                                                                      tatTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGcc                          >  2:376889/1‑151 (MQ=255)                                                                                                                            accccggcg‑cCCA‑GTTCTTTGGTCCAAAAAAAAACGGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCAtt                       <  1:720233/142‑1 (MQ=255)                                                                                                                            aGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGcc                    >  2:807888/1‑151 (MQ=255)                                                                                                                                    g‑cCCA‑GGTCATTGGTGCGAAAATAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGa             >  2:810463/1‑150 (MQ=255)                                                                                                                                    g‑cCCA‑GGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGa             >  2:567693/1‑150 (MQ=255)                                                                                                                                       ccaaggttttgggTGCAAAAAAAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAAc           <  2:53014/140‑1 (MQ=255)                                                                                                                                         a‑GGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGcc        >  2:99462/1‑151 (MQ=255)                                                                                                                                              cATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTg    >  1:603640/1‑151 (MQ=255)                                                                                                                                                ttGGTACAAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGa   <  1:744825/150‑1 (MQ=255)                                                                                                                                                 tGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGAc  >  1:452456/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                AACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTTAATATTGTAGTCCTGCG‑CCCA‑GGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGAC  >  NC_000913/1222449‑1222739

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