Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,913	0	G	C	53.2%	54.1 / 76.0	47	G44G (GGC→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (17/5);  new base C (0/25);  total (17/30)
Fisher's exact test for biased strand distribution p-value = 9.61e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.57e-07

AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCT  >  NC_000913/993770‑994057                                                                                                                                                 |                                                                                                                                                  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)                                                                                                                                                 |                                                                                                                                                  AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCT  >  NC_000913/993770‑994057

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