Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,913	0	G	C	42.0%	76.2 / 16.7	51	G44G (GGC→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (24/5);  new base C (0/21);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 1.01e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTC  >  NC_000913/993766‑994059                                                                                                                                                      |                                                                                                                                                   cGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGccc                                                                                                                                                 >  1:653891/1‑151 (MQ=255)  ggTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGccc                                                                                                                                                 >  2:118262/1‑150 (MQ=255)        aGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCg                                                                                                                                          <  2:610674/151‑1 (MQ=255)           cgcgCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCca                                                                                                                                       <  2:266844/151‑1 (MQ=255)               cATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCcaca                                                                                                                                     >  2:595600/1‑149 (MQ=255)                   cTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTg                                                                                                                                >  1:339950/1‑150 (MQ=255)                                        ttCCTAATTTCCGGCCAGGGGGGTGGCCCCCCCAACCACACCCCGTGGGGTGGGGTTTC‑CCAGCCCGGCCAGCGGGCGCAGCAGGTTTTTTTCC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTaaa                                                                                                             <  1:608589/148‑1 (MQ=255)                                        tccataaaTTCCGGCCGGGGGTGTGGCGCCGGTAAACACTCCGCTTGC‑GTTGGCCTTTCCCACCCTTGCCCGCCGGCGAAGGAGGGTTTTTTCG‑CCACCACCGGTGGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                           <  1:410205/144‑1 (MQ=255)                                          ccgaaatataagccgggGGGGTGGGCCCCGCAAACCCACCGCCGGGGTTTGGGTTTC‑CCGGCCCGGCCGGCGGCCGCAGCGGGTTTTTTTGG‑CCACCACCTCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTaa                                                                                                              <  2:381059/135‑1 (MQ=255)                                          ccTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                          >  1:627716/1‑149 (MQ=255)                                          ccTAAATTCCGGCCGGGGGTGTGGCGCCCGCTAACACATCCCCTGCGGTTGGTTTTC‑CCAGCCCGGCCAGCAGGCGCAGCGGGGT‑CTTTTGCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCa                                                                                                         <  1:556845/150‑1 (MQ=255)                                              attttCAGCCGGGGGTGTGGTGCCCGCAAACACCTCGCCTGGGTTTGGCTTTC‑CCAGCCCTCCCAGCGGGCGCAGCGGGTTTTTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGt                                                                                                       <  1:327844/146‑1 (MQ=255)                                              atttcccgccaggGGTGTGGCCCCCGCAAACACCTCGCCGGGGTTTGGCTTTC‑CCAGCCCTCCCAGCAGGCGCAGCGGGTTTTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                      <  1:938972/142‑1 (MQ=255)                                                tttCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGtt                                                                                                  >  1:1145338/1‑151 (MQ=255)                                                ttccagccgggGGTGTGGTCCCCGCAAACACCCCGCCGGCGTTGGGCTTTC‑CCAGCCCTCCCAGCAGGCGCAGCGGGTTTTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                      <  1:794155/144‑1 (MQ=255)                                                ttccacccgggggTGTGGTGCCCGCAAACCCCCCCCCGGGGTTGGGCTTTT‑CCGGCCCTCCCAGCGGGCGCAGCAGGGTTTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                <  1:611783/126‑1 (MQ=255)                                                ttccacccggcGTTGTGGTCCCCGCTAACACACCGCGTGGGTTGGGCTTTC‑CCAGCCCTGCCAGCGGGCGCAGCAGGTTTTTTTTC‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt                                                                                                   <  1:918691/141‑1 (MQ=255)                                                    agccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg                                                                                              >  1:700446/1‑151 (MQ=255)                                                      ccgggggTGGGGTGCCCCCAACCCCGTCGCCGGGGTTGGGCTTTC‑CCACCCCCGCCGGCGGGCGCAGCGGGTTTCTTTTG‑CCACCCCCGGTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGtt                                                                                                  <  2:648464/140‑1 (MQ=255)                                                      ccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGAc                                                                                            >  2:565272/1‑151 (MQ=255)                                                       caggGGGGTGTTGCCCGCTACCACATCGCCTGCGTTTGGCTTTC‑CCAGCCCCGCCAGCAGGCGCAGCGGGGTTCTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTgg                                                                                                    <  1:84779/142‑1 (MQ=255)                                                          cGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGat                                                                                         >  1:278596/1‑150 (MQ=255)                                                           gggggggtcCCCGCTAACACACCGCCGGCGGTGGGCTTTC‑CCAGCCCTGCCAGCAGGCGCAGCGGGTTATTTTTG‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGata                                                                                        <  1:902858/145‑1 (MQ=255)                                                               tggtgCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTc                                                                                   >  2:1052885/1‑151 (MQ=255)                                                                          aaCACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                        >  2:682532/1‑151 (MQ=255)                                                                               aTCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                   >  2:663437/1‑151 (MQ=255)                                                                                tCGCGTGTGGTTGGGGTTC‑CGAGCCCGGCAGGCAGGCGCAGGTGGGTCCTTGGG‑CCCCCACGTCCGCCACCAGCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTACCGCGTATTGTTTGCtt                                                                   <  2:736902/150‑1 (MQ=255)                                                                                    ctgggttgggCGTTC‑CCGAACCCGCTGGCAGGCGCAGTGGGTTACTTTTG‑CCACCCCCGGGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTAc                                                                 <  1:791423/140‑1 (MQ=255)                                                                                    cTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                              >  1:106566/1‑151 (MQ=255)                                                                                     tGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCCGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                              >  1:215729/1‑150 (MQ=255)                                                                                                 ttcccAGCCCTGCCGGCAGGCCAAGCGGGGTACTTTTC‑CCACCCCCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGca                                                      <  2:559732/144‑1 (MQ=255)                                                                                                     cAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGc                                              >  1:570473/1‑151 (MQ=255)                                                                                                      aaaCCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATCTTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                             >  2:633734/3‑151 (MQ=255)                                                                                                       gcccTGCCAGCGGGCGCAGCGGGTTATTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                               <  1:1178490/146‑1 (MQ=255)                                                                                                         ccgcccAGCGGCCGCAGCGGGGTTTTTTCC‑CCACCCCCGCGCCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGc                                          <  2:1172581/147‑1 (MQ=255)                                                                                                         ccTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGc                                          >  1:207849/1‑151 (MQ=255)                                                                                                              cagcggGCGCAGCAGGTTCTTTTTC‑CCCCCCCCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                       <  1:17402/149‑1 (MQ=255)                                                                                                               agcagGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                                    >  2:266358/1‑151 (MQ=255)                                                                                                                   ggCGCAGCGGGTTTTTTTCC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                                 <  1:399955/150‑1 (MQ=255)                                                                                                                   ggCGCAGCGGGGTTTTTTCG‑CCCCCCCCGGTGCCGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTg                                      <  2:1171422/145‑1 (MQ=255)                                                                                                                   ggCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                  >  1:810532/1‑149 (MQ=255)                                                                                                                          gctgggtcTTTTGCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                              <  1:1004536/140‑1 (MQ=255)                                                                                                                           agGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                          >  1:1177052/1‑149 (MQ=255)                                                                                                                            ggttttttttC‑CCCCCCCCGCTGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGa                            <  2:1048283/140‑1 (MQ=255)                                                                                                                                 cTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                   >  1:418590/1‑150 (MQ=255)                                                                                                                                   tttccccccccccGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGt                       <  1:1127002/134‑1 (MQ=255)                                                                                                                                        ccccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATa                <  1:863694/141‑1 (MQ=255)                                                                                                                                        ccccccccGCTCCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt            <  1:671649/145‑1 (MQ=255)                                                                                                                                           ccaccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGc         >  2:739380/1‑151 (MQ=255)                                                                                                                                                  tGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTc  >  2:1117819/1‑151 (MQ=255)                                                                                                                                                  tGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTc  >  2:1056797/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTC  >  NC_000913/993766‑994059

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