Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,575	0	A	C	55.4%	62.2 / 88.9	56	V555G (GTG→GGG)	ybaL	inner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (22/3);  new base C (0/31);  total (22/34)
Fisher's exact test for biased strand distribution p-value = 1.07e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-10

GTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTT  >  NC_000913/501436‑501724                                                                                                                                            |                                                                                                                                                      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)                                                                                                                                 cccccGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCaa             <  1:1166533/149‑1 (MQ=255)                                                                                                                                   cccGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAatat          >  2:1126605/1‑151 (MQ=255)                                                                                                                                     cGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATcc        >  1:1192090/1‑151 (MQ=255)                                                                                                                                           caccTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGAttt  >  2:1199924/1‑151 (MQ=255)                                                                                                                                            |                                                                                                                                                      GTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTT  >  NC_000913/501436‑501724

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