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

ACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCC  >  NC_000913/1222450‑1222733                                                                                                                                                   |                                                                                                                                          aCGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGaaa                                                                                                                                        >  1:248110/1‑150 (MQ=255)   gCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAaaac                                                                                                                                      >  2:31757/1‑150 (MQ=255)               cattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGcc                                                                                                                         >  2:792802/1‑151 (MQ=255)                    ccGCAACAATTTCTGCAACCTGTGCGCTGTCATTTTTTGCCTTTTCGTTGGCTTCTGCGAAAGTAATATAGGTTTCTGACGGCAACCCTTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCAAAAACAAACAGGCCCGCCAGTa                                                                                                                     <  2:345791/150‑1 (MQ=255)                         acatTTTCTGCAACCTGTGCGCTGTCAGTTTTTCCCATTTCGTTGGCTTCTGCGCAAGAAATATAGGTTTCGGACGGCAAACCTTTTTTAATATTGTAGTCCGGCGCCCAGTTCTTTGGTGCAAAAACAAACAGGCCCGCCAGTAAAg                                                                                                                  <  1:602762/148‑1 (MQ=255)                             tttCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAAttt                                                                                                           >  1:906987/1‑151 (MQ=255)                                   cacccTGTGCGCTGTCGTTTTTCCCCTTTTCGTTGGTTTCTGCAAAAGTAATAGAGGTTTGCGAGGACAACCCTTTTTTAATATTGTAGCCCTGCGCCCATGTCTTTGGGGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtc                                                                                                       <  1:220981/146‑1 (MQ=255)                                    aaCCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtcatc                                                                                                    >  2:277533/1‑151 (MQ=255)                                         gttcgctttattttttgGCAATTTCGTTGGTTTCTGCGCAAGTAATATAGGTTTCTCACGGCACACCGTTTTAAATATGGTAGCCCGGCGCCCAGTTCTTTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtca                                                                                                      <  1:809664/132‑1 (MQ=255)                                         gtgCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAACCCTTTTTTAATATTGTAGTCCTGCCCCCAGGTCATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtcatcatca                                                                                                <  2:39391/150‑1 (MQ=255)                                         gggcgcTGTAATTTTTTGCCTTTTCGTGGTTTTCTCCGCAAGTAATATAGGTTTCTGCCGGAAACCCTTTTTAAATTTTGTAGCCCGCCGCCCAGTTAATTGGTGAAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTcatcatcat                                                                                               <  2:809820/149‑1 (MQ=255)                                             gcTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCt                                                                                           >  1:498547/1‑151 (MQ=255)                                                     tttttGCCATTTCGTTGGTTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTAAATATTGTAGCCCTGCCCCCAGTTCATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCtt                                                                                          <  1:312907/144‑1 (MQ=255)                                                      ttttGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtt                                                                                  >  1:720940/1‑151 (MQ=255)                                                       tttCCCTTTTCTTTGTCTTCTGCGAAATTATTATGGTTTTCTGACGGAAACCTTTTTTTAATATTGTGGTCCGGCCCCCAGGTTTTTGGTGCAAAAACAACCAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtt                                                                                  <  1:556079/150‑1 (MQ=255)                                                           ctttttGGTTGGTTTCGCCGAAAGTAATATAGTTTTCTGACGGAAACCCTTTTTAAATATTGTAGTCCTCCCCCCAGGTTTTTGGTGAAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTa                                                                               <  1:811608/146‑1 (MQ=255)                                                              tttCGTTGGTTTCTGCGCAAGTAATATAGGTTTCTGACGGCAACCCTTTTTTAATATTGTAGTCCTGCGCCCAGGTCTTTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAccc                                                                            <  2:591415/149‑1 (MQ=255)                                                                   ttGGTTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCTTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAAt                                                                       <  1:419018/149‑1 (MQ=255)                                                                   ttGGTTTCTGCGCAAGTAATATAGGTTTCTGACGGAAACCCGTTTTTAATTTTGTAGCCCTGCGCCCAGGTTTTTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAAtt                                                                      <  1:475708/150‑1 (MQ=255)                                                                    tGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGAAAACCTTTTTTAATATTGTAGTCCTGCCCCCAGGTCATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTg                                                                     <  2:183214/150‑1 (MQ=255)                                                                     ggtttCTGCGCAAGTAATATAGGTTTCTGACGGCAACCCTTTTTAAATATTGTAGTCCTGCGCCCAGTTTTTTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGca                                                                   <  1:387175/148‑1 (MQ=255)                                                                      gCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGcac                                                                  >  1:511831/1‑151 (MQ=255)                                                                         tCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGcaccac                                                               >  1:463040/1‑151 (MQ=255)                                                                            gcgcAAGTAATATAGTTTTCTGACGGCAACCCTTTTTTAATATTGTTGTCCTGCCCCCAGGTCTTTGGTGCAAAAACAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCg                                                            <  2:372229/151‑1 (MQ=255)                                                                             cgcAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGt                                                           >  2:719928/1‑151 (MQ=255)                                                                                aaGAAATATTGTTTTCTGACGTTAACCTTTTTTTAATTTTGTATCCCGGCGCCTGTGTCTTTGGGGAAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCaccacc                                                              <  2:423407/145‑1 (MQ=255)                                                                                   tAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGcc                                                     >  1:203977/1‑151 (MQ=255)                                                                                    aaTATAGGTTTCTGACGGCAACCCTTTTTTAATATTGTAGTCCTGCGCCCAGGTCTTTGGTGAAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGcc                                                     <  1:914371/150‑1 (MQ=255)                                                                                    aaTATAGGTTTCTGACGGCAAACCTTTTTTAATATTGTAGCCCTGCGCCCAGGTCATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGc                                                      <  2:314186/149‑1 (MQ=255)                                                                                      tataGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAAt                                                  >  2:269393/1‑151 (MQ=255)                                                                                         aGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCAc                                               >  2:899889/1‑151 (MQ=255)                                                                                          ggTTTCTGACGGCAAACCTTTTTTAATATTGTATCCCTGCCCCCAGTTCATTGGTGCAAAAAAAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATc                                                 <  2:930178/148‑1 (MQ=255)                                                                                           gTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACgg                                             >  2:621043/1‑151 (MQ=255)                                                                                                tGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCt                                        >  1:15977/1‑151 (MQ=255)                                                                                                tGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCt                                        >  1:305947/1‑151 (MQ=255)                                                                                                 gACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCtt                                       >  2:817138/1‑151 (MQ=255)                                                                                                       aacccttttttAATATTGTAGCCCGGCGCCCAGGTATTTGGGGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATc                                  <  2:528407/144‑1 (MQ=255)                                                                                                       aacccttttttAATATTGTAGCCCGGCCCCCAGGTAATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATc                                  <  2:824142/144‑1 (MQ=255)                                                                                                       aaccctttttaAATATTGTGGCCCGGCCCCCAGTTATTGGGTGCAAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCg                                 <  2:248138/145‑1 (MQ=255)                                                                                                       aaccctttttaAATATTGTAGCCCTGCCCCCAGGTCATTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCg                                 <  2:552131/145‑1 (MQ=255)                                                                                                       aaaCCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCg                                 >  2:305192/1‑151 (MQ=255)                                                                                                        aaCCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGc                                >  2:799075/1‑151 (MQ=255)                                                                                                         aCCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGcc                               >  2:904715/1‑151 (MQ=255)                                                                                                            gTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCt                            >  1:394496/1‑151 (MQ=255)                                                                                                             tttttAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCtt                           >  1:826563/1‑151 (MQ=255)                                                                                                              ttttAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTc                          >  1:576866/1‑151 (MQ=255)                                                                                                                   atatTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGc                     >  2:641006/1‑151 (MQ=255)                                                                                                                   atatTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGc                     >  2:469505/1‑151 (MQ=255)                                                                                                                         tAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGc               >  1:636560/1‑151 (MQ=255)                                                                                                                           gTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCg             >  1:102659/1‑151 (MQ=255)                                                                                                                            tCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCgg            >  2:794903/1‑151 (MQ=255)                                                                                                                             ccTGCCCCGAGGTCATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCg             <  1:595501/149‑1 (MQ=255)                                                                                                                                  gcCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGa       >  2:128507/1‑150 (MQ=255)                                                                                                                                   cccGGGTCTTTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGaa      <  2:826925/150‑1 (MQ=255)                                                                                                                                   cccAGGTCATTGGTGAAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGaa      <  2:125290/150‑1 (MQ=255)                                                                                                                                   cccAGGTCATTGGGGCAAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGaa      <  2:64433/150‑1 (MQ=255)                                                                                                                                        gTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGcc  >  1:506324/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                          ACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCC  >  NC_000913/1222450‑1222733

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