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%	188.9 / NA	52	T416A (ACC→GCC)	ynfE	putative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/24);  total (28/24)

AGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGC  >  NC_000913/1659166‑1659458                                                                                                                                                     |                                                                                                                                                 aGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTAtgcc                                                                                                                                                <  1:475036/151‑1 (MQ=255)    aaaaCCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTAtgccatt                                                                                                                                             >  2:530471/1‑151 (MQ=255)            cTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTAtgccattaccat                                                                                                                                        >  1:104225/1‑148 (MQ=255)             ttATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATa                                                                                                                                       >  1:504440/1‑148 (MQ=255)              tataTCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGa                                                                                                                                     <  1:708858/149‑1 (MQ=255)               atatCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACg                                                                                                                                  >  2:129404/1‑151 (MQ=255)                      ccAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCg                                                                                                                           <  1:334027/151‑1 (MQ=255)                        aGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCgg                                                                                                                          >  1:580600/1‑150 (MQ=255)                         gggCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTg                                                                                                                        >  1:377878/1‑151 (MQ=255)                         gggCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTg                                                                                                                        >  1:562444/1‑151 (MQ=255)                              gggggCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTg                                                                                                                     <  2:377858/149‑1 (MQ=255)                                   ccACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATc                                                                                                              <  2:133902/151‑1 (MQ=255)                                   ccACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATc                                                                                                              >  1:656895/1‑151 (MQ=255)                                    cacaACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATcc                                                                                                             <  2:300710/151‑1 (MQ=255)                                       aaCGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGTCGGTGCTGGATAAt                                                                                                               >  2:152826/1‑146 (MQ=255)                                            cAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCaaaa                                                                                                     <  1:530419/151‑1 (MQ=255)                                                    cGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAAtc                                                                                             <  2:572221/151‑1 (MQ=255)                                                        gAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCt                                                                                         >  2:652918/1‑151 (MQ=255)                                                         aaCTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTg                                                                                        >  2:734787/1‑151 (MQ=255)                                                         aaCTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTg                                                                                        <  1:106264/151‑1 (MQ=255)                                                             gactgCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCtt                                                                                     <  2:193499/150‑1 (MQ=255)                                                             gactgCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTc                                                                                    <  1:43922/151‑1 (MQ=255)                                                                 gCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGc                                                                                       >  2:477337/1‑144 (MQ=255)                                                                     gcgcTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGAc                                                                            <  1:585824/151‑1 (MQ=255)                                                                      cgcTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGa                                                                             >  1:746637/1‑149 (MQ=255)                                                                         tATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGAc                                                                            >  1:394191/1‑147 (MQ=255)                                                                          aTTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATg                                                                       >  1:31151/1‑151 (MQ=255)                                                                              cTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTAt                                                                   >  2:508610/1‑151 (MQ=255)                                                                              cTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCGTCTTCAGCTGGACAGATGCTAt                                                                   >  1:650300/1‑151 (MQ=255)                                                                                aTGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCg                                                                 >  1:728895/1‑151 (MQ=255)                                                                                    tACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATc                                                              >  2:230254/1‑150 (MQ=255)                                                                                      ccTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATc                                                              <  1:532892/148‑1 (MQ=255)                                                                                      ccTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATg                                                           <  2:532012/151‑1 (MQ=255)                                                                                       cTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATgg                                                          >  2:181447/1‑151 (MQ=255)                                                                                          ttttGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTcc                                                       >  1:250956/1‑151 (MQ=255)                                                                                            ttGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCg                                                      <  1:266783/150‑1 (MQ=255)                                                                                             tGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGc                                                     >  1:471466/1‑150 (MQ=255)                                                                                              gACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCa                                                    <  1:209977/150‑1 (MQ=255)                                                                                                    cAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGAc                                              >  2:12482/1‑150 (MQ=255)                                                                                                     aaTGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACgg                                            >  1:260187/1‑151 (MQ=255)                                                                                                        gTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCa                                          >  1:300675/1‑150 (MQ=255)                                                                                                          cGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCaa                                         <  1:229681/149‑1 (MQ=255)                                                                                                                 aGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACg                                <  1:748206/151‑1 (MQ=255)                                                                                                                   cggcggAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGTCAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGc                              <  2:658635/151‑1 (MQ=255)                                                                                                                     gcggAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGt                            <  1:504645/151‑1 (MQ=255)                                                                                                                       ggAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTcc                          <  2:185978/151‑1 (MQ=255)                                                                                                                          aaTAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCcgcg                       >  1:215298/1‑151 (MQ=255)                                                                                                                                ggCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCaaa                  <  2:230848/150‑1 (MQ=255)                                                                                                                                      cgTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAAc           >  1:732621/1‑151 (MQ=255)                                                                                                                                      cgTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAAc           <  1:57668/151‑1 (MQ=255)                                                                                                                                         gAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTgg        >  1:727924/1‑151 (MQ=255)                                                                                                                                             cGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGAtg     >  1:119241/1‑150 (MQ=255)                                                                                                                                              ggccTATGCCATTTCCATAGAACGCCTGCCGGTGCTGGATAATCCTGTCAAAACGTCAATCTCCTGCTACAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGAtgtg   <  1:260233/149‑1 (MQ=255)                                                                                                                                               aCCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGc  >  2:274864/1‑151 (MQ=255)                                                                                                                                                ccTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGAtgtg   <  2:521533/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 AGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGC  >  NC_000913/1659166‑1659458

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