Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,901	0	A	C	54.4%	64.4 / 75.7	46	G48G (GGT→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (16/5);  new base C (0/25);  total (16/30)
Fisher's exact test for biased strand distribution p-value = 2.05e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.22e-10

CCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTG  >  NC_000913/993756‑994050                                                                                                                                                   |                                                                                                                                                       ccGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCaccacc                                                                                                                                                   <  2:302861/151‑1 (MQ=255)  cGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCg                                                                                                                                                  >  2:85467/1‑151 (MQ=255)  cGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCg                                                                                                                                                  <  1:15213/151‑1 (MQ=255)       tCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGc                                                                                                                                              <  1:867354/150‑1 (MQ=255)        ccACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCg                                                                                                                                             <  2:812681/150‑1 (MQ=255)        ccACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCg                                                                                                                                             <  1:497455/150‑1 (MQ=255)               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:674478/1‑151 (MQ=255)                                                                     ggTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                              >  2:487145/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)                                                                                                                                                  caccacCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg  >  2:723479/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                       CCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTG  >  NC_000913/993756‑994050

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