Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,658,755	T→G	G229G (GGT→GGG)	ynfE →	putative selenate reductase, periplasmic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,658,755	0	T	G	100.0%	207.2 / NA	55	G229G (GGT→GGG)	ynfE	putative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (27/28);  total (27/28)

GAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGTAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTAT  >  NC_000913/1658614‑1658901                                                                                                                                              |                                                                                                                                                   gaACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCg                                                                                                                                           <  1:702709/151‑1 (MQ=255)   aCTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGc                                                                                                                                         <  1:881269/151‑1 (MQ=255)      gCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCa                                                                                                                                        <  2:14531/149‑1 (MQ=255)         aCGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCCAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCCCGGGTATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAAACCGGCAGaaa                                                                                                                                    >  2:532391/1‑150 (MQ=255)               gTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCAt                                                                                                                             >  2:825480/1‑151 (MQ=255)               gTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCAt                                                                                                                             <  2:177354/151‑1 (MQ=255)                   aCTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGc                                                                                                                         <  2:332796/151‑1 (MQ=255)                   aCTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGc                                                                                                                         <  2:871410/151‑1 (MQ=255)                        aCCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCggtgg                                                                                                                    <  2:495378/151‑1 (MQ=255)                        aCCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCggtg                                                                                                                     >  2:399440/1‑150 (MQ=255)                              aTGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCa                                                                                                               <  2:507187/150‑1 (MQ=255)                               tGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCAt                                                                                                              >  2:489310/1‑150 (MQ=255)                               tGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATc                                                                                                             >  2:242529/1‑151 (MQ=255)                                ggCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATc                                                                                                             >  1:151998/1‑150 (MQ=255)                                    ccTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACtt                                                                                                         >  1:433016/1‑150 (MQ=255)                                     cTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTAt                                                                                                       >  2:233213/1‑151 (MQ=255)                                      tACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTAtc                                                                                                      >  1:467182/1‑151 (MQ=255)                                       aCAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTAtc                                                                                                      <  1:473549/150‑1 (MQ=255)                                       aCAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTAtc                                                                                                      <  1:421572/150‑1 (MQ=255)                                          gCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTAtcttct                                                                                                  <  1:592988/151‑1 (MQ=255)                                             cTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTg                                                                                                >  1:212310/1‑150 (MQ=255)                                               gcgcAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGa                                                                                               <  2:467170/149‑1 (MQ=255)                                                 gcAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGaaaaa                                                                                           <  1:781552/151‑1 (MQ=255)                                                  cAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAg                                                                                          <  1:266142/151‑1 (MQ=255)                                                   aaaTTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAgc                                                                                         >  1:615298/1‑151 (MQ=255)                                                   aaaTTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTATTGaaaaa                                                                                           >  2:283845/1‑149 (MQ=255)                                                       ttCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTATGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAgcgcg                                                                                      <  1:298866/150‑1 (MQ=255)                                                        tCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAgcgcgc                                                                                     <  2:812518/150‑1 (MQ=255)                                                                 aTGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGa                                                                           <  1:758209/151‑1 (MQ=255)                                                                    cccTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACg                                                                        >  2:454726/1‑151 (MQ=255)                                                                        acacCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCaa                                                                    <  2:264818/151‑1 (MQ=255)                                                                         cacCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCaa                                                                    >  2:863707/1‑150 (MQ=255)                                                                          acCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCaaa                                                                   >  1:30870/1‑150 (MQ=255)                                                                            cTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAAt                                                                 >  2:268092/1‑150 (MQ=255)                                                                             tACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGTGCGCGAGAAATCGAACGCCAAAATg                                                                >  1:654368/1‑150 (MQ=255)                                                                              aCGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATg                                                                <  2:730354/149‑1 (MQ=255)                                                                               cGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGAtt                                                             >  1:247540/1‑151 (MQ=255)                                                                                      aaTGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCg                                                      >  1:508970/1‑151 (MQ=255)                                                                                         gatgGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATc                                                   >  1:550729/1‑151 (MQ=255)                                                                                            ggCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCcg                                                 <  2:913256/150‑1 (MQ=255)                                                                                                 cAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGatat                                           <  1:476515/151‑1 (MQ=255)                                                                                                   gCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGAtata                                          <  1:438321/150‑1 (MQ=255)                                                                                                       cacGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGAt                                     <  2:433006/151‑1 (MQ=255)                                                                                                        acGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGAt                                     >  1:572094/1‑150 (MQ=255)                                                                                                            atatTGAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGc                                >  2:187530/1‑151 (MQ=255)                                                                                                                 gAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCg                           >  2:18625/1‑151 (MQ=255)                                                                                                                 gAAAACAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCg                           <  2:890300/151‑1 (MQ=255)                                                                                                                    aaCAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTc                        <  1:291795/151‑1 (MQ=255)                                                                                                                      cAGCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGt                      <  1:808281/151‑1 (MQ=255)                                                                                                                        gCAAGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGa                    <  1:522211/151‑1 (MQ=255)                                                                                                                          aaGCTGGTCGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGa                    >  2:296947/1‑149 (MQ=255)                                                                                                                                  cGTGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATgg          <  1:751232/151‑1 (MQ=255)                                                                                                                                    tGATGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCt        >  1:590212/1‑151 (MQ=255)                                                                                                                                      aTGTTTGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTcc      >  1:510950/1‑151 (MQ=255)                                                                                                                                          ttGGGAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTAt  >  1:194317/1‑151 (MQ=255)                                                                                                                                              |                                                                                                                                                   GAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATGATGGCAACAGCACCACGGATATTGAAAACAGCAAGCTGGTCGTGATGTTTGGTAACAACCCGGCAGAAACCCGCATGAGCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTAT  >  NC_000913/1658614‑1658901

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