Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,917	0	A	C	61.7%	23.4 / 57.5	47	V43G (GTG→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (16/2);  new base C (0/29);  total (16/31)
Fisher's exact test for biased strand distribution p-value = 1.02e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.04e-08

AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATGAT  >  NC_000913/993770‑994066                                                                                                                                                      |                                                                                                                                                      aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGCGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑Cccacc                                                                                                                                                     >  1:685847/1‑150 (MQ=255)       cgcgCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGcc                                                                                                                                               >  2:626072/1‑150 (MQ=255)       cgcgCATCCTGAAACATCATTCGCGTGTCTGCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGcc                                                                                                                                               <  2:1001022/150‑1 (MQ=255)                        catTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCgg                                                                                                                              >  2:26816/1‑150 (MQ=255)                          ttCGCTTTTCTTCCTAATTTCCCGCCGGGGGTGTGGTGCCCGCTAACACCTCGCCGGCGTTTGGCTTTTCC‑AGCCCGGCCAGCGGGCGCGGCGGGTTTTTTTTCCCCCCCCCGCGGCGC‑CCCCCCACCGCCACAAACTGACCTGcc                                                                                                                                <  1:666772/146‑1 (MQ=255)                          tgcgctTGTCTTCCTAATTTTCAGCCGGCGGTGGGGTGCCCGAAACCACCCCGCCGGGGGTGGGCTTTCCC‑AGCCCTGCCAGCGGGCGCAGCGGGTTTTTTTTCCCCCCCCCGCGGCGC‑CCCCCCACCGCCACAAACTGACCTGCCg                                                                                                                               <  1:1044361/145‑1 (MQ=255)                            cgggtttctcCCTAAATTTCACCCAGCGGTGTGGTCCCCGCTAACACATCGCCGGGGGTTGGCTTTTCC‑AGCCCTGCCAGCAGGCGCAGCGGGTTCTTTTTCCCCCCCCCGCTGCGC‑CCCCCCCCCGCCACAAACTGACCTGCCGGAat                                                                                                                           <  1:977922/143‑1 (MQ=255)                               gtgtCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                       >  1:566476/1‑150 (MQ=255)                                ttccTTCCTGAGTTCCAGCGAGCGGTGTGGTGCCGGCTACCGCCCCGCGTGCGGTTGGCTTTCCC‑ACCCCAGCCAGCAGGCGCAGCGGGGTTCTTTGGCCCCCACCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTa                                                                                                                      <  1:745707/147‑1 (MQ=255)                                     tCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                                 >  2:657683/1‑150 (MQ=255)                                         gAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                            >  1:239183/1‑151 (MQ=255)                                            ttccacccgggggTGTGGTCCCCGCAACCACACCCCCTGGGTTTGGCTTTCCC‑AGCCCGGCCAGCAGGCGCAGCGGGGTTCTTTTCCCCCCCCCGCGGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                            <  1:587303/137‑1 (MQ=255)                                             ttCACCCGGGGGTGGGGTCCCCGCTAACACCCCGCCTGGGGTTGGGTTTCCC‑AGCCCTGCCAGCGGGCGCAGCGGGTTTTTTTTGCCCCCACCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGtt                                                                                                         <  1:46944/150‑1 (MQ=255)                                                  ccaggGGTGTGGTCCCCGCAAACCCACCGCCTGGGGTTGGCGTTCCC‑AGCCCTGCCAGCGGGCGCAGCAGGGTACTTTTCCCCCCCCCGCGGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGa                                                                                                    <  1:31263/150‑1 (MQ=255)                                                     ggggTGTGGTCCCCGCAAACCCACCGCCTGCGGTTGGCTTTTCC‑AGCCCCGCCAGCAGGCGCAGCGGGTTTCTTTTGCCCCCCCCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg                                                                                                     <  1:665156/144‑1 (MQ=255)                                                       ggTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                              >  2:487145/1‑151 (MQ=255)                                                       ggTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                              >  2:674478/1‑151 (MQ=255)                                                                cccGCAAACGCACCGCCTGGGGTGGGCGTTCCC‑ACCCCCGCCAGCGGGCGCAGGAGGTTCTTTTTGCCCCCCCCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgc                                                                                       <  2:373663/149‑1 (MQ=255)                                                                 ccGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                                     >  2:770741/1‑150 (MQ=255)                                                                      aaCACCCCGCCGGGGGTTGGCTTTCCC‑ACCCCTGCCAGCGGGCGCAGCGGGGTTTTTTTCCCCCCCCCCCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                                 <  1:398239/149‑1 (MQ=255)                                                                         acaTCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc                                                                            >  1:753668/1‑151 (MQ=255)                                                                               cccgggttgggCTTTCCC‑AGCCCCCCCAGCAGGCGCGGCAGTTTTTTTTCCCCCCCACCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAAt                                                                                  <  2:613608/130‑1 (MQ=255)                                                                               ccTGGGGTGGGCTTTCCC‑AGCCCTCCCGGCAGGCGCAGCAGGTTTTTTTTGCCCCCCCCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                          <  1:1035566/147‑1 (MQ=255)                                                                                 tGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                                     >  1:1042088/1‑150 (MQ=255)                                                                                        ggCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                              >  2:122686/1‑150 (MQ=255)                                                                                         gCGTTTCCAAAACCTGCCGGCAGGCGCAGAGGGGTCCTTTTCCCCCCCCCGGTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                              <  1:822985/150‑1 (MQ=255)                                                                                            ttccc‑AGCCCTGCCGGCGGGCGCGGAGGGGTTTTTTTCCCCCCCCCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                                           <  2:1044382/146‑1 (MQ=255)                                                                                            ttccc‑AGCCCCGCCAGCGGGCGCAGCGGGTTTTTTTCCCCCCCCCCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAg                                                               <  2:1084108/142‑1 (MQ=255)                                                                                             ttcccagcccgcccAGCAGGCGCAGCAGGTTCTTTTTCCCCCCCCCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaa                                                            <  1:506732/136‑1 (MQ=255)                                                                                               cc‑AGCCCTGCCGGCAGGCGCAGCAGGTTACTTTTCCCCCCCCCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                        <  1:176426/149‑1 (MQ=255)                                                                                                  agcccTGCCAGCGGGCGCAGCGGGTTTTTTTCCCCCCCCCCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGc                                                     <  2:701341/147‑1 (MQ=255)                                                                                                    aCCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                                  >  1:845576/1‑151 (MQ=255)                                                                                                     cccgCCAGCAGGCGCAGCAGGGTACTTTTCCCCCCCCCGCTGCGG‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                                  <  1:352485/147‑1 (MQ=255)                                                                                                     ccccccAGCGGGCGCAGCGGGTTCTTTTTCCCCCCCCCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                                    <  2:71551/144‑1 (MQ=255)                                                                                                            gcagGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGtt                                          >  1:368001/1‑151 (MQ=255)                                                                                                               ggggcagggggttttttttCCCCCCCCCGCTGGGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                                        <  1:678960/135‑1 (MQ=255)                                                                                                                 cgcAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                                     >  2:542270/1‑151 (MQ=255)                                                                                                                   actggaggTACTTTGTCCCTCCCCGCCGCGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGa                                      <  2:1116178/143‑1 (MQ=255)                                                                                                                          ttatttttCCCCCCCCCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTa                             <  1:1064645/146‑1 (MQ=255)                                                                                                                           ttcTTTTGCCCCCCCCGCTGCGC‑CCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTa                             <  1:2900/147‑1 (MQ=255)                                                                                                                              tttcGCCCCCCCCGCTGCGG‑CCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt                            <  1:41770/147‑1 (MQ=255)                                                                                                                               ttccccccccccGCGGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                        <  2:228714/140‑1 (MQ=255)                                                                                                                                   ccccccccGCTGCGG‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt                   <  2:678830/145‑1 (MQ=255)                                                                                                                                    caccacCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg                  >  2:723479/1‑151 (MQ=255)                                                                                                                                      ccaccGCCGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc                     <  2:669202/146‑1 (MQ=255)                                                                                                                                           gCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct           >  1:634938/1‑151 (MQ=255)                                                                                                                                            cTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcc            <  1:312617/149‑1 (MQ=255)                                                                                                                                                 gcccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt          <  1:289794/145‑1 (MQ=255)                                                                                                                                                  ccccACCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct           <  2:79307/144‑1 (MQ=255)                                                                                                                                                   ccccccGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTc         <  1:194301/146‑1 (MQ=255)                                                                                                                                                      accaccGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATGAt  >  2:799919/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATGAT  >  NC_000913/993770‑994066

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