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

CGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATT  >  NC_000913/1620976‑1621266                                                                                                                                               |                                                                                                                                                     cGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCt                                                                                                                                              <  1:512003/151‑1 (MQ=255)  ggTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCtt                                                                                                                                             >  1:983250/1‑151 (MQ=255)   gTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCttt                                                                                                                                            >  2:29102/1‑151 (MQ=255)      gTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTc                                                                                                                                          >  1:236592/1‑150 (MQ=255)       tGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTcgc                                                                                                                                        >  2:717489/1‑151 (MQ=255)         cGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTcgcg                                                                                                                                       <  2:279669/150‑1 (MQ=255)             aaaaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACAtt                                                                                                                                  <  1:675767/151‑1 (MQ=255)               aaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCt                                                                                                                                <  2:613738/151‑1 (MQ=255)                       gCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCctc                                                                                                                        <  1:316001/151‑1 (MQ=255)                        ccGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCctc                                                                                                                        >  2:459876/1‑150 (MQ=255)                                    ccAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGCCCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCt                                                                                                            <  2:993078/150‑1 (MQ=255)                                    ccAGCATAAAGCGAAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGCCCAAAAAATTTAGCCCCCATACGCCCACTCCCAGTAGCCCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATc                                                                                                             <  1:200108/149‑1 (MQ=255)                                     cAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCt                                                                                                            >  1:269821/1‑149 (MQ=255)                                            aaGCGCAACCCGGCCACCATCAGCGGTGACATTTTATGAACCCCCACTTTGACGGCCACAAATTTTAACCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCTCGGCATTCTACCGCCTCTGAATTTCATCTTTTGTaa                                                                                                     <  1:871671/149‑1 (MQ=255)                                            aaGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGc                                                                                                   <  2:256723/151‑1 (MQ=255)                                                 cAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATc                                                                                               >  1:1123559/1‑150 (MQ=255)                                                  aacccGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATc                                                                                               <  2:856263/146‑1 (MQ=255)                                                  aaaCCGGCCAGCATCAGCGGTGGCATGTTATGAACCCCCCCTTTGATGCCCACAAAATTTACCCCCCATACGACCCCTCCCAGTAGCGCCACCCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCaa                                                                                             <  1:1048613/151‑1 (MQ=255)                                                   accgggCCAGCATCGGCGGTGCCATTTTAGGAAGCCCCCCTTGAATACCCAAAAAATTTGCCCCCCACACGCCCACTGCCGGAAGCCCCACTCCCCCTTCTTTCCGCGAAATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCa                                                                                              <  1:624371/145‑1 (MQ=255)                                                       ggCCGGCATCAGCGGTGCCATTTTAAGAAGCCCCACTTTAATCCCCACAAATTTTAGCCCCCATACGCCCTCTACCATAAGCGCCAACCCCCCATCTTTTCCCGACATTCTTCCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTa                                                                                         <  2:1074791/150‑1 (MQ=255)                                                       ggCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACtt                                                                                          >  2:266372/1‑149 (MQ=255)                                                         ccaccaTCAGGGGTGGATTGTTATGAACCCCCACTTTGAGGCCCAAAAAATTTAGCCCCCATCCACCCACTACCAGTAGCGCCAACCCCCCATCTTTCCGCGACATTCCACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTa                                                                                         <  2:869981/148‑1 (MQ=255)                                                                cAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAttt                                                                                >  2:719941/1‑150 (MQ=255)                                                                 aGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTAc                                                                              >  1:421554/1‑151 (MQ=255)                                                                  gCGGTGGCATTTTATGACGCCCCCCTTTGATGCCCACAAAATTTAGCCCCCATACGCCCACTCCCAGTAGCGCCAACCCCCCTTCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACt                                                                             <  1:480259/151‑1 (MQ=255)                                                                     gTGGTATTTTATGAAGCCCCCCTTTGATGACCACAAAATTTAGCCCCCATCCGACCACTACCAGTAGCGCCAACCCCCCATCTTTCCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtttt                                                                          <  1:367008/151‑1 (MQ=255)                                                                       ggCATTTTATGAACCCCCCCTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCAAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtttt                                                                          <  2:899238/149‑1 (MQ=255)                                                                             ttATGAACCCCCACTTTGATGCCCAAAAAATTTACCCCCCATACGCCCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTaa                                                                    <  2:956392/149‑1 (MQ=255)                                                                                 gAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTg                                                              >  1:1000512/1‑151 (MQ=255)                                                                                  aacccccACTTTGATAACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGa                                                             <  2:750717/148‑1 (MQ=255)                                                                                      cccACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCCCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCg                                                         <  2:517138/151‑1 (MQ=255)                                                                                       ccACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGCAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTc                                                          >  1:22733/1‑149 (MQ=255)                                                                                             tgatgaCCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCg                                                  >  1:131072/1‑151 (MQ=255)                                                                                                   ccACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGtt                                            >  1:907384/1‑151 (MQ=255)                                                                                                   ccACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGtt                                            >  1:924906/1‑151 (MQ=255)                                                                                                       aaaaTTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacg                                         >  2:216474/1‑150 (MQ=255)                                                                                                         aaTTTAGCCCCCATACGACCCCTCCCAGTAGCGCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgacg                                      <  1:305779/151‑1 (MQ=255)                                                                                                         aaTTTACCCCCCATACGACCCCTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTTTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgacg                                      <  1:614991/151‑1 (MQ=255)                                                                                                               gCCCCCATACGACCATTACCAGTACCGCCAACACCCCATCTTTGCGCGCCATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGc                                    <  2:434959/147‑1 (MQ=255)                                                                                                                cccccATACGCCCCCTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTaa                               <  2:435806/151‑1 (MQ=255)                                                                                                                 ccccATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGCATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAAt                              >  1:55445/1‑151 (MQ=255)                                                                                                                       aCGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCa                         >  1:1025041/1‑150 (MQ=255)                                                                                                                        cGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCaa                        >  1:398486/1‑150 (MQ=255)                                                                                                                           ccACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAAt                      <  2:33072/149‑1 (MQ=255)                                                                                                                            cACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATa                     >  2:72325/1‑149 (MQ=255)                                                                                                                                 ccAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGc               <  1:449003/150‑1 (MQ=255)                                                                                                                                   aGTGGCGCCACCCCCCCATCTTTCCCCAACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCtat            <  2:494655/151‑1 (MQ=255)                                                                                                                                     tggcgcCACCCCCCCATCTTTCCGCGACATTCTGCCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCtata           <  1:209582/148‑1 (MQ=255)                                                                                                                                     tAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCACAAATAAGTCGCtata           <  2:263637/150‑1 (MQ=255)                                                                                                                                        cgcAAACCCCCAACCTTTTCCCGACATTCAACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTc       <  1:725126/151‑1 (MQ=255)                                                                                                                                          ccAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTc       <  2:820615/149‑1 (MQ=255)                                                                                                                                            aaCACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAt   >  1:309504/1‑151 (MQ=255)                                                                                                                                             acacCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAtt  >  1:247687/1‑151 (MQ=255)                                                                                                                                               |                                                                                                                                                     CGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATT  >  NC_000913/1620976‑1621266

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