Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,913	0	G	C	44.2%	89.4 / 11.6	52	G44G (GGC→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (22/7);  new base C (0/23);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 8.45e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAA  >  NC_000913/993768‑994061                                                                                                                                                   |                                                                                                                                                      gTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCccac                                                                                                                                                  <  2:307028/151‑1 (MQ=255)   aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCccacc                                                                                                                                                 <  1:58646/150‑1 (MQ=255)    gCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCcaccac                                                                                                                                               >  2:147053/1‑151 (MQ=255)        aCGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGcc                                                                                                                                           <  1:233545/151‑1 (MQ=255)         cgcgCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCca                                                                                                                                          <  2:423101/151‑1 (MQ=255)         cgcgCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCca                                                                                                                                          <  1:601111/151‑1 (MQ=255)             cATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACaa                                                                                                                                       >  1:308578/1‑150 (MQ=255)                     aaCATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTg                                                                                                                              >  2:211219/1‑151 (MQ=255)                       catcatTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGcc                                                                                                                            >  2:422654/1‑151 (MQ=255)                          cttccgcttttCTTCCTAAATTC‑CAGCCGGGGGTGGGGTCCCCCCTACCCCCCCCCCTGGGTTTGGCTTTCCCACCCCCCCCAGCAGGCCCAGCAGGGTTTTTTTCCCCCCCCCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGcc                                                                                                                            <  2:149927/138‑1 (MQ=255)                          catTCGCGTTTCTTCCTGATTTT‑CAGCCGGGGGTGGGGTCCCCGCTAACACATCGCCTGGGTTTGGCTTTCCCAGCCCGGCCAGCAGGCGCAGCAGGGTTTTTTTGCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGa                                                                                                                         <  1:101579/151‑1 (MQ=255)                              cgcgTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAatat                                                                                                                     >  2:423270/1‑151 (MQ=255)                                  tgtCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTa                                                                                                                  >  1:288337/1‑150 (MQ=255)                                        ccggatttt‑CACCCGGCGGGGGGGTCCCCGCTACCCCCCCGCCGGCGGTTGGCTTTCCCACCCCTGCCAGCGGCCGCAGCAGGTTATTTTTCCCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                   <  2:258650/137‑1 (MQ=255)                                         cctaaTTTCCAGCCGGGGGTGGGGGCCCCGCTAACACACCGCCGGGGGTGGGCTTTCCCAGCCCTGCCAGCGGGCGCAGCGGGTTACTTTTCCCCCCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTaa                                                                                                                 <  1:761260/142‑1 (MQ=255)                                         cTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑Gt                                                                                                         >  1:891665/1‑151 (MQ=255)                                            aaTTC‑CGGCCAGGGGTGGGGTCCCCGCTACCACCCCGCCGGCGTTTGGCTTTCCCAGCCCCGCCAGCGGGCGCAGCAGGTTCTTTTTCCCCCCACCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑Gtt                                                                                                        <  2:320795/149‑1 (MQ=255)                                             aTTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGt                                                                                                     >  1:441396/1‑151 (MQ=255)                                             aTTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGG‑AATTTTGCCACCCCCCCCGCGGCCCAACACCGCCCCAAAATGACCTTCCGGGATAATTAAATCCAGGTTGGt                                                                                                     >  1:493765/1‑151 (MQ=255)                                                tccagccggcGGTGTGGTCCCCGCTAACACACCGCCTGCGGTTGGCTTTCCCAGCCCTCCCAGCAGGCGCAGCAGGGTACTTTTGCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGtt                                                                                                    <  1:532777/148‑1 (MQ=255)                                                     ccaggGGTGTGGGCCCCGCTAACACACCGCCTGCGGTTGGCTTTCCCAGCCCTGCCAGCAGGCGCAGCGGGTTACTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGAc                                                                                              <  1:617959/151‑1 (MQ=255)                                                     ccagcGGTGTGGTCCCCGCTAACCCATCCCCTGGGGTTGGCTTTCCCAGCCCTGCCAGCGGGCGCAGCGGGTTCCTTTGCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGa                                                                                               <  1:548361/150‑1 (MQ=255)                                                         cGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGat                                                                                           >  2:728068/1‑150 (MQ=255)                                                          ggTGTGGTGCCCGCAAACAAACCGCCGGGGGTTGGCGTTTCCAGCCCCGCGAGCGGGCGGAGCGGGTTATTTTTGCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACg                                                                                             <  1:104803/147‑1 (MQ=255)                                                          ggTGGGGTCCCCCCAAACCCACCGCCGGGGTTGGGCTTTCCCAGCCCCCCCAGCGGCCGCAGAGGGGTCTTTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGata                                                                                          <  2:138663/150‑1 (MQ=255)                                                             ggggtcCCCGCAACCCCACCGCCTGCGTTTGGCGTTCCCAGCCCTGCCAGCAGGCGCAGCAGGGTTCTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATAttt                                                                                       <  2:700005/148‑1 (MQ=255)                                                                   cccGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGt                                                                                >  1:804339/1‑151 (MQ=255)                                                                   cccGCTAACACACCGCCTGCGGTTGGCTTTCCCACCCCGGCCAGCAGGCGCAGCGGGTTTTTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCcg                                                                                   <  2:723146/148‑1 (MQ=255)                                                                     cGCTAACACATCGCCTGGGGTTGGTTTTTCCACCCCTGCCAGCGGGCGCAGCGGGGTACTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTaa                                                                              <  1:757709/151‑1 (MQ=255)                                                                         aaCACACCCCGTGGGTTTGGCGTTCCCAGCCCTGCAAGCGGGCGCAGCGGGTTACTTTTGCCCCCACCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGt                                                                           <  2:172778/150‑1 (MQ=255)                                                                             caTCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCt                                                                      >  1:149388/1‑151 (MQ=255)                                                                                       ggtgggtgTTCCCGGTCCTCCCAGCAGCCGCGGAAGGTTATTTTGCCCACCCCCGCGGCGCCCCACCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCa                                                               <  1:386860/148‑1 (MQ=255)                                                                                         ttGGCGTTTCCAGGCCCCCCAGAGGCCGCGGCGGGGTTTTTTTGCCCCCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGa                                                           <  2:924175/150‑1 (MQ=255)                                                                                         ttGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                            >  1:315157/1‑149 (MQ=255)                                                                                              gTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                                      >  1:842018/1‑150 (MQ=255)                                                                                                  ccAGCCCTGCCAGCGGGCGCAGCGGGTTTTTTTTCCCCCCCCCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                   <  1:439728/149‑1 (MQ=255)                                                                                                  ccAGCCCTGCCAGAAGGCGCGGAAGGTTACTTTTGCCCCCACCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa                                                    <  2:486271/148‑1 (MQ=255)                                                                                                  ccAGCCCGGCCAGCAGGCGCAGCGGGTTCCTTTCCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                   <  2:281638/149‑1 (MQ=255)                                                                                                      aCCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                             >  1:497697/1‑151 (MQ=255)                                                                                                           ccaccaGGCGCGGCGGGGTTTTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                         <  2:635391/150‑1 (MQ=255)                                                                                                               cgggCGCGGGAGGGTCTTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                                      <  2:507706/147‑1 (MQ=255)                                                                                                                 ggCGCAGCGGGGTCCTTTGCCCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                                   <  1:471213/150‑1 (MQ=255)                                                                                                                        cagGGTTTTTTCCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGa                              <  2:32218/148‑1 (MQ=255)                                                                                                                          ggctccTTTTCCCCCCCCCGAGGCGCCCCGCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCa                           <  2:726368/144‑1 (MQ=255)                                                                                                                             tACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                      >  2:248709/1‑151 (MQ=255)                                                                                                                              aCTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                      >  1:601081/1‑150 (MQ=255)                                                                                                                                      caccacCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg             >  2:108245/1‑151 (MQ=255)                                                                                                                                       accaccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg            >  2:701643/1‑151 (MQ=255)                                                                                                                                       accaccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg            >  1:648292/1‑151 (MQ=255)                                                                                                                                        cccccGCAGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt              <  2:224093/145‑1 (MQ=255)                                                                                                                                              cTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt     >  1:212331/1‑151 (MQ=255)                                                                                                                                                 cGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCa   >  1:266570/1‑150 (MQ=255)                                                                                                                                                  ggCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCaa  >  2:901138/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCA‑GTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAA  >  NC_000913/993768‑994061

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