Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,659,873	A→G	Y602C (TAC→TGC)	ynfE →	putative selenate reductase YnfE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,659,873	0	A	G	100.0%	118.3 / NA	34	Y602C (TAC→TGC)	ynfE	putative selenate reductase YnfE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/19);  total (15/19)

GGATCCTGAGTGAAGTCGCGAAACGTCTTGGACCTGACGTCTATCAAAAGTTCACAGAAGGTCGCACGCAGGAACAATGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTACGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCT  >  NC_000913/1659741‑1659982                                                                                                                                     |                                                                                                              ggATCCTGAGTGAAGTCGCGAAACGTCTTGGACCTGACGTCTATCAAAAGTTCACAGAAGGTCGCACGCAGGAACAATGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGaa                                                                                                        <  2:6238/140‑1 (MQ=255)       tGAGTGAAGTCGCGAAACGTCTTGGACCTGACGTCTATCAAAAGTTCACAGAAGGTCGCACGCAGGAACAATGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGaaa                                                                                                  >  2:933816/1‑140 (MQ=255)          gTGAAGTCGCGAAACGTCTTGGACCTGACGTCTATCAAAAGTTCACAGAAGGTCGCACGCAGGAACAATGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGaaaaaa                                                                                               <  1:542994/140‑1 (MQ=255)            gAAGTCGCGAAACGTCTTGGACCTGACGTCTATCAAAAGTTCACAGAAGGTCGCACGCAGGAACAATGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAAt                                                                                             <  2:44747/140‑1 (MQ=255)            gAAGTCGCGAAACGTCTTGGACCTGACGTCTATCAAAAGTTCACAGAAGGTCGCACGCAGGAACAATGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAAt                                                                                             <  2:732537/140‑1 (MQ=255)             aaGTCGCGAAACGTCTTGGACCTGACGTCTATCAAAAGTTCACAGAAGGTCGCACGCAGGAACAATGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATg                                                                                            >  1:908280/1‑140 (MQ=255)              aGTCGCGAAACGTCTTGGACCTGACGTCTATCAAAAGTTCACAGAAGGTCGCACGCAGGAACAATGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATgg                                                                                           <  2:908280/140‑1 (MQ=255)              aGTCGCGAAACGTCTTGGACCTGACGTCTATCAAAAGTTCACAGAAGGTCGCACGCAGGAACAATGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATgg                                                                                           <  1:57251/140‑1 (MQ=255)              aGTCGCGAAACGTCTTGGACCTGACGTCTATCAAAAGTTCACAGAAGGTCGCACGCAGGAACAATGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATgg                                                                                           >  2:962078/1‑140 (MQ=255)              aGTCGCGAAACGTCTTGGACCTGACGTCTATCAAAAGTTCACAGAAGGTCGCACGCAGGAACAATGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATgg                                                                                           <  2:403086/140‑1 (MQ=255)               gTCGCGAAACGTCTTGGACCTGACGTCTATCAAAAGTTCACAGAAGGTCGCGCGCAGGAACAATGATTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATggg                                                                                          >  2:852667/1‑140 (MQ=255)                     aaaCGTCTTGGACCTGACGTCTATCAAAAGTTCACAGAAGGTCGCACGCAGGAACAATGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCta                                                                                    <  1:962078/140‑1 (MQ=255)                       aCGTCTTGGACCTGACGTCTATCAAAAGTTCACAGAAGGTCGCACGCAGGAACAATGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCtata                                                                                  <  1:623459/140‑1 (MQ=255)                                 ccTGACGTCTATCAAAAGTTCACAGAAGGTCGCACGCAGGAACAATGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGa                                                                        <  2:592316/140‑1 (MQ=255)                                  ctcaCGTCTATCAAAAGTTCACAGAAAGTCGCACGCAGGAACAATGGTTACAACATCTGTACGCCAAAAAGCTTGCCAAAGATCCGGCGATACCGTCTTGCGACGAACAGAAAAAAATGGGTATCTATAAGCGTAAAGAt                                                                       >  2:226964/4‑140 (MQ=255)                                        tCTATCAAAAGTTCACAGAAGGTCGCACGCAGGAACAATGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAAt                                                                 <  2:656928/140‑1 (MQ=255)                                                      acaGAAGGTCGCACGCAGGAACAATGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGc                                                   <  1:356707/140‑1 (MQ=255)                                                               cgcacgcaGGAACAATGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGc                                          >  1:604319/1‑140 (MQ=255)                                                                gcacgcaGGAACAATGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCa                                         >  2:690084/1‑140 (MQ=255)                                                                              tGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTgac                                >  1:580096/1‑135 (MQ=255)                                                                              tGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTgac                                <  2:580096/135‑1 (MQ=255)                                                                              tGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGAccc                           <  1:690084/140‑1 (MQ=255)                                                                              tGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGAccc                           <  2:608191/140‑1 (MQ=255)                                                                                gTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCg                         >  2:206461/1‑140 (MQ=255)                                                                                gTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCg                         >  1:549117/1‑140 (MQ=255)                                                                                 ttACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGa                        <  2:604319/140‑1 (MQ=255)                                                                                         tCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTGCCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCaa                   >  2:915131/1‑137 (MQ=255)                                                                                         tCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTGCCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCaa                   <  1:915131/137‑1 (MQ=255)                                                                                         tCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATc                >  1:491524/1‑140 (MQ=255)                                                                                             tACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACt            >  2:181254/1‑140 (MQ=255)                                                                                               cGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTa          >  2:278068/1‑140 (MQ=255)                                                                                                  cAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCATAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTaaaa       >  2:965647/1‑140 (MQ=255)                                                                                                   aaaaTGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAAc      >  2:412102/1‑140 (MQ=255)                                                                                                      aTGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGcc   <  2:426389/140‑1 (MQ=255)                                                                                                       tGCTTGCCAAAGATCCGGCGTTACCGTCTTGCGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCt  <  1:408002/140‑1 (MQ=255)                                                                                                                                     |                                                                                                              GGATCCTGAGTGAAGTCGCGAAACGTCTTGGACCTGACGTCTATCAAAAGTTCACAGAAGGTCGCACGCAGGAACAATGGTTACAACATCTGTACGCCAAAATGCTTGCCAAAGATCCGGCGTTACCGTCTTACGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCT  >  NC_000913/1659741‑1659982

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