Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	298,960	0	A	C	42.9%	91.3 / 15.1	56	Y659D (TAT→GAT)	paoC	PaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base A (29/3);  new base C (0/24);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 6.72e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACAC  >  NC_000913/298810‑299107                                                                                                                                                        |                                                                                                                                                    cGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATa                                                                                                                                                     <  1:335181/151‑1 (MQ=255)     ataCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATAcccc                                                                                                                                                 >  1:882420/1‑151 (MQ=255)       acacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGc                                                                                                                                               >  2:848657/1‑151 (MQ=255)          ccGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCTCCTCATCCCCCGCCat                                                                                                                                            <  1:863746/151‑1 (MQ=255)                  gCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGAt                                                                                                                                    >  2:396009/1‑151 (MQ=255)                    gATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGAtt                                                                                                                                   >  2:733771/1‑150 (MQ=255)                     aTAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTa                                                                                                                                  >  2:154641/1‑150 (MQ=255)                        gCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACg                                                                                                                               >  1:333657/1‑150 (MQ=255)                          cgcgCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGa                                                                                                                              >  2:824517/1‑149 (MQ=255)                              cTCCCGCCGAACGGGCCCAGCTCCCCAACCCCTTTGGCCTTCTGGGGGGAGGATATGGGGTCGGTATCATCCAGAAAAAC‑CCCCCCCTTTTTTGGGAAACCCGCATGACCCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAg                                                                                                                         <  1:503825/150‑1 (MQ=255)                                  cgccgcACAGGCCCGGCTCCCCGACCCCTTTGGCCTTCATGGGGGAGGATATGGGGTCGGTATCACCCAGAAAAAT‑CCCCTCCTTTTTTGGGATACCCGCATGAACCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCa                                                                                                                     <  2:275109/150‑1 (MQ=255)                                    ccgcACAGGCCCAGCTCCCCGCCACCTTTGGCCTTCACCGGGGAGGATATGGGGTCGGTATCACCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCaaa                                                                                                                   <  2:474448/150‑1 (MQ=255)                                    ccgcACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAAc                                                                                                                  >  1:72064/1‑151 (MQ=255)                                    ccgcACAGCCCCAGCTCACCGACACCTTGGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGAAAAAT‑CCCCTCCTGTTTGGGGATATCCGCAGGAACCGCCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCaaa                                                                                                                   <  2:749548/150‑1 (MQ=255)                                       cacaGCCCCGCCACACCGCCCCCTTTGGCTTTCCGGGGGGAGAATAGGGGGTCGGTATCCGCCGGAAAAAT‑CCCCTCCTGTTTGGGAATATCCGCAAGACCCGGCCCCCCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGt                                                                                                               <  1:255609/151‑1 (MQ=255)                                        acaGGCCCAGCTCCCCGACCCTTTTGGCCTTCATGGGGGAGGATAGGGGGTCGGTATCATCCAGAAAAAT‑CCCCTCCTGTTTTGGGATACCCGCATGAACCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTc                                                                                                              <  2:962457/151‑1 (MQ=255)                                         caGACCCAGCTCACCGACACCTTTGGCCTTCAGCGGGGAGGATATGGGGTCGGCATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCaaa                                                                                                                   >  1:187809/1‑145 (MQ=255)                                                 gCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGc                                                                                                     >  1:121149/1‑151 (MQ=255)                                                   tCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGcc                                                                                                    >  1:424233/1‑150 (MQ=255)                                                       cGACACCTTTGGCCTTCATCGGGGAGGATATGAGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctc                                                                                               >  2:564154/1‑151 (MQ=255)                                                         acacCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcc                                                                                              >  1:720779/1‑150 (MQ=255)                                                          cacCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcctc                                                                                            >  2:92833/1‑151 (MQ=255)                                                            ctttggccctccaggggggAGGATGTGGGGTCGGTATCACCCAGGAAAAC‑CCCCCCTTGTTTTGGAATACCCACGTAACCCGCCACCCAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCtcctcc                                                                                           <  1:271831/135‑1 (MQ=255)                                                                 ggCCTTCATGGGGGAGGATAGGGGGTCGGTATCACCCGGAAAAAT‑CCCCTCCTGTTTTGGAATATCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGt                                                                                     <  2:42356/151‑1 (MQ=255)                                                                   ccTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGc                                                                                   >  1:593664/1‑151 (MQ=255)                                                                     ttcaggggggAGGATAGGGGGTCGTTACCACCAAAAAAAAC‑CCCCTCCTTTTTGGGGATACCCACATGACCCGGCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTg                                                                                    <  1:559488/142‑1 (MQ=255)                                                                       cacgggggAGAATAGGGGGTGGTTTTCA‑CCAAGAAAATCCCCCTCCTGTTTGGGGATACCCGAATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcgcg                                                                               <  1:51868/147‑1 (MQ=255)                                                                            gggAGGATGGGGGGTCGTTATCACCCAGAAAACT‑CCCCTCTTTTTTTGGTATACCCACAGGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccat                                                                          <  2:254961/151‑1 (MQ=255)                                                                                    atGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat                                                                    >  2:1070587/1‑149 (MQ=255)                                                                                    atGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat                                                                    >  2:926974/1‑149 (MQ=255)                                                                                      ggggTGGGTATCACCCAGAAAAAT‑CCCCTCCTTTTTTGGGATACCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAg                                                                  <  2:525142/149‑1 (MQ=255)                                                                                      ggggTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                                >  2:58732/1‑151 (MQ=255)                                                                                       gggTCGTTATCACCCAGAAAACC‑CCCCCCCTTTTTGGGAATTCCCGCATGACCCGCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCa                                                               <  1:993175/151‑1 (MQ=255)                                                                                           cGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgc                                                           >  1:382337/1‑151 (MQ=255)                                                                                            ggTATCATCCAGAAAACC‑CCCCCCCTGTTTGGGGATACCCGCATGAACGGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTg                                                            <  1:496945/149‑1 (MQ=255)                                                                                            ggTATCACCCAGAAAAAC‑CCCCTCCTGTTTGGGGATACCCGAATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgc                                                           <  2:141017/150‑1 (MQ=255)                                                                                                    cagagAAAAT‑CACCCCCTTTTGTGGGATACCCGCAGAACCCGGCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAc                                                  <  1:358316/147‑1 (MQ=255)                                                                                                     cAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAcc                                                 >  2:630925/1‑151 (MQ=255)                                                                                                      aGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCt                                                >  1:171470/1‑151 (MQ=255)                                                                                                         aaaaT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGc                                             >  1:892054/1‑151 (MQ=255)                                                                                                         aaaaT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGc                                             >  1:527552/1‑151 (MQ=255)                                                                                                           aaacccccccctGTTTGGGGATACCCGCATGACCCGGCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCt                                            <  2:333661/142‑1 (MQ=255)                                                                                                               cccctcctTTTTTGGGATACCCGCATAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgc                                        <  2:103255/149‑1 (MQ=255)                                                                                                                     ctGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtt                                  >  2:73956/1‑151 (MQ=255)                                                                                                                       gTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtttt                                >  1:611904/1‑151 (MQ=255)                                                                                                                        ttttGGGATACCCGCATAACCCGCCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt                                 <  2:862667/149‑1 (MQ=255)                                                                                                                        tttgGGAATACCCGCAGGAACCGCCCCCTCATCCCCCGCCATAGCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGGCCTGGCTGCGCGCAGtttt                                <  2:135926/150‑1 (MQ=255)                                                                                                                         ttggggATACCCGCATGACCCGCCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTc                               <  1:510863/147‑1 (MQ=255)                                                                                                                         ttggggATACCCGCATACCCCGCCCCCTAACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTc                               <  1:653713/147‑1 (MQ=255)                                                                                                                            gggATATCCGCATGACCCGGCCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGa                            <  2:629698/150‑1 (MQ=255)                                                                                                                            gggATAGCGGCATGACCCGCCACCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGa                            <  2:340205/150‑1 (MQ=255)                                                                                                                             ggATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGAtt                          >  2:202512/1‑151 (MQ=255)                                                                                                                                tacccGCAGGAACCGGCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCa                        <  2:766768/147‑1 (MQ=255)                                                                                                                                          accccacACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCt             <  2:882439/145‑1 (MQ=255)                                                                                                                                                 aCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcac      >  2:852282/1‑151 (MQ=255)                                                                                                                                                    tCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacac    >  1:82493/1‑150 (MQ=255)                                                                                                                                                       tACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacacac  >  1:212720/1‑149 (MQ=255)                                                                                                                                                        |                                                                                                                                                    CGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAT‑CACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACAC  >  NC_000913/298810‑299107

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