Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,917	0	A	C	45.3%	73.3 / 12.2	53	V43G (GTG→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (22/7);  new base C (0/24);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 3.50e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993768‑994063                                                                                                                                                       |                                                                                                                                                   gTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCccac                                                                                                                                                   >  1:344386/1‑151 (MQ=255)   aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCccacc                                                                                                                                                  >  2:625289/1‑150 (MQ=255)       gACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGc                                                                                                                                             >  2:325112/1‑151 (MQ=255)          gcgcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCca                                                                                                                                           <  1:706198/150‑1 (MQ=255)             cATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACaa                                                                                                                                        >  2:491338/1‑150 (MQ=255)                 cTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTg                                                                                                                                    >  1:132696/1‑150 (MQ=255)                   gAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGAcc                                                                                                                                 >  1:655853/1‑151 (MQ=255)                      aCATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGc                                                                                                                              >  1:326731/1‑151 (MQ=255)                                  tttctcCCTGATTTCCAGCGGGGGGTGGGGCCCCCGCAAACCCCCCCCCGGGGGTGGGCTTTCCCAGCCCGGCCAGCAGGCGCAGCAGGGTACTTTCC‑CCACCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTa                                                                                                                   <  1:923948/148‑1 (MQ=255)                                       tCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                              >  1:127687/1‑150 (MQ=255)                                        ccTGATTTTCAGCCGGCGGTGGGGGCCCCGCAAACCCTCCGCCTGGGTTGGGCTTTTCCAGCCCTCCCAGCGGGCGCAGCAGGTTATTTTTC‑CCCCCACCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                               <  2:762537/148‑1 (MQ=255)                                              ttccagccgggGGTGGGGTCCCCGCAAACCCACCGCCGGGGGTGGGCGTTCCCAGCCCGGCCAGCAGGCGCAGCAGGGTTTTTTTC‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTgg                                                                                                        <  2:650254/146‑1 (MQ=255)                                                 cagccgggGGTGTGGTGCCGCCAACCACACGGTGTGCGTTTGGCTTTCCCACCCCCCCAAGCGGGCCAAGCGGGTTCTTTTTC‑CCACCCCCCCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAg                                                                                                   <  2:673359/151‑1 (MQ=255)                                                  agccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg                                                                                                  >  2:878056/1‑151 (MQ=255)                                                    ccggggttggggtcCCCGCAACCCCACCGCCGGCGGTGGGCTTTCCCAGCCCTCCCAGCAGGCGCAGCGGGTTTCTTTTC‑CCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAg                                                                                                   <  2:100225/138‑1 (MQ=255)                                                         ggTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                           >  2:1028333/1‑151 (MQ=255)                                                            gtggtgCCCGCTAACACATCGCCGGGGTTTGGCTTTCCCAGCCCGCCCAGCAGGCGCAGGGGGGT‑CTTTTGCCCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt                                                                                         <  2:432138/150‑1 (MQ=255)                                                              ggcccccGCAAACACTCCGCCTGCGTTGGGCTTTCCCACCCCGCCCAGCGGGCGCAGCAGGTTTTTTTTG‑CCCCCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt                                                                                         <  1:179258/144‑1 (MQ=255)                                                                    cGCTACCACATCGCCGGCGGTTGGCGTTCCCAGCCCTGCCAGCGGGCGCAGCAGGGTCTTTTTG‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTaa                                                                                <  2:327784/151‑1 (MQ=255)                                                                        accactcccccgggggtgggtgTTCCCACCCCGGCCAGCGGGCCCAGCAGGGTTTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                            <  1:1043936/137‑1 (MQ=255)                                                                            caccGCCTGGGGTGGGCTTTCCCAGCCCGGCCAGCGGGCGCAGCAGGGTATTTTTG‑CCACCCCCGCTGCGGCCTACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                           <  2:89784/145‑1 (MQ=255)                                                                               ccccggggtttGGCGTTCCCGGACCTGCCAGCGGGCGCGGCGGGGTTTTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTa                                                                      <  2:574052/141‑1 (MQ=255)                                                                               cGCCGGGGGTGGGCTTTCCCAGCCCGGCCAGCAGGCGCGGCAGGTTCTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                       <  2:300685/149‑1 (MQ=255)                                                                                    gggtttGGCGTTCCCAGCCCGGCCAGCGGGCGCGGCAGGGTCCTTTTG‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                                  <  2:282478/145‑1 (MQ=255)                                                                                      gtttGGCTTTCCGAGCCCGGCCAGCGGGCGCAGCGGGTTATTTTCC‑CCCCCCCCGTTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACt                                                                    <  2:169315/143‑1 (MQ=255)                                                                                      ggTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                              >  1:686460/1‑151 (MQ=255)                                                                                                cccaccccggCCAGCAGGCGCAGCGGGTTTCTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGa                                                             <  1:140841/133‑1 (MQ=255)                                                                                                 ccAGCCCTGCCAGAGGGCGCAGCGGGTTTTTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATg                                                    <  2:399329/150‑1 (MQ=255)                                                                                                      cccgCCAGCGGCCGCGGGGGGGTTTTTTTC‑CCCCCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGc                                              <  1:369816/148‑1 (MQ=255)                                                                                                       cTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGcc                                             >  2:978121/1‑151 (MQ=255)                                                                                                       cTGCAAGCGGGCGCGGAAGGTTATTTTGC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                                 <  2:143130/147‑1 (MQ=255)                                                                                                          ccAGCAGGCGCAGCAGGTTATTTTTG‑CCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                           <  2:166664/150‑1 (MQ=255)                                                                                                           cagcggGCGCAGCGGGTTCTTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTg                                          <  1:778617/150‑1 (MQ=255)                                                                                                               agGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                      >  2:522338/1‑150 (MQ=255)                                                                                                               agGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                      >  1:876511/1‑150 (MQ=255)                                                                                                                ggCGCAGCAGGGTTTTTTCG‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAg                                    <  2:81561/151‑1 (MQ=255)                                                                                                                ggCGCAGCAGGGTCTTTTTC‑CCCCCCCCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                                     <  2:30009/150‑1 (MQ=255)                                                                                                                ggCGCAGCAGGGTCTTTTCC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAg                                    <  2:122948/151‑1 (MQ=255)                                                                                                                  cgcagggggttctttttC‑CCACCCCCGCTGCGCCCCACCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGa                                   <  1:267277/137‑1 (MQ=255)                                                                                                                   gcagcaGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                 >  1:146347/1‑151 (MQ=255)                                                                                                                    caccaggttttttttC‑CCCCCCCCGCTGCGCCCCCCCTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                      <  1:625179/134‑1 (MQ=255)                                                                                                                      gcagattACTTTTC‑CCACCACCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGa                                   <  2:132716/140‑1 (MQ=255)                                                                                                                        gggttctttttG‑CCCCCCCCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCa                             <  2:806137/143‑1 (MQ=255)                                                                                                                          ggTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGt                           >  1:585674/1‑150 (MQ=255)                                                                                                                           gTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt                         >  2:264974/1‑151 (MQ=255)                                                                                                                              cttttccccccccccGCTGCCCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCa                             <  1:221492/132‑1 (MQ=255)                                                                                                                                tttG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                     >  2:940520/1‑150 (MQ=255)                                                                                                                                        ccaccGCGCCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGc             <  1:976893/151‑1 (MQ=255)                                                                                                                                          accGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct           >  2:810549/1‑151 (MQ=255)                                                                                                                                           ccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctc          >  2:1022357/1‑151 (MQ=255)                                                                                                                                            cGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcc         >  1:656718/1‑151 (MQ=255)                                                                                                                                               tGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTc      >  2:392593/1‑151 (MQ=255)                                                                                                                                                gcgcCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcc         <  2:90673/147‑1 (MQ=255)                                                                                                                                                 cgcccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt       <  2:496473/145‑1 (MQ=255)                                                                                                                                                  ggCCCACCACCGCCACAAACCGCCCTGCCATAATAAGTAAAACGAGTTTGTGCACGACGATATTCTCGGCGTAATGTCTGCCTGCTGCGTTGAGCAACAATGTCCTTCACTGCTTCAGACGAGCAGTATTCACACCGAGGCCTCCTTc      >  1:653249/1‑148 (MQ=255)                                                                                                                                                   gCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAt  >  2:121312/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                   GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993768‑994063

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: ‑