Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,660,102	C→T	N678N (AAC→AAT)	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,102	0	C	T	100.0%	91.4 / NA	29	N678N (AAC→AAT)	ynfE	putative selenate reductase YnfE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/13);  total (16/13)

ACCGTCTTACGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTT  >  NC_000913/1659865‑1660339                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               aCCGTCTTACGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGt                                                                                                                                                                                                                                    >  1:109383/1‑250 (MQ=255)            acgaACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCa                                                                                                                                                                                                                        <  1:205549/251‑1 (MQ=255)                           ggggTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGt                                                                                                                                                                                                        <  1:139886/250‑1 (MQ=255)                                    ataAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTGAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCGTCCCACTGCAAC‑TGTTTGGTTTCCAtt                                                                                                                                                                                                 <  1:361855/249‑1 (MQ=255)                                            aaaGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCg                                                                                                                                                                                       >  1:290914/1‑251 (MQ=255)                                                   ccAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCGTactcac                                                                                                                                                                                <  1:504098/251‑1 (MQ=255)                                                            aTTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGTGGCAACTCCACTTAACACGCCTTCCGTTAAGATTGAATTTTATTCCCGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGAGGAAGTGATAGGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACT‑CAACCTGTTTGGTTTCCATTCCAAATCCCGTACTCCCTCgccccac                                                                                                                                                                       >  2:33870/1‑245 (MQ=255)                                                                    gCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAatat                                                                                                                                                               >  1:442947/1‑251 (MQ=255)                                                                     gcTACAAAGCGTTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGTTGAAGTGATAAGCCCCTTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCAACTGCAAC‑TGATTGGTTTCCTTTAAAAATACCGTACTCAATCGGCCTACGGCAATAtt                                                                                                                                                              >  2:104217/2‑251 (MQ=255)                                                                        aCAAAGCATGTCGTTACGACCCCGAGGCAAATCCCCTTAAAACGCCTTCCGGTAAGATGGAAATTTATTCCAGCAGGCTGGCGGAAAGTGCCCGTACCTGGGAACTTGAAAAAGATGAAGTGATAAGCCCACTGCCAGTTTATGCCTCACCCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGTTTTCCATATCAAATCCCGTACTCCCTCGCCCTACGGCAATACTGAt                                                                                                                                                           <  2:299713/251‑1 (MQ=255)                                                                             gCATTTCGAGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGATGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGAGATAAGCCCATTGCCAGTTTAAGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCt                                                                                                                                                      >  1:499304/1‑251 (MQ=255)                                                                                tttCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGa                                                                                                                                                    >  1:65626/1‑250 (MQ=255)                                                                                       acgacCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGGTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGctgcc                                                                                                                                            <  2:191469/251‑1 (MQ=255)                                                                                           ccccGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTGTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTGCAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCtgcctgcc                                                                                                                                        <  1:473537/251‑1 (MQ=255)                                                                                                    aaaTCCCCTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGAGGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCaggagg                                                                                                                               <  2:329546/251‑1 (MQ=255)                                                                                                                      gtCCGGGAAGATTGAAATTTATTACAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAAATGGTAAAAGATGACGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACGGCAAC‑TGTTTGGTTTCCCTTACAAATCCCGTACTCACTCGACCTCCGGCAATATTTATCTCCTGCAGGCTGCCTACCGTCAGGAGGTGTGGAACAACCCtat                                                                                                               >  2:315831/2‑249 (MQ=255)                                                                                                                       tCCGGTGAGAGTGAAATTTATTCCAGCAGGCTGGCGGGAATTGCCCGTACCTGGGAACTAGAAAAAGATGAAGTGATAAGCCCATGGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCACCTGAACGAAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTTCAAATCCCGTACTCACTCGCCCTACGGCAATATTGGTCTCATGAAGGCTGCCTGCCGTCAGACGGTGGGGATCAACCCTATAgg                                                                                                            >  2:278414/1‑250 (MQ=255)                                                                                                                       tCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGCTCCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGATCCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAg                                                                                                             <  2:388056/250‑1 (MQ=255)                                                                                                                                    aaaTTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGt                                                                                               <  2:250370/251‑1 (MQ=255)                                                                                                                                        ttATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAAc                                                                                                                    >  2:269530/1‑226 (MQ=255)                                                                                                                                        ttATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAAc                                                                                                                    <  1:269530/226‑1 (MQ=255)                                                                                                                                               agcagGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGGAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCa                                                                                     >  1:456339/1‑250 (MQ=255)                                                                                                                                                             aaTTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGc                                                                      <  1:227909/251‑1 (MQ=255)                                                                                                                                                                   ccGTACCTGGGAACTTGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGACCGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGAGCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGgtgt                                                                <  2:265804/251‑1 (MQ=255)                                                                                                                                                                             gAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCg                                                      >  2:457935/1‑251 (MQ=255)                                                                                                                                                                                    aaaaaGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGa                                                >  1:200930/1‑250 (MQ=255)                                                                                                                                                                                           tGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGc                                        >  1:367192/1‑251 (MQ=255)                                                                                                                                                                                                                 gTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGCGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAacacaac                   >  1:460548/1‑247 (MQ=255)                                                                                                                                                                                                                                 ttGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGCGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTCCCAGCCAAAGTAACACCACGTATTCTCCCTGGAGtt  >  2:73931/1‑251 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               ACCGTCTTACGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAAC‑TGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTT  >  NC_000913/1659865‑1660339

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