Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,662,571	T→C	G672G (GGT→GGC)	ynfF →	putative selenate reductase YnfF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,662,571	0	T	C	100.0%	119.3 / NA	36	G672G (GGT→GGC)	ynfF	putative selenate reductase YnfF
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/21);  total (15/21)

TCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATC  >  NC_000913/1662439‑1662706                                                                                                                                     |                                                                                                                                        tCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATg                                                                                                                                  >  1:1052279/1‑140 (MQ=255) tCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATg                                                                                                                                  >  1:1338221/1‑140 (MQ=255)    acaGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCt                                                                                                                               <  1:730165/140‑1 (MQ=255)                       acacCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGaa                                                                                                            <  1:157200/140‑1 (MQ=255)                          ccTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAAc                                                                                                         <  1:1376481/140‑1 (MQ=255)                           cTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCgg                                                                                                              >  1:1317513/1‑134 (MQ=255)                           cTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCgg                                                                                                              <  2:1317513/134‑1 (MQ=255)                            ttCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCt                                                                                                       >  2:170646/1‑140 (MQ=255)                            ttCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCt                                                                                                       <  2:349145/140‑1 (MQ=255)                              cGGGGAAAATCGAATTTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATcccc                                                                                                                   <  1:210474/126‑1 (MQ=255)                              cGGGGAAAATCGAATTTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATcccc                                                                                                                   >  2:210474/1‑126 (MQ=255)                              cGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCt                                                                                                       <  2:972598/138‑1 (MQ=255)                              cGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCt                                                                                                       >  1:972598/1‑138 (MQ=255)                              cGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTAt                                                                                                     >  1:278472/1‑140 (MQ=255)                               ggggAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATc                                                                                                    <  1:744253/140‑1 (MQ=255)                                 ggAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCa                                                                                                  <  2:738366/140‑1 (MQ=255)                                    aaaTCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTg                                                                                               <  2:571937/140‑1 (MQ=255)                                                aTTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGc                                                                                   <  2:278472/140‑1 (MQ=255)                                                aTTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGc                                                                                   <  2:915292/140‑1 (MQ=255)                                                aTTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATCCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGc                                                                                   <  2:79941/140‑1 (MQ=255)                                                    tGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTcc                                                                               <  2:1338221/140‑1 (MQ=255)                                                       aCGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCAtt                                                                            >  1:975854/1‑140 (MQ=255)                                                           cTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACaa                                                                        <  2:997589/140‑1 (MQ=255)                                                           cTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACaa                                                                        <  1:1183945/140‑1 (MQ=255)                                                                    aTTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTAc                                                               >  2:28106/1‑140 (MQ=255)                                                                     ttGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTAcc                                                              <  2:1085290/140‑1 (MQ=255)                                                                     ttGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTAcc                                                              <  2:166012/140‑1 (MQ=255)                                                                         aGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACt                                                          <  2:107792/140‑1 (MQ=255)                                                                             aCCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATACCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAg                                                      >  2:405894/1‑140 (MQ=255)                                                                              ccTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGc                                                     >  1:1180342/1‑140 (MQ=255)                                                                                   gAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACgg                                                <  2:975854/140‑1 (MQ=255)                                                                                                              ccccTTCCTGCGTATACCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTg                     >  2:1061441/1‑140 (MQ=255)                                                                                                                          tataCCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGgt         >  1:206869/1‑140 (MQ=255)                                                                                                                             aCCCCAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGtgtg      >  1:133563/1‑140 (MQ=255)                                                                                                                                ccAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGAt   <  1:28106/140‑1 (MQ=255)                                                                                                                                 cAGGCTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATc  >  2:767523/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                                        TCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATC  >  NC_000913/1662439‑1662706

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: ‑