Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,624,357	(T)8→7	coding (141/393 nt)	ydeI ←	hydrogen peroxide resistance OB fold protein, putative periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,624,350	0	T	.	100.0%	143.3 / NA	34	coding (148/393 nt)	ydeI	hydrogen peroxide resistance OB fold protein, putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (14/20);  total (14/20)

CCGCAGGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAA  >  NC_000913/1624215‑1624468                                                                                                                                        |                                                                                                                       ccGCAGGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCAtttt                                                                                                                     >  1:279831/1‑139 (MQ=255)    cAGGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCAttttttt                                                                                                                  <  1:376227/139‑1 (MQ=255)                          aCCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGc                                                                                           >  2:58847/1‑139 (MQ=255)                             gCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAt                                                                                        <  1:227496/139‑1 (MQ=255)                                 ttATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAg                                                                                            >  1:675255/1‑131 (MQ=255)                                 ttATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAg                                                                                            <  2:675255/131‑1 (MQ=255)                                   aTCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtct                                                                                      >  1:204950/1‑135 (MQ=255)                                   aTCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtct                                                                                      <  2:204950/135‑1 (MQ=255)                                    tCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTAt                                                                                 >  1:591843/1‑139 (MQ=255)                                             aCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAgg                                                                        <  2:623489/139‑1 (MQ=255)                                              cGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAggc                                                                       <  2:349484/139‑1 (MQ=255)                                                 aaCGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAggcggc                                                                    >  2:656125/1‑139 (MQ=255)                                                      tCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTgg                                                               <  2:279831/139‑1 (MQ=255)                                                       ccTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGgc                                                              <  1:348152/139‑1 (MQ=255)                                                             cTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCa                                                        <  1:292783/139‑1 (MQ=255)                                                             cTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCa                                                        <  2:219461/139‑1 (MQ=255)                                                                tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATcc                                                     >  2:311097/1‑139 (MQ=255)                                                                  tgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                                       >  2:519761/1‑103 (MQ=255)                                                                  tgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                                       <  1:519761/103‑1 (MQ=255)                                                                        aTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCcgcg                                                                                                <  2:198808/88‑1 (MQ=255)                                                                        aTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCcgcg                                                                                                >  1:198808/1‑88 (MQ=255)                                                                            aaTTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCg                                         >  2:447759/1‑139 (MQ=255)                                                                               tACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGccc                                      <  1:114379/139‑1 (MQ=255)                                                                                 ccGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTg                                    >  2:74648/1‑139 (MQ=255)                                                                                                gAAGCGCCATCGTGCATGGTTTTGGCGAAATGAAAGGTCA‑TTTAATCGCATCGTCGGTTCCGGGATAGCCATCTCCTCTTGACTGAGGCGCCGGTGTCCCTGCATCCTGTTTGAAACCGCCCTGGTAATATGCCACCGc                     >  1:15753/1‑139 (MQ=255)                                                                                                 aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCa                    >  1:799113/1‑139 (MQ=255)                                                                                                 aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCa                    <  1:656125/139‑1 (MQ=255)                                                                                                 aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCa                    <  1:447759/139‑1 (MQ=255)                                                                                                      ccATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAgg               <  2:22466/139‑1 (MQ=255)                                                                                                                 ggTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCAt            <  2:280421/131‑1 (MQ=255)                                                                                                                 ggTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCAt            >  1:280421/1‑131 (MQ=255)                                                                                                                   tttttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGaaaa  <  1:74648/139‑1 (MQ=255)                                                                                                                   tttttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGaaaa  <  1:58847/139‑1 (MQ=255)                                                                                                                   tttttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGaaaa  <  2:799113/139‑1 (MQ=255)                                                                                                                     tttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGc                                        <  2:403818/99‑1 (MQ=255)                                                                                                                     tttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGc                                        >  1:403818/1‑99 (MQ=255)                                                                                                                                        |                                                                                                                       CCGCAGGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAA  >  NC_000913/1624215‑1624468

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