Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,917	0	A	C	45.8%	65.7 / 15.7	49	V43G (GTG→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (22/4);  new base C (0/22);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 5.67e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATGAT  >  NC_000913/993773‑994066                                                                                                                                                   |                                                                                                                                                      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)                                                                                                                                                ccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt          <  1:1108319/145‑1 (MQ=255)                                                                                                                                                  caccacCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATGAt  >  2:931875/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                      AGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATGAT  >  NC_000913/993773‑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: ‑