Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,660,686	G→T	R44L (CGG→CTG)	ynfF →	putative selenate reductase YnfF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,660,686	0	G	T	100.0%	121.6 / NA	35	R44L (CGG→CTG)	ynfF	putative selenate reductase YnfF
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (16/19);  total (16/19)

AGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACACCGTGTACTGGGTGGAGTCTGATACGACAGGTG  >  NC_000913/1660549‑1660817                                                                                                                                          |                                                                                                                                    aGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCtgg                                                                                                                                   >  2:395332/1‑140 (MQ=255)  gTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTggc                                                                                                                                  >  1:155395/1‑140 (MQ=255)      tgatgaAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCa                                                                                                                              <  2:86860/140‑1 (MQ=255)             aaTCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTc                                                                                                                       >  1:310906/1‑140 (MQ=255)             aaTCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTc                                                                                                                       <  1:395332/140‑1 (MQ=255)                  aTACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTa                                                                                                                  <  2:310906/140‑1 (MQ=255)                         agagGCGCTCATGAAGGCTGAAATCAGCCGCAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAgaga                                                                                                           >  1:72775/1‑140 (MQ=255)                            ggCGCTGATGAAGGCTGAAATCAGCAGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCGCTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAg                                                                                                        <  2:72775/140‑1 (MQ=255)                                cTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGt                                                                                                    >  1:19483/1‑140 (MQ=255)                                cTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGt                                                                                                    >  2:225960/1‑140 (MQ=255)                                         gCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTc                                                                                           >  2:329977/1‑140 (MQ=255)                                                       tAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACt                                                                             >  2:273322/1‑140 (MQ=255)                                                                  tgaAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGc                                                                  <  2:19483/140‑1 (MQ=255)                                                                  tgaAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGc                                                                  <  1:225960/140‑1 (MQ=255)                                                                   gaAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGct                                                                 <  1:308699/140‑1 (MQ=255)                                                                           ccGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGTGCCGCTGTCTGTTa                                                         >  1:338557/1‑140 (MQ=255)                                                                           ccGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTa                                                         >  1:97622/1‑140 (MQ=255)                                                                                 ttGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTg                                                   >  1:129641/1‑140 (MQ=255)                                                                                   ggCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCa                                                 >  2:152485/1‑140 (MQ=255)                                                                                          ttGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGaaa                                          >  1:365612/1‑140 (MQ=255)                                                                                              gcTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGTGCCGCTGTCTGTTACGTTTGCATGTGAAAGATg                                      <  2:338557/140‑1 (MQ=255)                                                                                              gcTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATg                                      <  2:280818/140‑1 (MQ=255)                                                                                              gcTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATg                                      <  1:329977/140‑1 (MQ=255)                                                                                                  gCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATg                                      <  1:401493/136‑1 (MQ=255)                                                                                                  gCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATg                                      >  2:401493/1‑136 (MQ=255)                                                                                                   cAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACAcc                                 <  2:35634/140‑1 (MQ=255)                                                                                                      gCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACACCgtg                              <  2:68511/140‑1 (MQ=255)                                                                                                         gTGCTATCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACACCGTGTAc                           <  2:171636/140‑1 (MQ=255)                                                                                                               tCACTCTGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACACCGTGTACTGGGTg                     <  2:155395/140‑1 (MQ=255)                                                                                                                  ctctGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACACCGTGTACTGGGTGGAg                  <  1:152485/140‑1 (MQ=255)                                                                                                                    ctGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACACCGTGTAc                           >  2:388858/1‑129 (MQ=255)                                                                                                                    ctGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACACCGTGTAc                           <  1:388858/129‑1 (MQ=255)                                                                                                                     tGCCATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACACCGTGTACTGGGTGGAGTCt               <  2:365612/140‑1 (MQ=255)                                                                                                                        cATTTTCCCAAATGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACACCGTGTACTGGGTGGAGTCTGAt            >  2:273476/1‑140 (MQ=255)                                                                                                                                  aaTGGTCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACACCGTGTACTGGGTGGAGTCTGATACGACAGGTg  >  2:365007/1‑140 (MQ=255)                                                                                                                                       tCCTGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTAACGTTTGCATGTGAAAGATGACACCGTGTACTGGGTGTAGTCTGATACGa        <  2:324831/129‑1 (MQ=255)                                                                                                                                          |                                                                                                                                    AGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTGCGGGAGCCGCTGTCTGTTACGTTTGCATGTGAAAGATGACACCGTGTACTGGGTGGAGTCTGATACGACAGGTG  >  NC_000913/1660549‑1660817

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