Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,629,247	C→T	I11I (ATC→ATT)	ydfZ →	putative selenoprotein YdfZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,629,247	0	C	T	100.0%	55.9 / NA	16	I11I (ATC→ATT)	ydfZ	putative selenoprotein YdfZ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/8);  total (8/8)

TCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATCACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTT  >  NC_000913/1629001‑1629471                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                 tCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATtacc                                                                                                                                                                                                                               >  1:118912/1‑250 (MQ=255)    gCCAGGAAAACAGCGACTGGTTTAGCGCAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTg                                                                                                                                                                                                                           >  2:79883/1‑251 (MQ=255)     ccAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCCTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTgg                                                                                                                                                                                                                          <  1:64865/251‑1 (MQ=255)                          tAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGCTCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCgg                                                                                                                                                                                                     >  2:185952/1‑251 (MQ=255)                                     aaTTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCGTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTcac                                                                                                                                                                                           <  2:152078/250‑1 (MQ=255)                                                 ctcATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGAt                                                                                                                                                                                <  2:140647/249‑1 (MQ=255)                                                        catccGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGAt                                                                                                                                                                       >  1:185233/1‑251 (MQ=255)                                                                          tttCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGa                                                                                                                                                      <  2:237739/250‑1 (MQ=255)                                                                                                gcatgaatCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTGATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCTCACGGGCATGATCCTGTCGATTGATACTGAATGTCTGTCCTCTGAGCAAATCCGCCGCgg                                                                                                                                <  2:139771/250‑1 (MQ=255)                                                                                                           aaaaGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAgt                                                                                                                    >  1:208998/1‑251 (MQ=255)                                                                                                                                       aaaTATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACTCCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTc                                                                                                                                                                          <  2:210688/169‑1 (MQ=255)                                                                                                                                       aaaTATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTc                                                                                                                                                                          >  1:210688/1‑169 (MQ=255)                                                                                                                                              tCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTg                                                                                  <  1:217264/250‑1 (MQ=255)                                                                                                                                                              tttATATTGCGACCGTTACGTGGGATTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCgg                                                                 >  1:206476/1‑251 (MQ=255)                                                                                                                                                                            gTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTCGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGAATCGTCTCGGCATGAACTCAGCgt                                                   >  2:129602/1‑251 (MQ=255)                                                                                                                                                                                                         aGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCttt                      >  2:135077/1‑251 (MQ=255)                                                                                                                                                                                                                             acCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATAtttt  <  1:9655/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                 TCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATCACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTT  >  NC_000913/1629001‑1629471

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