Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,621,118	0	A	C	44.4%	77.0 / 12.6	54	V7G (GTG→GGG)	eamA	cysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base A (23/7);  new base C (0/24);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 2.69e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTT  >  NC_000913/1620973‑1621267                                                                                                                                                  |                                                                                                                                                       tttCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACAcccc                                                                                                                                                    >  1:685623/1‑150 (MQ=255)   tCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATAAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATc                                                                                                                                                 >  1:363450/1‑151 (MQ=255)    cGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCt                                                                                                                                                <  2:741479/151‑1 (MQ=255)         gTGCGACAAAAAAGATAGCGGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTc                                                                                                                                            >  1:1014408/1‑150 (MQ=255)                   aaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTa                                                                                                                                <  2:218359/151‑1 (MQ=255)                       aTAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGc                                                                                                                            >  1:898735/1‑151 (MQ=255)                           ccGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACC‑TTCTACCGCctct                                                                                                                        >  1:542953/1‑151 (MQ=255)                               aaaaGCGCCCAGCATAAGGCCAAAACCGCCCGGCATCAGGGGTGCCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCTTACGACCCCTACAAGTAGCGCCAACCCCCCCTCTTTTCG‑GAGATTTATACCGCCTCTGag                                                                                                                     <  1:534272/150‑2 (MQ=255)                                  aGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTc                                                                                                                 >  2:30382/1‑151 (MQ=255)                                            aTAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCtttt                                                                                                          <  2:718561/148‑1 (MQ=255)                                            aTAAAGCGCAAACCGGCCAGCATCAGCGGTGCCATGTTATGAAGCCCCACTTTGATGCCCACAAAATTTAGCCCCCATACGACCCCTACCAGTAGCGCCAACCCCCCATCTTTCCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTa                                                                                                       <  1:146144/151‑1 (MQ=255)                                              aaaGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTa                                                                                                       >  2:421720/1‑149 (MQ=255)                                              aaaGCGCAAACCGGCCAGCATAAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAg                                                                                                     <  2:291567/151‑1 (MQ=255)                                                 gcgcAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCACCACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCa                                                                                                   >  1:593057/1‑150 (MQ=255)                                                  cgaaaaCCGGCCAGCATCAGCGGGGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGCCCCCTACCAGTAGCGCCACCCCCCCATCTTTGCGCGAAA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCa                                                                                                   <  1:57279/146‑1 (MQ=255)                                                     aaaCCGGCCAGCATCGGCGGTGGCATGTTATGAAGCCCCATTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCTCCAACCCCCCATCTTTTCGCGACT‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATc                                                                                                <  1:974178/149‑1 (MQ=255)                                                     aaaCAGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAt                                                                                                 >  1:861710/1‑148 (MQ=255)                                                             cAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTg                                                                                      >  2:252175/1‑151 (MQ=255)                                                             cAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACCACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTg                                                                                      <  1:779407/151‑1 (MQ=255)                                                              aGCATCAGGGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCACCCCCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTg                                                                                      <  2:257930/150‑1 (MQ=255)                                                              aGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTg                                                                                      >  1:645579/1‑150 (MQ=255)                                                                cATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGaa                                                                                    <  2:646100/150‑1 (MQ=255)                                                                 aTCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGaa                                                                                    <  2:364481/149‑1 (MQ=255)                                                                 aTCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAtt                                                                                  >  1:612397/1‑151 (MQ=255)                                                                     gCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTa                                                                                <  2:165295/149‑1 (MQ=255)                                                                     gCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACt                                                                              >  1:579977/1‑151 (MQ=255)                                                                       ggTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCAAACCCCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTAc                                                                               <  2:326252/148‑1 (MQ=255)                                                                          ggCATGTTATGAAGCCCCACTTTTATGCCCAAAAAATTTAGCCCCCATCCGCCCCCTCCCATTACCGCCATCCCCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCt                                                                         <  1:1021150/151‑1 (MQ=255)                                                                           acATTTTAGGAAGCCCCTATGTGACAACCAAAAAATTTACCCCCCAGACGCCCATTACTAGTAGCGCAACCCCCCCGTCTTTTCCGGACG‑TTCTCCCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCtt                                                                        <  2:745476/150‑1 (MQ=255)                                                                                ttATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCCCTCCCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTa                                                                      <  2:651320/148‑1 (MQ=255)                                                                                  aTGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAg                                                                  <  2:155764/150‑1 (MQ=255)                                                                                   tGAACCCCCCTTTTGATAACCACAAAATTTAGCCCCCATACGACCTCTACCAGTAGCCCCAACCCCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACa                                                                   <  2:866816/148‑1 (MQ=255)                                                                                     aaGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTg                                                               >  1:58387/1‑150 (MQ=255)                                                                                        ccccACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTc                                                           >  1:183512/1‑151 (MQ=255)                                                                                        ccccACTTTGATGACCAAAAAATTTAGCCCCCATACGACCACTCCAAGTAGCCCCAACCCCCCTTCTTTCCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTg                                                               <  1:591682/147‑1 (MQ=255)                                                                                          ccACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGt                                                         >  2:111704/1‑151 (MQ=255)                                                                                              tttGATGCCCACAAAATTTAGCCCCCATACGCCCACTCCCATTAGCGCCAACCCCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTa                                                       <  1:508870/149‑1 (MQ=255)                                                                                               ttGATGACCAAAAAATTTAGCCCCAATACGACCACTCCAAGTCGCGCCAACCCCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGt                                                     <  2:380409/150‑1 (MQ=255)                                                                                                tgatgaCCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCg                                                   >  1:201862/1‑151 (MQ=255)                                                                                                tgatgaCCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCCTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCg                                                   >  1:772290/1‑151 (MQ=255)                                                                                                      ccACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGtt                                             <  2:120075/151‑1 (MQ=255)                                                                                                              tttAGCCCCCATCCGCCCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACgg                                      <  2:538058/150‑1 (MQ=255)                                                                                                                tAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCAt                                   >  1:195728/1‑151 (MQ=255)                                                                                                                tAGCCCCCAGGCCACCACGCCCAGTGGCGCCATCCCCCCATCCTTGCGCCACA‑TTTTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCa                                    <  1:919274/150‑1 (MQ=255)                                                                                                                  gCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCAtt                                  >  1:1005868/1‑150 (MQ=255)                                                                                                                    ccccATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTGCCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAAt                               <  1:779185/151‑1 (MQ=255)                                                                                                                        aTACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgcgc                           >  2:470973/1‑151 (MQ=255)                                                                                                                           ccaCCCCTCCCAGTAGCGCCACCCCCCCATCTTTTCGCGACA‑TTCAACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCaa                         <  2:1014783/148‑1 (MQ=255)                                                                                                                              ccACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAAt                       <  2:1024853/149‑1 (MQ=255)                                                                                                                                aCTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGt                   >  2:721699/1‑151 (MQ=255)                                                                                                                                    ccATTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCg                 <  2:295336/149‑1 (MQ=255)                                                                                                                                      aGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCtat             >  1:487517/1‑151 (MQ=255)                                                                                                                                           cgacAACACCCCATCTTTTCGGGACT‑TTCTACCGCCTCTGAATTTCATCTGTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAAGGCGCAAATAAGTCGCTATAc           <  1:575945/145‑1 (MQ=255)                                                                                                                                               aaCCCCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAt    <  2:418347/151‑1 (MQ=255)                                                                                                                                                 cacCCCATCTTCTGGCGACA‑TTCTACCCTCTCTGTATTTCCTCTTTTGTAAGCAATCCACTTAGCTGCATTTACTTTTCTGTAACAGTTGATTTGCTCGTCGCCGGTTACGACCGCATTAATGCGCAAATAAGTCGCTATACTTCGGAttt  >  1:157307/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                       TTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACA‑TTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTT  >  NC_000913/1620973‑1621267

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