Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,660,030	T→C	R654R (CGT→CGC)	ynfE →	putative selenate reductase YnfE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,660,030	0	T	C	100.0%	128.2 / NA	35	R654R (CGT→CGC)	ynfE	putative selenate reductase YnfE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/15);  total (20/15)

TAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCAC  >  NC_000913/1659907‑1660165                                                                                                                            |                                                                                                                                        tAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGaaa                                                                                                                         >  1:544867/1‑140 (MQ=255) tAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGaaa                                                                                                                         >  1:75311/1‑140 (MQ=255) tAAAGATCCCAATGGCCATTTTGTCG‑CTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGaaaa                                                                                                                        >  2:146133/1‑140 (MQ=255)   aaGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCGAATCCACTTAAAACGCCTTCCCGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGaaaaa                                                                                                                       >  2:319514/1‑140 (MQ=255)   aaGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGaaaaa                                                                                                                       >  2:501867/1‑140 (MQ=255)   aaGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGaaaaa                                                                                                                       >  1:627502/1‑140 (MQ=255)   aaGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGaaaaa                                                                                                                       >  1:708689/1‑140 (MQ=255)         ccAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGaa                                                                                                                 >  2:761084/1‑140 (MQ=255)             tGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGa                                                                                                             <  2:218511/140‑1 (MQ=255)                        tCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTg                                                                                                  >  1:504399/1‑140 (MQ=255)                         cGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGc                                                                                                 <  1:469598/140‑1 (MQ=255)                            cTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAg                                                                                              >  1:1085562/1‑140 (MQ=255)                                             acgacCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCttt                                                                             >  1:51060/1‑140 (MQ=255)                                              cgacCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTg                                                                            <  1:1067117/140‑1 (MQ=255)                                                 ccccGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAg                                                                         >  2:573759/1‑140 (MQ=255)                                                    cGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTg                                                                            >  1:68880/1‑134 (MQ=255)                                                    cGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTg                                                                            <  2:68880/134‑1 (MQ=255)                                                          aaaTCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACt                                                                <  2:400253/140‑1 (MQ=255)                                                            aTCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTcc                                                              >  1:497789/1‑140 (MQ=255)                                                               cACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCt                                                           <  1:573759/140‑1 (MQ=255)                                                                     aaaCGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGt                                                     <  1:1055709/140‑1 (MQ=255)                                                                         gCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGaa                                                 <  2:75311/140‑1 (MQ=255)                                                                            ttCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCt                                              >  1:164886/1‑140 (MQ=255)                                                                            ttCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCt                                              >  2:575796/1‑140 (MQ=255)                                                                             tCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCtt                                             <  2:497789/140‑1 (MQ=255)                                                                                  tAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCAc                                        <  2:51060/140‑1 (MQ=255)                                                                                       ttGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCaa                                   <  2:627502/140‑1 (MQ=255)                                                                                         gAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCaa                                   <  1:845081/138‑1 (MQ=255)                                                                                         gAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCaa                                   >  2:845081/1‑138 (MQ=255)                                                                                                       cagGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCAtt                   >  1:871653/1‑140 (MQ=255)                                                                                                       cagGCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCAtt                   >  1:544314/1‑140 (MQ=255)                                                                                                         ggCTGGCGGAAATTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTAc                 <  2:1085562/140‑1 (MQ=255)                                                                                                                  aaaTTGCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGt        <  2:708689/140‑1 (MQ=255)                                                                                                                       gCCCGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTactca   <  2:871653/140‑1 (MQ=255)                                                                                                                        cccGCACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTactcac  >  1:73852/1‑140 (MQ=255)                                                                                                                            |                                                                                                                                        TAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCAC  >  NC_000913/1659907‑1660165

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