Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,564	0	A	C	53.6%	66.0 / 92.7	56	*559E (TAA→GAA)	ybaL	inner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (22/4);  new base C (0/30);  total (22/34)
Fisher's exact test for biased strand distribution p-value = 6.98e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.71e-09

GGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCA  >  NC_000913/501419‑501712                                                                                                                                                  |                                                                                                                                                     gggATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCg                                                                                                                                                 <  2:491074/151‑1 (MQ=255)    aTCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCg                                                                                                                                                 >  2:601694/1‑148 (MQ=255)     tCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcac                                                                                                                                             >  1:651242/1‑151 (MQ=255)                  gtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACccgcc                                                                                                                                >  1:841407/1‑151 (MQ=255)                     ggTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCgg                                                                                                                              >  2:899191/1‑150 (MQ=255)                     ggTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGt                                                                                                                             >  2:19687/1‑151 (MQ=255)                          gTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGt                                                                                                                        >  2:498224/1‑151 (MQ=255)                          gTAGGCCAGATAAGGCGTCCACGCCGCATCGGGCATTTGGCTCTCGATGCCTGATGCGACGCGGGCGCGTTTTATCATCCCTCCATTTTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGTCACCCCCCCCCCCGCCGGTGGc                                                                                                                          <  1:248419/149‑1 (MQ=255)                            aGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGttt                                                                                                                      >  2:106182/1‑151 (MQ=255)                             ggCCAGAAAGGGCTTTCACGCCGCATCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCGGGCGCGTCTTATCATCCCTACATATTTTTCATTTTTTACACCCGGCACCCACCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGtt                                                                                                                       <  1:550544/149‑1 (MQ=255)                             ggCCAAATAGGGCTTCCCCGCCGCATCTGGCTTTTGGCCCTCGATGCCTGATGCGCCGCTGGCGCGTCTTATCATCCCTACATATTTTTCATATTTTACATCCGGCACCCCCCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTgg                                                                                                                           <  1:345462/145‑1 (MQ=255)                                    tAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGt                                                                                                              >  2:292316/1‑151 (MQ=255)                                     agggCGTCCCCCCCGCATCTGGCTTTTGGCTCTCGATGCCTGAGGCGACGCTGGCGCTTTTTATCATCCCTACATTTTTTTCATATTTTACATCCGGCACCCACCGTTTACCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTCca                                                                                                                   <  1:107208/143‑1 (MQ=255)                                      gggcttacccgccgcATGGGTTATTTGGCTCTCGATGCCGGAGGCGCCGCGGGCGCGTTTTATCATGCCTACATTTTTTTCATTTTTTACACCGGCCAACCACCTTTTCCCCCTCCACCCCCCCACCCGCCGGTGGCGTTTCcagcag                                                                                                               <  1:229299/139‑1 (MQ=255)                                       ggCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTcc                                                                                                           >  2:898583/1‑151 (MQ=255)                                            tCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt                                                                                                      >  2:569860/1‑151 (MQ=255)                                               cgccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                                   >  1:125557/1‑151 (MQ=255)                                                      tCTGGCTTTTGGCTCCCGATGCCTGAGGCGACGCTGGCGCGTTTTACCAGGCCTAAATTTTTTTCATATTTTACATCCGGCAACCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACggg                                                                                              <  1:549384/149‑1 (MQ=255)                                                         ggtttttGGCTCTGGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCACCCACCTTTTCCCCCGTCCCCACCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAt                                                                                          <  2:580424/146‑1 (MQ=38)                                                         ggCTTTTGGCTCTCGAGCCCTGAGGCGCCGCTGGCGCTTTTTACCATGCCTCCTTATTTTCCATTTTTTCCACCCGGAACCCACCGTTTCCCCCCTCACCACCTCTCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACggg                                                                                              <  1:1207439/146‑1 (MQ=255)                                                          gCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAttt                                                                                        >  2:685833/1‑151 (MQ=255)                                                           cATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTc                                                                                       >  1:1057121/1‑151 (MQ=255)                                                             tttGCCTCTCGATGCCTGATGCGCCGCGGGCGCGTCTTTTCATCCCTACATTTTTTTCATATTTTACACCCGGCACCCACCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCa                                                                                      <  1:138689/150‑1 (MQ=255)                                                              ttGGCTCTCGATGCTTGATGCGACGCGGGCGCTTTTTATCATGCCTACATTTTTTTCATTTTTTCCACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCAc                                                                                     <  1:924002/150‑1 (MQ=255)                                                                  ctctCGATCCCTGATGCGCCGCTGGCGCTTCTTATCATGCCTAAATTTTTTTCATATTTTACACCCGACAACCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTc                                                                                 <  1:806960/150‑1 (MQ=255)                                                                        aTGCCTGTTGCGGCGCTGGGGCTTTTTATTAGCCCTACATTTTTTTCATTTTTTACATCCGAAACCCACGTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATc                                                                          <  1:301079/151‑1 (MQ=255)                                                                           ccTGATGCGCCGCTGGCGCGTTTTATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCACCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                                       <  1:1082041/151‑1 (MQ=255)                                                                           ccTGATGCGACGCTGGCGCGTCTTATCATTCCTACATTTTTTTCTTTTTTTCCACCCGGCAACCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCa                                                                         <  1:82267/149‑1 (MQ=255)                                                                              gATGCGACGCGGGCGCGTTTTATCATGCCTACATTTTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCAc                                                                        <  1:392683/147‑1 (MQ=255)                                                                                 gcgcCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGCAACCCCCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAcc                                                                    <  2:754232/148‑1 (MQ=255)                                                                                   gACGCTGGCGCGTCTTATCATGCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAt                                                                <  1:430260/150‑1 (MQ=255)                                                                                      gcgggCGCGTCTTATCATCCCTACATTTTTTTCATTTTTTACACCCGGCAACCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCg                                                             <  2:627918/147‑1 (MQ=255)                                                                                         ggCGCTTTTTATCATGCCTATATTTTTTTCATTTTTTACACCCGCCAACCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTg                                                                  <  1:452636/142‑1 (MQ=255)                                                                                         ggCGCTTTTTATCATCCCTACATTTTTTTCATTTTTTACACCCGCCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAt                                                                <  2:752769/144‑1 (MQ=255)                                                                                           ccgattttACCATGCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCGTTTCCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacca                                                        <  1:362598/144‑1 (MQ=255)                                                                                                  tATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacca                                                        >  2:362601/1‑143 (MQ=255)                                                                                                      aTGCCTAATTATTTTTCATATTTTACAGCCGGAACCCCCCTTTTCCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGa                                              <  2:341793/149‑1 (MQ=255)                                                                                                       tGCCTACATTTTTTTTATATTTTCCATCCGGCACCCCCCGTTTCCCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGat                                             <  2:835127/149‑1 (MQ=255)                                                                                                         ccTACATTTTTTTCATATTTTACACCCGACAACCCCCGTTTCCCCCGCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTg                                               <  1:562607/145‑1 (MQ=255)                                                                                                           tACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGc                                       >  2:905850/1‑151 (MQ=255)                                                                                                           tACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGc                                       >  1:307709/1‑151 (MQ=255)                                                                                                           tACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGc                                       >  1:887602/1‑151 (MQ=255)                                                                                                                 tttttCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACt                                   <  1:414862/149‑1 (MQ=255)                                                                                                                           tttACGCCCGGCACCCCCGGTTACCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa                        <  2:547877/150‑1 (MQ=255)                                                                                                                           tttACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa                        >  2:314829/1‑150 (MQ=255)                                                                                                                              acacCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATg                     <  1:1217610/150‑1 (MQ=255)                                                                                                                               cATCCGCCACCCACCTTTTACCCCGTCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATggg                   <  2:768904/151‑1 (MQ=255)                                                                                                                                aTCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATgg                    <  2:1236776/149‑1 (MQ=255)                                                                                                                                aTCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgc                  >  2:1059669/1‑151 (MQ=255)                                                                                                                                  ccGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcgc                >  2:791305/1‑151 (MQ=255)                                                                                                                                    ggCACCCCCCGTTTCCCCCGCCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg              <  2:1102122/151‑1 (MQ=255)                                                                                                                                       acccacgttttCCCCCTTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCggg             <  1:401144/141‑1 (MQ=255)                                                                                                                                         ccaccGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaat         >  1:638730/1‑151 (MQ=255)                                                                                                                                          caccGTTTCCCCCGCCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaata        <  1:608374/151‑1 (MQ=255)                                                                                                                                               tttACCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCG‑CCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtct    <  1:969358/149‑1 (MQ=255)                                                                                                                                                 tACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCa  >  2:416830/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCA  >  NC_000913/501419‑501712

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