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

GCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993771‑994063                                                                                                                                                     |                                                                                                                                                   gCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCcaccac                                                                                                                                                <  2:879662/151‑1 (MQ=255)       gcgcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGTTTGGCGTTTCCAGACCTGCCAGCAGGCGCCGCAGGGTACTTTT‑GCCTCCCCCGCTGCGGccctccccctcc                                                                                                                                            <  1:997792/149‑1 (MQ=255)             ccTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTg                                                                                                                                    >  1:540914/1‑151 (MQ=255)                  aaCATCATTCGCGTTTCTTCCTAATTTC‑CAGCCGGCGGTGGGGTCCCCGCTAACACACCGCTGGGGGTTGGCGTTCCCAGCCCTGCCAGCAGGCGCAGCAGGTTACTTTT‑GCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTg                                                                                                                               <  2:136604/151‑1 (MQ=255)                  aaCATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGAc                                                                                                                                  >  2:181568/1‑148 (MQ=255)                   aCATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGc                                                                                                                              >  1:711928/1‑151 (MQ=255)                    catcatTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGcc                                                                                                                             >  1:206209/1‑151 (MQ=255)                       catTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCgg                                                                                                                           >  2:66587/1‑150 (MQ=255)                        atTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGa                                                                                                                          >  2:614609/1‑150 (MQ=255)                         ttCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAat                                                                                                                        >  1:1001524/1‑151 (MQ=255)                           cgcgTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAatat                                                                                                                      >  1:1092966/1‑151 (MQ=255)                                   ttcttaatttt‑CACCCGGGGGTTGGGTCCCCGCTACCCCACCGCCTGGGTTTGCCTTTCCCAGCCCTGCCAGCGGGCGCAGCGGGTTTTTTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                               <  2:941444/142‑1 (MQ=255)                                   ttCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                              >  2:594357/1‑151 (MQ=255)                                    tccctaaTTTCCAGCCGGGGGTGGGGTCCCCGCAAACACCTGCCCTGGGGTTGGCTTTCCCAGCCCGGCCGGCGGGCGCGGC‑GGGTTCTTTTGGCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAatat                                                                                                                      <  2:398697/138‑1 (MQ=255)                                     ccTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                              >  2:1191679/1‑149 (MQ=255)                                          atttccgcccggcGGTGTGGCCCCCAAAAACACTTCGCCGGGGGTGGGCTTCCCGAGCCCTCCCAGCGGCGGCAGCAGGTTTCTTTC‑CCCACCCGCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                                <  1:397924/132‑1 (MQ=255)                                          atttccagccgggGGTGGGGGCCCCGCAAACCCCCCGCCGGGGGTGGGCTTTCCCAACCCGGCCAGCGGGCGCGGCGGGTTCTTTTT‑CCCCCCCCCGCGGGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                          <  2:225893/144‑1 (MQ=255)                                           tttccagccgggGGTGTGGCCCCCTCTAACACACCGCGGGGGGTGGGCTTTCCCAGCCCTCCCAGCGGGCGCAGGAGGTTTTTTTC‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                         <  1:819456/145‑1 (MQ=255)                                               cagccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                     >  2:810508/1‑149 (MQ=255)                                               cagccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                     >  1:201526/1‑149 (MQ=255)                                                  ccgggggTGTGGTGCCCCTTAACCCACCGCCGGGGGTTGGGTTTTCCAGCCCGCCCGGCAGGCGCAGAGGGGTTTTTTT‑GCCCCCACCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                     <  2:107324/141‑1 (MQ=255)                                                   cagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                               >  1:1015628/1‑151 (MQ=255)                                                       ggtggggccccctgtAACACTTCTCCGGCGGTTGGTTTTTCAGGCCCTGCCGGGAGGCGCAGCAGGTTACTTTT‑CCCCCCCCCCCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGAc                                                                                                <  2:318096/132‑1 (MQ=255)                                                            ggtcCCCGCAAACACACCGCTGGCGGTGGGCTTTCCCAGCCCTGCCGGCGGGCGCGGCGGGTTATTTTT‑CCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtttt                                                                                        <  2:131950/149‑1 (MQ=255)                                                                ccgccaaaccccccGCCTGGGTTGGGCTTTCCCAGCCCGCCCAGCGGGCGCAGAAGGTTTTTTTC‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgcg                                                                                   <  1:456484/137‑1 (MQ=255)                                                                cccGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                                  >  2:1010886/1‑151 (MQ=255)                                                                      aaCACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                            >  2:819067/1‑151 (MQ=255)                                                                      aaCACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                            >  1:87976/1‑151 (MQ=255)                                                                           aTCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGTCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                       >  2:629615/1‑151 (MQ=255)                                                                            tCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                       >  2:982763/1‑150 (MQ=255)                                                                              gCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACt                                                                    >  2:100188/1‑151 (MQ=255)                                                                                               ccAGACCTCCCAGCAGGCGCAGCGGGTTATTTTT‑CCCACCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATg                                                    <  2:112797/150‑1 (MQ=255)                                                                                                cAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGc                                                  >  2:461717/1‑151 (MQ=255)                                                                                                 aGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                                 >  1:219319/1‑151 (MQ=255)                                                                                                   cccTGCCAGCGGGCGCAGCGGGTTTTTTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGt                                                <  2:536645/149‑1 (MQ=255)                                                                                                   cccTGCCAGCAGGCGCAGCGGGTTACTTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                               <  2:72237/150‑1 (MQ=255)                                                                                                    cctgcaggggggCGCGGCGGGTTCCTTTT‑CCCCCCTCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGc                                              <  1:516781/141‑1 (MQ=255)                                                                                                    cctccacgagggCGCAGCGGGTTATTTTT‑GCCCCCCCCGTGGCCGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGc                                                  <  1:1030306/137‑1 (MQ=255)                                                                                                    ccTCCCAGCGGGCGCAGCGGGGTCTTTTT‑CCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGt                                                <  2:499799/149‑1 (MQ=255)                                                                                                     cTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGcc                                             >  1:956289/1‑151 (MQ=255)                                                                                                      tGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGccc                                            >  2:591008/1‑151 (MQ=255)                                                                                                              ggCGCAGCGGGGTTCTTTT‑CCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTgg                                         <  2:949447/146‑1 (MQ=255)                                                                                                              ggCGCAGCGGGGTTCTTTC‑CCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTgg                                         <  1:966798/146‑1 (MQ=255)                                                                                                                cgcAGGGGGGTCTTTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                                     <  2:422521/148‑1 (MQ=255)                                                                                                                             tttt‑GCCACCCCGCCGGCCGCCCCCACCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCa                      <  1:421640/150‑1 (MQ=255)                                                                                                                               tt‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                   >  2:304856/1‑151 (MQ=255)                                                                                                                                   ccccccccccTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                   <  1:3280/139‑1 (MQ=255)                                                                                                                                     accaccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg              >  2:338497/1‑151 (MQ=255)                                                                                                                                         ccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct           <  2:655426/150‑1 (MQ=255)                                                                                                                                             tGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt       >  1:692644/1‑150 (MQ=255)                                                                                                                                                 cccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAt  <  1:54334/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993771‑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: ‑