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%	156.2 / NA	38	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 . (18/20);  total (18/20)

GGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAACTT  >  NC_000913/1624220‑1624471                                                                                                                                   |                                                                                                                          ggAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCAtttttttgc                                                                                                                   <  2:1153068/139‑3 (MQ=255)     tAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGcatca                                                                                                              <  2:1175592/139‑1 (MQ=255)      aaCGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGcatcat                                                                                                             <  2:566160/139‑1 (MQ=255)         gACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCgg                                                                                                          >  1:206607/1‑139 (MQ=255)            aTTAATTTCACCGCTCTTATCGCGGAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTc                                                                                                       >  1:612381/1‑139 (MQ=255)                    cACCGCTCTTATCGCGCACAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCGATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAg                                                                                               <  1:976467/139‑1 (MQ=255)                     aCCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGc                                                                                              >  1:471513/1‑139 (MQ=255)                     aCCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGc                                                                                              <  2:1209909/139‑1 (MQ=255)                         ctctTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                                          <  2:206607/139‑1 (MQ=255)                         ctctTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                                          <  1:740577/139‑1 (MQ=255)                          tctTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtct                                                                                         >  2:660368/1‑139 (MQ=255)                                    aaaaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTa                                                                                     <  1:1282653/133‑1 (MQ=255)                                    aaaaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTa                                                                                     >  2:1282653/1‑133 (MQ=255)                                                          ttGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATc                                                         <  2:670451/139‑1 (MQ=255)                                                          ttGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATc                                                         <  2:471513/139‑1 (MQ=255)                                                           tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCGTCGGTTCCGCGATAGCCAtc                                                                                          <  1:722672/105‑1 (MQ=255)                                                           tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                                          >  2:722672/1‑105 (MQ=255)                                                           tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATcc                                                        >  1:603824/1‑139 (MQ=255)                                                             tgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGgcg                                                                >  2:1031806/1‑129 (MQ=255)                                                             tgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGgcg                                                                <  1:1031806/129‑1 (MQ=255)                                                                gaAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGttt                                                   >  2:541989/1‑139 (MQ=255)                                                                 aaaTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGAt                                                                                                 <  1:749740/92‑1 (MQ=255)                                                                 aaaTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGAt                                                                                                 >  2:749740/1‑92 (MQ=255)                                                                  aaTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCa                                                           >  1:734537/1‑129 (MQ=255)                                                                  aaTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCa                                                           <  2:734537/129‑1 (MQ=255)                                                                    tCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTaaa                                               <  1:660368/139‑1 (MQ=255)                                                                                            aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCa                                                                                >  2:937346/1‑82 (MQ=255)                                                                                            aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCa                                                                                <  1:937346/82‑1 (MQ=255)                                                                                            aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCa                       <  2:603824/139‑1 (MQ=255)                                                                                             aGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCAt                      >  2:1326265/1‑139 (MQ=255)                                                                                                    tCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCAt               >  2:712001/1‑139 (MQ=255)                                                                                                         cATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACaa          >  1:587667/1‑139 (MQ=255)                                                                                                              tttttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCAt               <  2:439000/129‑1 (MQ=255)                                                                                                              tttttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCAt               >  1:439000/1‑129 (MQ=255)                                                                                                              tttttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGaaaa     <  1:541989/139‑1 (MQ=255)                                                                                                              tttttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGaaaa     <  1:712001/139‑1 (MQ=255)                                                                                                                 ttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGtt                                                                                                        >  1:1293224/1‑37 (MQ=37)                                                                                                                 ttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGtt                                                                                                        <  2:1293224/37‑1 (MQ=37)                                                                                                                 ttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAACtt  >  2:543549/1‑139 (MQ=255)                                                                                                                                   |                                                                                                                          GGAATAACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAACTT  >  NC_000913/1624220‑1624471

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