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

AGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAA  >  NC_000913/993773‑994061                                                                                                                                                |                                                                                                                                                     aGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCg                                                                                                                                            >  2:386743/1‑151 (MQ=255) aGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCg                                                                                                                                            >  1:373861/1‑151 (MQ=255)   aCGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGcc                                                                                                                                          >  2:898390/1‑151 (MQ=255)    cgcgCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCca                                                                                                                                         <  1:633211/151‑1 (MQ=255)      cgcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCcac                                                                                                                                        >  1:471190/1‑150 (MQ=255)      cgcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCcac                                                                                                                                        >  2:788069/1‑150 (MQ=255)         aTCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAAc                                                                                                                                    >  1:866036/1‑151 (MQ=255)                             tgtCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTa                                                                                                                 >  1:844089/1‑150 (MQ=255)                                  ttcTGAATTTCCACCCGGGGGTGTGGCCCCCAATAACCCCCCGCCTG‑GGGTGGGGTTTCCAACCCCTCCCAGAGGGCGCAGGGGGGTCTTTTTC‑CCCCCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                  <  1:708078/143‑1 (MQ=255)                                   ccTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                            >  2:678642/1‑149 (MQ=255)                                       aaTTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTgg                                                                                                      >  2:254872/1‑151 (MQ=255)                                       aaTTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTgg                                                                                                      >  1:10344/1‑151 (MQ=255)                                         ttt‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGtt                                                                                                    >  2:28551/1‑151 (MQ=255)                                          ttccacccgggggTGTGGTCCCCGCAAACACACCGCCGGGGTTTGGGTTTCCC‑AGCCCGCCCAGCGGGCGCAGCGGGGTCTTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                        <  1:972690/136‑1 (MQ=255)                                           tccggccgggggTGTGGTCCCCGCAACCACATCGCCTGGGGTTGGCTTTCCC‑AGCCCTGCCAGCGGGCGCAGCGGGTTCTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGtt                                                                                                    <  1:413657/140‑1 (MQ=255)                                           tccagccgggGGTGGGGTCCCCGCAAACACACCCCCGGGGGTTGGCTTTCCC‑AGCCCTGCCAGCAGGCGCAGCAGGGTATTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGtt                                                                                                    <  1:279307/148‑1 (MQ=255)                                            ccacccgggggTGGGGTCCCCGCTAACACATCCCCGGGGGTGGGCTTTCCC‑AGCCCTGCCAGCAGGCGCAGCGGGTTACTTTTG‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCa                                                                                                  <  1:213324/142‑1 (MQ=255)                                             cagccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                   >  1:129110/1‑149 (MQ=255)                                                   gggggggggtgCCCGCTAACACACCGCGTGGGGTGGGCGTTTCC‑AGCCCGGCCAGCAGGCGCAGCGGGTTACTTTTG‑CCCCCACCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGAc                                                                                              <  1:288126/141‑1 (MQ=255)                                                    cGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGat                                                                                           >  1:347192/1‑150 (MQ=255)                                                        gtggtgCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtttt                                                                                      >  1:113851/1‑151 (MQ=255)                                                        gtggtgCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtttt                                                                                      >  1:902285/1‑151 (MQ=255)                                                        gtggtgCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtttt                                                                                      >  2:831331/1‑151 (MQ=255)                                                        gtggtgCCCGCAAACCCTCCGCCGGGGGTTGGCGTTCCC‑ACCCCTGCCAGCGGGCGCAGCGGGTTTCTTTCC‑CCCCCCCCGCTGCGCCCCACCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt                                                                                       <  1:186762/150‑1 (MQ=255)                                                          ggtgCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc                                                                                    >  2:1032941/1‑151 (MQ=255)                                                               cctcaaaCAAATCGCCGGGGTTGGGCTTTCCC‑AGCCCTCCCAGCAGGCGCAGGAGGGTTTTTTCG‑CCCCCCCCGCTGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTa                                                                               <  1:423942/146‑1 (MQ=255)                                                               ccGCAAACACCCCGCCTGGGTTGGGCTTTCCC‑AGCCCTGCCAGCGGGCGCAGCGGGTTTTTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTa                                                                               <  1:677218/151‑1 (MQ=255)                                                                 gCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAAt                                                                             >  2:572029/1‑151 (MQ=255)                                                                    aacccttcccttggggTTGGCTTTCCC‑AGCCCTCCCGGCGGGCCCAGAAGGTTTTTTTTC‑CCCCCCCGGGTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                            <  2:988775/135‑1 (MQ=255)                                                                    aaCACCCCGCCGGGGGTTGGCTTTCCC‑ACCCCTCCCGGCGGGCGCAGGGGGTTATTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                          <  1:100970/151‑1 (MQ=255)                                                                     acacATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                         >  1:888685/1‑151 (MQ=255)                                                                            gCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCCGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTAc                                                                   >  2:122459/1‑150 (MQ=255)                                                                                  ggTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAtt                                                             >  2:743593/1‑150 (MQ=255)                                                                                    ttGGCTTTCCC‑CGCCCTCCCGGCGGGCGCAGAGGGGTTTTTTTC‑CCCCCCCCGCTGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGt                                                                           <  2:425436/134‑1 (MQ=255)                                                                                    tgggCTTTCCC‑AGCCCAGCCAGCGGGCCCAGGAGTGTTCTTTCC‑CCCCCCCCGCGGGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                         <  2:1040250/134‑1 (MQ=255)                                                                                           ttcccagcccgcccAGCGGGCGCAGCGGGTTATTTTTC‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGa                                                           <  1:603381/132‑1 (MQ=255)                                                                                              c‑AGACCTGCAATCAAGCGCAGCAGGGTACTGTTG‑CCACCACCGCTGCTCACCACCACCGCCCCAAACTGACCTGCCGGAATATGTAAATCCAGTTGGGTCAGGACGATATTTTCCGCGCAATGTTTGCATGAAGCAACTAGCAACAATgg                                                 >  2:749816/1‑150 (MQ=255)                                                                                                agcctgcccAGCAGGCCCGGCAGGTTTCTTTCG‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                               <  2:138075/144‑1 (MQ=255)                                                                                                accccTGCCGGCAGGCGCAGCAGGTTTCTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                 <  1:1010608/146‑1 (MQ=255)                                                                                                aGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGCCGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                               >  2:11725/1‑151 (MQ=255)                                                                                                      gCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGCCGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGccc                                          >  2:1020296/1‑150 (MQ=255)                                                                                                          gcagGCGAAGC‑GGGTTTTTTTGCCCCCCACCGCGGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                                      <  1:419300/150‑1 (MQ=255)                                                                                                             ggCGCAGCAGGTTTTTTTGG‑CCCCCCCCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTg                                        <  1:742726/145‑1 (MQ=255)                                                                                                                 cagcagGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                               >  1:729016/1‑150 (MQ=255)                                                                                                                      ggttttttttC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAg                             <  1:788365/141‑1 (MQ=255)                                                                                                                         tACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                      >  2:844734/1‑151 (MQ=255)                                                                                                                            ttttC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                     <  1:386087/149‑1 (MQ=255)                                                                                                                             tttgcccacccCCGCTGCCCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                      <  2:359509/143‑1 (MQ=255)                                                                                                                                  ccaccaCCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg               >  2:259825/1‑150 (MQ=255)                                                                                                                                    accaccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg            >  1:461198/1‑151 (MQ=255)                                                                                                                                       accGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct         >  2:130724/1‑151 (MQ=255)                                                                                                                                            tGCGGCCCACCACCGCCACCAACTGACCTGCCGGAATATGTAAGTCCAGTTGGTTCCGCACGACATTTTCCGCGTACTGTTTGCTCCCTGCATTGAGCAACAATGGCGTCCCCTGGTTCCGTCGGGCAGTATTCATCCCGTGGCCTCCTTc    >  2:117689/1‑151 (MQ=255)                                                                                                                                               ggCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCaa  >  1:149388/1‑150 (MQ=255)                                                                                                                                               ggCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCaa  >  2:997598/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                     AGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAA  >  NC_000913/993773‑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: ‑