Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,659,314	A→G	T416A (ACC→GCC)	ynfE →	putative selenate reductase, periplasmic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,659,314	0	A	G	100.0%	173.9 / NA	49	T416A (ACC→GCC)	ynfE	putative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/23);  total (26/23)

CAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATG  >  NC_000913/1659168‑1659455                                                                                                                                                   |                                                                                                                                              cAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTAtgccat                                                                                                                                           >  1:139601/1‑151 (MQ=255)        cGCATATATCAGCCAGGGCTGGGGGCCACAACGCCAGGCAACCGGCGAACTGACTGCACGCGCCATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCATGAATCGACCTATGCCATTCCCATa                                                                                                                                    <  1:758769/151‑1 (MQ=255)          cTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTAtgccattacca                                                                                                                                      <  1:758387/147‑1 (MQ=255)          cTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTAtgccattacca                                                                                                                                      >  2:758444/1‑147 (MQ=255)             atatCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACg                                                                                                                               <  2:674191/151‑1 (MQ=255)                tCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCt                                                                                                                            >  2:623516/1‑151 (MQ=255)                 cTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTg                                                                                                                           >  2:775638/1‑151 (MQ=255)                      aGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCATTGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCgg                                                                                                                       >  2:11423/1‑150 (MQ=255)                             ggggCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGAt                                                                                                               >  2:804477/1‑151 (MQ=255)                              gggCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGAt                                                                                                               >  2:35986/1‑150 (MQ=255)                                 ccACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCTGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATc                                                                                                           >  2:142567/1‑151 (MQ=255)                                  cacaACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATcc                                                                                                          <  1:806445/151‑1 (MQ=255)                                    caACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCgg                                                                                                        <  2:423844/151‑1 (MQ=255)                                            ggCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACg                                                                                                <  1:245716/151‑1 (MQ=255)                                              cAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGt                                                                                               <  1:343788/150‑1 (MQ=255)                                               aaaCGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTc                                                                                              >  2:802225/1‑150 (MQ=255)                                                 aCGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTAGTTTGACGGACAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACAATAGAACGCCTGCCGGTGCTGGATAAACCGGTCAAAACGTCAAt                                                                                           <  1:172795/151‑1 (MQ=255)                                                              tgCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGc                                                                              >  1:36211/1‑151 (MQ=255)                                                              tgCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGc                                                                              >  2:320389/1‑151 (MQ=255)                                                                  cgcgCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTgg                                                                           >  1:481901/1‑150 (MQ=255)                                                                           gCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCt                                                                  <  1:329249/150‑1 (MQ=255)                                                                                gTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGa                                                             >  1:306208/1‑150 (MQ=255)                                                                                 ttACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATc                                                           >  2:19073/1‑151 (MQ=255)                                                                                  tACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATc                                                           >  1:504820/1‑150 (MQ=255)                                                                                   aCCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCATGCTTCAGCTGGACAGATGCTATCGAACAt                                                         >  2:523147/1‑151 (MQ=255)                                                                                      tATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAAAGCCTGCCGGCGCTGGATAATCCGGTCAACACGGCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGt                                                      >  1:183544/1‑151 (MQ=255)                                                                                       aTTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGt                                                      <  1:593913/150‑1 (MQ=255)                                                                                             aCGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCaa                                                >  1:792576/1‑150 (MQ=255)                                                                                               gggCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATg                                             >  1:804165/1‑151 (MQ=255)                                                                                                 gCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGAc                                           <  1:134687/151‑1 (MQ=255)                                                                                                  cAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACg                                          >  1:191077/1‑151 (MQ=255)                                                                                                  cAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACg                                          <  1:654546/151‑1 (MQ=255)                                                                                                   aaTGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACgg                                         >  2:505312/1‑151 (MQ=255)                                                                                                    aTGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGc                                        <  2:460922/151‑1 (MQ=255)                                                                                                     tGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGc                                        >  2:331523/1‑150 (MQ=255)                                                                                                       tCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAAt                                     >  2:532227/1‑151 (MQ=255)                                                                                                        cGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAAt                                     <  1:379942/150‑1 (MQ=255)                                                                                                           cATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCcgc                                 >  2:409295/1‑151 (MQ=255)                                                                                                                 cggcggAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGc                           <  1:140397/151‑1 (MQ=255)                                                                                                                 cggcggAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGc                           >  1:381086/1‑151 (MQ=255)                                                                                                                     ggAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTcc                       <  1:748732/151‑1 (MQ=255)                                                                                                                        aaTAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCcgcg                    <  1:423674/151‑1 (MQ=255)                                                                                                                         aTAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCgg                   <  1:335296/151‑1 (MQ=255)                                                                                                                           aGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGc                  >  1:705926/1‑150 (MQ=255)                                                                                                                           aGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGc                  >  1:72367/1‑150 (MQ=255)                                                                                                                                cgcgcgTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGAt            <  1:177542/151‑1 (MQ=255)                                                                                                                                 gcgcgTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATa           <  2:93049/151‑1 (MQ=255)                                                                                                                                  cgcgTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATaa          <  1:519922/151‑1 (MQ=255)                                                                                                                                     gTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACt       <  2:147412/151‑1 (MQ=255)                                                                                                                                     gTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACt       <  1:819544/151‑1 (MQ=255)                                                                                                                                       gAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTgg     >  2:775245/1‑151 (MQ=255)                                                                                                                                           cGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGAtg  <  1:212810/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                              CAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATG  >  NC_000913/1659168‑1659455

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