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%	203.3 / NA	49	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 . (23/26);  total (23/26)

AACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAACTTGCTAATA  >  NC_000913/1624225‑1624478                                                                                                                              |                                                                                                                                 aaCGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGcatcat                                                                                                                    <  2:972366/139‑1 (MQ=255)     aCATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGt                                                                                                                <  1:1423039/139‑1 (MQ=255)      cATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGtt                                                                                                               <  2:372423/139‑1 (MQ=255)        ttAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTcc                                                                                                             >  1:1085913/1‑139 (MQ=255)            tttCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGa                                                                                                         <  1:569145/139‑1 (MQ=255)                    ctctTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGa                                                                                                         <  1:75221/131‑1 (MQ=255)                    ctctTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGa                                                                                                         >  2:75221/1‑131 (MQ=255)                    ctctTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                                                 >  2:94845/1‑139 (MQ=255)                         aTCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTa                                                                                            >  2:1527057/1‑139 (MQ=255)                               aaaaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCAt                                                                                      >  2:1451133/1‑139 (MQ=255)                                 aaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                                                 <  1:449081/126‑1 (MQ=255)                                 aaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                                                 >  2:449081/1‑126 (MQ=255)                                           gTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTg                                                                          <  1:1296879/139‑1 (MQ=255)                                                 tcCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                                                 >  2:232031/1‑110 (MQ=255)                                                 tcCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                                                 <  1:232031/110‑1 (MQ=255)                                                     ttGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATc                                                                <  2:1377921/139‑1 (MQ=255)                                                     ttGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATc                                                                <  1:1527057/139‑1 (MQ=255)                                                      tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCcg                                                                                                            <  2:1232553/94‑1 (MQ=255)                                                      tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCcg                                                                                                            >  1:1232553/1‑94 (MQ=255)                                                      tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTg                                                                          >  1:1485626/1‑128 (MQ=255)                                                      tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTg                                                                          <  2:1485626/128‑1 (MQ=255)                                                      tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATcc                                                               >  2:881659/1‑139 (MQ=255)                                                                    ttACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGcc                                                 <  2:1085913/139‑1 (MQ=255)                                                                        cgcgTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGcatcat                                                                                                                    >  2:803545/1‑68 (MQ=255)                                                                        cgcgTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGcatcat                                                                                                                    <  1:803545/68‑1 (MQ=255)                                                                         gcgTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGt                                            >  1:922845/1‑139 (MQ=255)                                                                                 aaaCCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGc                                     <  2:242172/138‑1 (MQ=255)                                                                                 aaaCCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGc                                     >  1:242172/1‑138 (MQ=255)                                                                                  aaCCGAAGCGCCATCGGGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCa                                   <  1:929246/139‑1 (MQ=255)                                                                                    ccGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAAc                                 >  1:1107803/1‑139 (MQ=255)                                                                                       aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGcc                                                                                                    <  2:112717/69‑1 (MQ=255)                                                                                       aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGcc                                                                                                    >  1:112717/1‑69 (MQ=255)                                                                                       aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCa                              <  1:94845/139‑1 (MQ=255)                                                                                        aGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCAtttttttgc                                                                                                                          <  2:1450094/47‑3 (MQ=255)                                                                                        aGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCAtttttttgc                                                                                                                          >  1:1450094/1‑45 (MQ=255)                                                                                            ccATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGa                                                                                                         <  1:1026854/59‑1 (MQ=255)                                                                                            acATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTACATCATCGGTTCCGCGa                                                                                                         >  2:1026854/2‑59 (MQ=38)                                                                                               tCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCAt                      <  1:881659/139‑1 (MQ=255)                                                                                               tCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCAt                      >  2:288314/1‑139 (MQ=255)                                                                                                  tGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtctt                                                                                               >  1:718294/1‑63 (MQ=255)                                                                                                  tGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtctt                                                                                               <  2:718294/63‑1 (MQ=255)                                                                                                    cATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACaa                 <  2:922845/139‑1 (MQ=255)                                                                                                         tttttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAAc                                                     >  2:291184/1‑98 (MQ=255)                                                                                                         tttttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAAc                                                     <  1:291184/98‑1 (MQ=255)                                                                                                         tttttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGaaaa            <  2:1356714/139‑1 (MQ=255)                                                                                                            ttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAACtt         >  2:253380/1‑139 (MQ=255)                                                                                                                   aTCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAg                                                                                                      <  2:1568157/39‑1 (MQ=37)                                                                                                                   aTCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAg                                                                                                      >  1:1568157/1‑39 (MQ=37)                                                                                                                   aTCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGcc                                                                                                    >  1:1440554/1‑41 (MQ=38)                                                                                                                   aTCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGcc                                                                                                    <  2:1440554/41‑1 (MQ=38)                                                                                                                   aTCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAACTTGCTAATa  >  1:23994/1‑139 (MQ=255)                                                                                                                              |                                                                                                                                 AACGACATTAATTTCACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAACTTGCTAATA  >  NC_000913/1624225‑1624478

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