Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,008,594	G→A	P106P (CCG→CCA)	yedE →	putative selenium transporter YedE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,008,594	0	G	A	100.0%	159.8 / NA	46	P106P (CCG→CCA)	yedE	putative selenium transporter YedE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/22);  total (24/22)

GAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTT  >  NC_000913/2008457‑2008729                                                                                                                                          |                                                                                                                                        gAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCcacg                                                                                                                                       <  1:663871/140‑1 (MQ=255) gAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCcacg                                                                                                                                       >  1:425717/1‑140 (MQ=255)       tGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCg                                                                                                                                 >  1:465056/1‑140 (MQ=255)         gggTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTa                                                                                                                               >  1:572364/1‑140 (MQ=255)         gggTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTa                                                                                                                               >  1:635047/1‑140 (MQ=255)          ggTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTAt                                                                                                                              >  2:255158/1‑140 (MQ=255)          ggTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTAt                                                                                                                              <  1:186625/140‑1 (MQ=255)               cTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCa                                                                                                                         >  1:79797/1‑140 (MQ=255)               cTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCa                                                                                                                         <  1:589545/140‑1 (MQ=255)               cTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCa                                                                                                                         <  2:271869/140‑1 (MQ=255)                   aaaaTTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTAt                                                                                                                     >  1:587326/1‑140 (MQ=255)                    aaaTTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATg                                                                                                                    <  2:34836/140‑1 (MQ=255)                      aTTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCa                                                                                                                  >  1:262770/1‑140 (MQ=255)                          tCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGcc                                                                                                              >  1:615813/1‑140 (MQ=255)                          tCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGcc                                                                                                              >  1:264435/1‑140 (MQ=255)                          tCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCT‑GGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTAt                                                                                                                              <  1:342032/123‑1 (MQ=255)                          tCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCT‑GGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTAt                                                                                                                              >  2:342032/1‑123 (MQ=255)                                        tCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCAt                                                                                                >  2:102849/1‑140 (MQ=255)                                         cGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCAtt                                                                                               <  2:465056/140‑1 (MQ=255)                                         cGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCAtt                                                                                               >  1:337611/1‑140 (MQ=255)                                           ccATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTAGGCAGGACATCATTGGCGGCATTAt                                                                                             >  2:726897/1‑140 (MQ=255)                                                      cATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTg                                                                                  >  2:568875/1‑140 (MQ=255)                                                               gatgatgatCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTc                                                                         <  1:312117/140‑1 (MQ=255)                                                                    tgatCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCa                                                                    <  1:568875/140‑1 (MQ=255)                                                                     gatCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTc                                                                         <  1:26932/134‑1 (MQ=255)                                                                     gatCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTc                                                                         >  2:26932/1‑134 (MQ=255)                                                                      atCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAAt                                                                  <  2:79797/140‑1 (MQ=255)                                                                      atCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAAt                                                                  <  1:255158/140‑1 (MQ=255)                                                                       tCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATg                                                                 <  2:425717/140‑1 (MQ=255)                                                                          tCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATgg                                                                >  2:104815/1‑138 (MQ=255)                                                                          tCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATgg                                                                <  1:104815/138‑1 (MQ=255)                                                                             gTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGt                                                           >  2:180295/1‑140 (MQ=255)                                                                                 gTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAAcc                                                       >  1:510563/1‑140 (MQ=255)                                                                                        ggcTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTggcggc                                                >  1:8612/1‑140 (MQ=255)                                                                                         gctgctTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGgcggcg                                               <  2:635047/140‑1 (MQ=255)                                                                                            gctTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTc                                            >  2:665231/1‑140 (MQ=255)                                                                                                 gCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCTTTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTAc                                       <  1:726897/140‑1 (MQ=255)                                                                                                 gCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTAc                                       <  2:615813/140‑1 (MQ=255)                                                                                                           tgtgGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCt                             <  2:337611/140‑1 (MQ=255)                                                                                                            gtgGGCCAACAATGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGTGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTc                            <  2:264435/140‑1 (MQ=255)                                                                                                               ggCCAACAATGTCAAACTACGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGt                         <  1:123343/140‑1 (MQ=255)                                                                                                                  caacaaTGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTtct                      <  1:43130/140‑1 (MQ=255)                                                                                                                  caacaaTGTCAAACTGCGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTtct                      <  2:262770/140‑1 (MQ=255)                                                                                                                                tGCGAATGCCACGCAGGCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTgg        >  2:642321/1‑140 (MQ=255)                                                                                                                                  cGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGtt      >  1:666026/1‑140 (MQ=255)                                                                                                                                  cGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGtt      >  1:341638/1‑140 (MQ=255)                                                                                                                                  cGAATGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTCGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGtt      >  2:21324/1‑140 (MQ=255)                                                                                                                                      tGCCACGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCttt  <  2:666026/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                        GAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCAGTTCTCGCTGCATGCCTGGTTCTTT  >  NC_000913/2008457‑2008729

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