Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,147,733	0	T	G	44.4%	90.0 / 13.5	54	G20G (GGT→GGG)	yegD	Hsp70 chaperone family protein
Reads supporting (aligned to +/- strand):  ref base T (3/27);  new base G (24/0);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 4.17e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTAC  >  NC_000913/2147586‑2147875                                                                                                                                                     |                                                                                                                                               gCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTaaa                                                                                                                                             <  1:342574/151‑1 (MQ=255)  cTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TTCGTGGCGGGaaa                                                                                                                                             >  2:402163/1‑150 (MQ=255)  cTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTaaa                                                                                                                                             <  2:491/150‑1 (MQ=255)   tGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGGGGCGGGAAAcc                                                                                                                                           >  1:786000/1‑151 (MQ=255)    gcgcTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGGaaa                                                                                                                                             >  2:863397/1‑148 (MQ=255)      gctttgctttCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCa                                                                                                                                        <  2:689581/151‑1 (MQ=255)       ctttgctttCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGGTA‑TGCGTGGCGGGAAAACGc                                                                                                                                         >  1:190688/1‑149 (MQ=255)         ttgctttCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCAt                                                                                                                                       <  2:728681/149‑1 (MQ=255)          tgctttCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCAttt                                                                                                                                     <  2:679901/150‑1 (MQ=255)             tttCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGGAGGGAAAACGCCTTTGCt                                                                                                                                  >  1:173151/1‑150 (MQ=255)                 atatTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTaaaaa                                                                                                                             <  2:594782/151‑1 (MQ=255)                   atTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGGAGGGAAAACGCATTTTGTaaaaa                                                                                                                             >  2:734483/1‑149 (MQ=255)                   atTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATg                                                                                                                           <  2:870092/151‑1 (MQ=255)                   atTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATg                                                                                                                           <  2:773299/151‑1 (MQ=255)                       cGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGaaa                                                                                                                       <  2:379845/151‑1 (MQ=255)                          ttGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TTCGTGGCGGGAAACCGCATTTGTTTAAAAAGGaaaa                                                                                                                      >  1:437453/1‑149 (MQ=255)                           tGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACg                                                                                                                    <  2:643378/150‑1 (MQ=255)                           tGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGGGACGGGAAAACGCCTTTTCTTAAAATGGaaaa                                                                                                                      >  2:593015/1‑148 (MQ=255)                              cgcgACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGGGACGGGAAAACGCCTTTTCTAAAAAAGGAAAAAGaaa                                                                                                                 >  2:700762/1‑148 (MQ=255)                               gcgACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCC‑TGCGTGACGGGAAAACCCACTTTCTAAAAAATGaaa                                                                                                                       >  1:542922/1‑143 (MQ=255)                                cgACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TTCGTGTAGGGAAACCCAATTTGCCTAAAATgggtaaacacagc                                                                                                               >  1:598821/1‑139 (MQ=255)                                  aCGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACg                                                                                                            <  1:779698/151‑1 (MQ=255)                                   cGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTTT‑TGCGTGTGGGGTAACCGCCCTTTCTTAAAAATGAAAAAGACa                                                                                                                 >  1:239162/1‑145 (MQ=255)                                             tCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCtt                                                                                                  <  1:371606/150‑1 (MQ=255)                                                    aGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCtt                                                                                          <  2:477755/151‑1 (MQ=255)                                                         cGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTgcg                                                                                      <  1:805699/150‑1 (MQ=255)                                                                tcttcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTGGGTG‑CGGGAAACACCATTTTGTTAAAATGGAAAACGACAGCACCCCGCTGCCCTCCCTACTTTTCGGGCCaa                                                                                >  1:395197/1‑149 (MQ=255)                                                                 cttcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCT‑TTGGTGGCGGGTAACCGCATTTTCTAAAAAAAGAAAAAGACAGCACGCTGCTGGCCTCTATTCTTTGTGCGccca                                                                                >  1:125838/1‑146 (MQ=255)                                                                  ttcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑AGCGTGACGGGAAAACGCCCTTGCTAAAAAAGGAAAAAGACAGCACGCTGCTGCCTTCAATTCTTTTCGCGCCCAcgcg                                                                            >  2:677487/1‑151 (MQ=255)                                                                    cACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTg                                                                          <  2:636727/151‑1 (MQ=255)                                                                    cACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGGGACGGGAAACCGCATTTTCTAAAAAAGGAAAACGACAACACGCTGCTGCCCTCAATTCTTTTCGCGccca                                                                                >  2:505070/1‑143 (MQ=255)                                                                     aCCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGGCGGTAAACCGCATTTTCTAAAAAAGGAAAAAGACAGCACACTGCTCCCCTCCATGCTTTGCGCGCCAACGCGTg                                                                          >  1:609942/1‑150 (MQ=255)                                                                              gggATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAg                                                                <  2:86414/151‑1 (MQ=255)                                                                                gATTATTTCGTGTTTATTGGTTTTGATTACCGTACAGCAAGCTGTTCAGCGGCGGTCA‑TGCGTGACGGTAAACCGCAGTTGCTAAAAATGGAACACTACAGCACGCTGCTGCCTTCGATGCTTTGCTCGCCGACGCGTGAATCGGTAAGCg                                                              <  1:585271/151‑1 (MQ=255)                                                                                  ttattTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGGGACGGGAAAACGCATTTTCTTAAAAAGGGAAAAGACAACACGCTGCTGCCTTCCATGCTTTTCGCGCCCACGCGTGAAGCGGGaa                                                                 >  2:18979/1‑146 (MQ=255)                                                                                   tattTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGaa                                                            <  1:169507/150‑1 (MQ=255)                                                                                    attTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATg                                                          <  2:231037/151‑1 (MQ=255)                                                                                      ttCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCC‑TGCGTGGAGGGAAACAGCCTTTTCTAAAAATAGAAAAAGACCGCACGCCGCTGCCTTCAATGCTTTTCGCGCCAACGCGAGAAGCGGGAAGCggaag                                                          >  1:438449/1‑145 (MQ=255)                                                                                      ttCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGGAAAACGCCTTTGCTAAAAAAGGAAAAAGAAAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCCACGCGTGAAGCGGTAAGCGaa                                                            >  2:324631/1‑147 (MQ=255)                                                                                       tCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGc                                                        <  2:385527/150‑1 (MQ=255)                                                                                          tgtTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCG‑TGCGTGACGGGAAAACGCATTTTCTAAAAATGGAAAAAGACAGCACGCCGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGAGGGAAGCGAATGGCTg                                                      >  2:375503/1‑149 (MQ=255)                                                                                           gtTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTAc                                                   <  2:819441/151‑1 (MQ=255)                                                                                              tATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGGAAACCGCATTTGCTTAAAATGGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCCACGCGTGAAGCGGTAAGCGAAATGCTGTACCGc                                                >  2:446987/1‑151 (MQ=255)                                                                                               cttGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGGAAAACGCAGTTGCTAAAAAAGGGAAAAGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGCCAACGCGTGAAACGGTGAACGAAAGGGCGTACCGcc                                               >  1:111346/2‑151 (MQ=255)                                                                                               aTTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGc                                                <  1:314012/150‑1 (MQ=255)                                                                                                        gATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCatgat                                       <  1:870459/150‑1 (MQ=255)                                                                                                        gATTACGGTACAGCAAACTGTTCAGTGGCGGCCA‑TGGGCGACGGGGAACCCCCTTTTCTGAAAAAGGAAAAAGACAACACGCCGCTGCCCTTAAAGCTTTGTGCGCCCACCCGTGAAGCGGGAAGCGAAAGGGTGGTCCGCcaccctgtt                                       >  1:727038/1‑143 (MQ=255)                                                                                                              ggTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGGCGGGAAAACGCATTTGCTAAAAAAGGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGCCCACGCGGGAAGCGGGAAGCGGAAGGCTGTACCGCCCTCATGATGTTCCgg                                >  2:22752/1‑151 (MQ=255)                                                                                                                   aGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAgac                           <  2:363190/151‑1 (MQ=255)                                                                                                                                gTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGc              <  2:786211/151‑1 (MQ=255)                                                                                                                                     ggTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGgcg         <  2:344942/151‑1 (MQ=255)                                                                                                                                      gTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGgcgc        >  2:207808/1‑151 (MQ=255)                                                                                                                                         a‑TGCGTGACGGTACACCGCCTTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGAGCCAACGCGTGAAGCGGTAAGCGAATGGCTGAACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGCCGCTGt     >  1:720660/1‑151 (MQ=255)                                                                                                                                             cGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTAc  <  2:703499/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                               GCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTAC  >  NC_000913/2147586‑2147875

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