Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	803,869	0	G	T	46.7%	63.4 / 24.1	45	V123L (GTG→TTG)	ybhJ	aconitase family protein
Reads supporting (aligned to +/- strand):  ref base G (4/20);  new base T (21/0);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 4.95e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGAT‑GGCAGGCGGC‑GGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑GTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAA  >  NC_000913/803720‑804013                                                                                                                                                      |                                                                                                                                                    aGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTtgt                                                                                                                                                    <  2:704993/151‑1 (MQ=255)  gACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTtgtg                                                                                                                                                   <  1:296610/151‑1 (MQ=255)  gACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGGTATGGCGGGCTTTTtttg                                                                                                                                                   >  1:630306/1‑151 (MQ=255)          aaGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCa                                                                                                                                           <  1:548804/151‑1 (MQ=255)            gCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTcc                                                                                                                                              >  1:931920/1‑146 (MQ=255)             cGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGGTTTTtttgtccccacgt                                                                                                                                          >  1:639855/1‑140 (MQ=255)                ccGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCGTAtt                                                                                                                                       >  1:401601/1‑149 (MQ=255)                          aCGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTcatccat                                                                                                                           <  1:262059/151‑1 (MQ=255)                           cGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTcatccatc                                                                                                                          <  2:444484/151‑1 (MQ=255)                                cccGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCATATTGCGGGCATCCATCAata                                                                                                                      >  1:524176/1‑150 (MQ=255)                                    cTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGCCTCCCCCTATTTCGGGCATCCATCCatatat                                                                                                                   >  1:40151/1‑149 (MQ=255)                                       ccGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTg                                                                                                              <  2:927843/151‑1 (MQ=255)                                       ccGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTg                                                                                                              <  1:611323/151‑1 (MQ=255)                                        cGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCCCCCCCTATTTGGGGCCTCCCACAATATATGCg                                                                                                                >  1:525177/1‑148 (MQ=255)                                           aTGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAgat                                                                                                          <  2:49447/151‑1 (MQ=255)                                              tgCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGGCCCCCCCTATTTGGGTCAACCCTCAATATATGCggggg                                                                                                            >  1:316205/1‑142 (MQ=255)                                              tgCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGGCCCCGCCTATTTCGGGCATCCATCAATATATGCg                                                                                                                >  1:422948/1‑142 (MQ=255)                                               gCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCCCCCCATATTGCGGTCATCCCCTAATATATGCGTGGgatg                                                                                                         >  1:317854/1‑148 (MQ=255)                                               gCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAgatgat                                                                                                       <  2:165561/150‑1 (MQ=255)                                                  gACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGAT‑GGc                                                                                                   <  1:345694/150‑1 (MQ=255)                                                  gACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGAT‑GGc                                                                                                   <  1:962689/150‑1 (MQ=255)                                                    ccAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGAT‑GGCAgg                                                                                                <  2:12885/151‑1 (MQ=255)                                                            cATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGTCTCCGCCTATTTCGGGCATTCATCAATATATGCGTGAGAAGAT‑GGCAggcgggg                                                                                           >  2:925881/1‑146 (MQ=255)                                                               aaCTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGAT‑GGCAGGCGGC‑GGCaaaa                                                                                    <  2:72035/151‑1 (MQ=255)                                                                actcactcTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGAT‑GGCAGGCGGC‑GGCAAAAt                                                                                   <  1:138157/151‑1 (MQ=255)                                                                  tcactcTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGGATTTTTTGTCTCCCCCTATTTCGGGCAACCCTCAATATATTGGTGAGAAGGT‑TGCAGGCGGG‑GGGaaaa                                                                                    >  1:663530/1‑148 (MQ=255)                                                                     ctctGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGAT‑GGCAGGCGGC‑GGCAAAATGATc                                                                               <  2:90182/150‑1 (MQ=255)                                                                          cgcCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGCCTCCCCCTATTGCGGGCATCCATCAATATATGCGGGGGGTGAT‑GGGAGGGGGCGGGCaaaaaataccccgg                                                                          >  2:655510/1‑141 (MQ=255)                                                                           gcCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCATATTTCGGGCATACCTCCATATATGCGTGGGGTGATGGGCGGGGGGG‑GGCAAAA‑GATTCTCggg                                                                         >  2:890969/1‑150 (MQ=255)                                                                            ccGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGTCTCCGCATATTTCGGGCATCCATCAATATATGCGTGGGATGAT‑GGCAGGCGGC‑GGCGAAAAGAACCTCggg                                                                         >  2:738742/1‑149 (MQ=255)                                                                            ccGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGAT‑GGCAGGCGGC‑GGCAAAATGATCCTCGGGTc                                                                       <  1:845796/151‑1 (MQ=255)                                                                             cGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGTCTCCGCCTATTGCGGGCATTCCTCAATATATGCGTGAGGTGAT‑GGGCGGCGGC‑GGGCAAAAAATCCCCGGGtg                                                                       >  2:663758/1‑149 (MQ=255)                                                                             cGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCCAATTTCGGGCATCCCTCAATAAAAGCGTGGGGTGGT‑GGGAGGGGGG‑GGGCAAAAGATTCTCggg                                                                         >  1:442118/1‑148 (MQ=255)                                                                                cggcggCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGAT‑GGCAGGCGGC‑GGCAAAATGATCCTCGGGTCAGAc                                                                   <  2:818645/151‑1 (MQ=255)                                                                                   cggcACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCCCCCCATATTTCGGGCATCCATCAAAATATGCGTGAGATGAT‑GGCAGGCGGC‑GGGAAAATGATCCTCCGGGCAGACa                                                                  >  1:747464/1‑149 (MQ=255)                                                                                         tATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTTTCTCCCCCTATTTCGGGCATCCCACAATATATGCGTGGGGTGAT‑GGCAGGGGGG‑GGCAAAAAGAATCTCGGGGCAGAAAGCCCCAcc                                                          >  2:114643/1‑151 (MQ=255)                                                                                               cGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGAT‑GGCAGGCGGC‑GGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑GTTAc                                                   <  2:816552/151‑1 (MQ=255)                                                                                                ggTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGAT‑GGCAGGCGGC‑GGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑GTTACg                                                  <  2:924315/151‑1 (MQ=255)                                                                                                 gTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGTCTCCGCCTATTGCGGGTATCCCACAAAATATGGGGGGGATGAA‑GGCAGGCGGC‑GGCAAAATGATCCTTGGGTCAGACAGCCACACCC‑GGTAggg                                                 >  1:879976/1‑151 (MQ=255)                                                                                                        ccATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCCTATTTCGGGCATCCCTCAATATATGCGTGAGATGAT‑GGCAGGCGGC‑GGCAAAATGATCCTCGGGGCAGACAGCCACACCC‑CTTACCGGGCAtt                                           >  1:535148/1‑150 (MQ=255)                                                                                                        ccATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGTCTCCGCAATTTGCGGGCATCCATCAATATATGGGTGGGATGAT‑GGGAGGCGGC‑GGGAAAAAGATCCTCGGGGCAGACAGCCACACCC‑CCTACAGGGCAtt                                           >  2:272713/1‑150 (MQ=255)                                                                                                            gTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGCATTTTTTTTCCTCCGCCTTTTTCGGGCCTCCATCAATATATGCGTGGGATGGT‑GGGAGGCGGG‑GGCCAAAAGATCCTCCGGGCAGACAGCCACACCC‑CGTACGGGGGATTAagg                                       >  2:971200/1‑150 (MQ=255)                                                                                                                       aTCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGCCCCCGCCAATTTCGGGCCTTCCTCAATATTTGCGGGGGGTGAT‑GGGAGGGGGC‑GGGAAAAAGATCCCCCGGTGAGACAGCCACCCCCGGTAAGGGGGCTTTAGGGCAAATgg                                >  2:234002/1‑147 (MQ=255)                                                                                                                          ggcggcCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGAT‑GGCAGGCGGC‑GGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑GTTACGGTGCATTAGGGACAATGGCAGTCGGt                        <  2:181358/151‑1 (MQ=255)                                                                                                                            cggcCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCCTATTTCGGTCATCCCTCAATATATGCGTGAGATGAT‑GGCAGGGGGC‑GGCAAAATGATCCTCGGGGCAGACAGCCACACCC‑GTTACGGTGCATTAGGGACAATGGCAGTCgggg                       >  2:200431/1‑148 (MQ=255)                                                                                                                                 cAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGAT‑GGCAGGCGGC‑GGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑GTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGggcg                 <  2:276722/151‑1 (MQ=255)                                                                                                                                 cAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGAT‑GGCAGGCGGC‑GGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑GCTACGGTGCATTAGGGACAATGGCAGTCGGTGAGggc                  >  1:652035/1‑150 (MQ=255)                                                                                                                                        aTGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGAT‑GGCAGGCGGC‑GGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑GTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTg           <  2:274346/150‑1 (MQ=255)                                                                                                                                                 tttttGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGAT‑GGCAGGCGGC‑GGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑GTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTaa  >  2:765290/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                    AGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGAT‑GGCAGGCGGC‑GGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑GTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAA  >  NC_000913/803720‑804013

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: ‑