Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,575	0	A	C	43.4%	85.3 / 13.2	53	V555G (GTG→GGG)	ybaL	inner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (24/6);  new base C (0/23);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 7.75e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTT  >  NC_000913/501431‑501724                                                                                                                                                 |                                                                                                                                                      ttATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCAc                                                                                                                                                 >  1:744814/1‑151 (MQ=255)   aTCGTCGGTTCGTAGGCCTGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCAccc                                                                                                                                               <  1:309211/151‑1 (MQ=255)    tcgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCAccc                                                                                                                                               <  2:758756/150‑1 (MQ=255)     cgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCACCCACCGTTTACCCCGTCACCCCCTCACccgc                                                                                                                                             <  2:475432/151‑1 (MQ=255)     cgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCAccc                                                                                                                                               <  1:559786/149‑1 (MQ=255)     cgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACccg                                                                                                                                              <  1:305898/150‑1 (MQ=255)           ttCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTgg                                                                                                                                       >  1:238373/1‑151 (MQ=255)             cGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCg                                                                                                                                     >  2:976878/1‑151 (MQ=255)              gTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGt                                                                                                                                    >  1:56549/1‑151 (MQ=255)               tAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGtt                                                                                                                                   >  2:1095/1‑151 (MQ=255)                  gCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGCGGCGTTTcc                                                                                                                                >  2:320680/1‑151 (MQ=255)                      gATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagc                                                                                                                             >  1:864900/1‑150 (MQ=255)                       aaagggcttccccgccgcATCTGGCTTTTGCCTCCCGTTGCCTGATGCGCCGCGGGCGCGTCTTATCATGCCTACTTTTTTTCCATATTTTACACCCGGCACCCACCTTTTCCCCCGCCCCCCCCTCACCCGCCGGTGGCGTTTCca                                                                                                                               <  2:642925/135‑1 (MQ=255)                         agggCGTTCCCGCCGCATCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATCCCTACATTTTTTTCATTTTTTCCACCCGGCACCCCCCGTTTCCCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGtt                                                                                                                         <  2:380410/149‑1 (MQ=25)                              gTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGc                                                                                                                    >  1:166736/1‑151 (MQ=255)                                tCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt                                                                                                                  >  1:416162/1‑151 (MQ=255)                                   cgcagcAGCTGGCTTTTGGCTCTCGACCCCTGAGGCGACGCTGGCGCGTTTTATCATCCTTACATTTTTTTCATTTTTAACACCCGGCACCAACCGTTTCCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATgg                                                                                                                <  2:915135/150‑1 (MQ=255)                                     ccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATCCCTACATATTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGCCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTa                                                                                                              <  1:990176/150‑1 (MQ=38)                                             ggCTTTTGGCTCTCGATGCCTGATGCGCCGCGGGCGCGTCTTATCATCCCTACTTTTTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCg                                                                                                        <  2:168048/148‑1 (MQ=255)                                               ctttgggCTCTCGAGCCCTGATGCGCCGCGGGCGCTTTTTATCATGCCTACATTTTTTTCATATTTAACACCCGGCACCCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCg                                                                                                        <  2:370550/141‑1 (MQ=255)                                               cATTTGGCTCTCGATGCCTGATCCGACGCTGGCGCGTCTTACCATCCCTACATATTTTTCATTTTTTACACCCGGCACCCCCCGTTTCCCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTc                                                                                                   <  1:396871/151‑1 (MQ=255)                                                    ggCCCTCGATCCCTGATGCGCCGCTGGCGCTTCTTATCATCCCTACATTTTTTTCATTTTTTACACCCGGCACCCACCTTTTCCCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                                                                               <  1:673993/150‑1 (MQ=255)                                                        cccGATGCCTGATGCGACGCTGGCGCGTTTTATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCACCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                                           <  1:860347/148‑1 (MQ=38)                                                          cgttGCCTGAGGCGACGCTGGCGCGTTTTATCATGCCTACATATTTTTCATTTTTTACACCCGCCACCCCCCGTTTCCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                                           <  1:1001038/145‑1 (MQ=255)                                                           gATGCCTGAGGCGCCGCGGGCGCGTCTTATCATCCCTACATTTTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCg                                                                                            <  2:817152/146‑1 (MQ=255)                                                              gCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGCGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCAc                                                                                    >  2:4629/1‑151 (MQ=255)                                                               cttaaTGAGACCCTGCCGCGTCTTATCACGTTTCCTTATTTTCCATTTTTCCCCTCCGGAAACCCCGTTTTACCCCGCCCCCCCCCCTCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATc                                                                                      <  2:434650/144‑1 (MQ=255)                                                                   aTGCGCCGTGGGCGCGTCTTATCATCCCTACATATTTTTCATATTTTACCCCCGGCACCCCCCGTTTCCCCCGTCACCCCCCCGCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAcc                                                                                <  1:323240/150‑1 (MQ=255)                                                                      cgccGCTGGCGCGTCTTACCGTCCCAACATTTTTTCTATATTTCCCACCCGACAACCCCCTTTCCCCCCATCCCCCCCACACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAt                                                                            <  1:885341/148‑1 (MQ=255)                                                                       gACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt                                                                           >  2:215354/1‑151 (MQ=255)                                                                             ggCGCGTTTTATCATCCCTACATTTTTTTCATTTTTTCCATCCGGCACCCCCCGTTTCCCCCGCCACCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcac                                                                      <  1:594741/150‑1 (MQ=255)                                                                               cgcgTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcaccac                                                                   >  2:566314/1‑151 (MQ=255)                                                                                    cTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                                              >  1:888738/1‑151 (MQ=255)                                                                                            ccctccttttttttCATATTTTCCACCCGGCACCCCCCTTTTCCCCCTTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGAtata                                                      <  2:118298/142‑1 (MQ=255)                                                                                             ccTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                                     >  1:149150/1‑151 (MQ=255)                                                                                                aCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGc                                                   >  1:138060/1‑150 (MQ=255)                                                                                                   tatTTTTCATATTTTCCATCCGGCACCCACCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCa                                                 <  1:299488/149‑1 (MQ=255)                                                                                                      ttttCATATTTTCCACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTc                                             <  1:581085/150‑1 (MQ=255)                                                                                                      ttttCAAATTTCACACCCGGCACCCCCCTTTTCCCCCTTCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTc                                             <  1:935736/150‑1 (MQ=255)                                                                                                            tatTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTc                                       >  2:312343/1‑150 (MQ=255)                                                                                                              ttttACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt                                     >  1:639417/1‑150 (MQ=255)                                                                                                                    aTCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgc                              >  1:435023/1‑151 (MQ=255)                                                                                                                     cccGCCACCCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg                             <  2:864879/150‑1 (MQ=255)                                                                                                                      cggtcaccccatgTTTCCCCCGTCCCCCCCTCTCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcgc                            <  1:678441/139‑1 (MQ=255)                                                                                                                        ggCACCCCCCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg                          <  1:20064/151‑1 (MQ=255)                                                                                                                        ggCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg                          >  2:639982/1‑151 (MQ=255)                                                                                                                           aaCCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGc                        >  2:228157/1‑150 (MQ=255)                                                                                                                              ccccttttCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaata                    <  1:396832/146‑1 (MQ=255)                                                                                                                               accGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaata                    >  1:76533/1‑150 (MQ=255)                                                                                                                                 cGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtc                 >  2:724249/1‑151 (MQ=255)                                                                                                                                   ttccccccGCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtct                <  2:319785/146‑1 (MQ=255)                                                                                                                                        cccGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAatat          >  1:617865/1‑151 (MQ=255)                                                                                                                                           gTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCg       >  2:373713/1‑151 (MQ=255)                                                                                                                                             caccccCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGa     <  1:317587/151‑1 (MQ=255)                                                                                                                                                caccTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGAttt  >  1:792285/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                      TTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTT  >  NC_000913/501431‑501724

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