Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,497,308	A→G	E151E (GAA→GAG)	opgD →	OPG biosynthetic periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,497,308	0	A	G	100.0%	137.2 / NA	39	E151E (GAA→GAG)	opgD	OPG biosynthetic periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/18);  total (21/18)

TGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTT  >  NC_000913/1497178‑1497436                                                                                                                                   |                                                                                                                                 tgggaatgggaTTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAg                                                                                                                          >  1:805246/1‑139 (MQ=255) tgggaatgggaTTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAg                                                                                                                          >  2:188867/1‑139 (MQ=255)               cGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTTGCCCGGAGTTGTTCAAATACAACGCTGCAGGTGTTGATACAATACAATTAGAGGGGCAAAGCGATCTCGGCTTTg                                                                                                            >  1:361431/1‑139 (MQ=255)                    ccgcGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGt                                                                                                       >  1:599194/1‑139 (MQ=255)                         ttCGTATGTTTTCTGTAGATCGAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCATATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTGTTcg                                                                                                  <  2:783468/139‑1 (MQ=255)                                           aTCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGa                                                                                >  2:736841/1‑139 (MQ=255)                                               aGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTg                                                                            >  1:772971/1‑139 (MQ=255)                                               aGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTg                                                                            >  1:497299/1‑139 (MQ=255)                                                  aaCACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGgcg                                                                         <  2:497299/139‑1 (MQ=255)                                                       aTCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTaa                                                                                         <  1:675045/118‑1 (MQ=255)                                                       aTCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTaa                                                                                         >  2:675045/1‑118 (MQ=255)                                                       aTCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCg                                                                    <  2:805246/139‑1 (MQ=255)                                                               cgTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATgtag                                                            >  2:569119/1‑139 (MQ=255)                                                                 tGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGTGCGCCGTGATgtagta                                                          >  1:891506/1‑139 (MQ=255)                                                                 tGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATgtagta                                                          >  1:876515/1‑139 (MQ=255)                                                                        cACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGTGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGc                                                                                      >  1:736372/1‑104 (MQ=255)                                                                        cACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGc                                                                                      <  2:736372/104‑1 (MQ=255)                                                                                             aaaTACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCcgcg                              <  2:599194/139‑1 (MQ=255)                                                                                               aTACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGcc                            <  2:772971/139‑1 (MQ=255)                                                                                                tACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCg                           <  1:569119/139‑1 (MQ=255)                                                                                                 aCAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAGAGCCCCCGAACTGGCGCGCGGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGt                          <  2:361431/139‑1 (MQ=255)                                                                                                    aCGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGc                                                                                      <  1:132148/76‑1 (MQ=255)                                                                                                    aCGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGc                                                                                      >  2:132148/1‑76 (MQ=255)                                                                                                    aCGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGTGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTtga                       <  2:891506/139‑1 (MQ=255)                                                                                                         gCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGTCCCCGAACTGGCGCGCCg                                                                    >  1:745746/1‑89 (MQ=255)                                                                                                         gCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGTCCCCGAACTGGCGCGCCg                                                                    <  2:745746/89‑1 (MQ=255)                                                                                                         gCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATGTCCGCGCCGTTGATGATa                  <  2:747468/139‑1 (MQ=255)                                                                                                             gtgtTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGcc                                                                     <  1:861119/84‑1 (MQ=255)                                                                                                             gtgtTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGcc                                                                     >  2:861119/1‑84 (MQ=255)                                                                                                                ttGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCg                                                                    <  2:155630/82‑1 (MQ=255)                                                                                                                ttGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCg                                                                    >  1:155630/1‑82 (MQ=255)                                                                                                                 tGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGgcg                                                                         <  1:1213017/76‑1 (MQ=255)                                                                                                                 tGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGgcg                                                                         >  2:1213017/1‑76 (MQ=255)                                                                                                                  gATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAAt         >  1:865185/1‑139 (MQ=255)                                                                                                                       aaaaCAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGc                                                                                      >  1:850038/1‑57 (MQ=255)                                                                                                                       aaaaCAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGc                                                                                      <  2:850038/57‑1 (MQ=255)                                                                                                                        aaaCAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGgcg                                             >  1:962509/1‑97 (MQ=255)                                                                                                                        aaaCAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGgcg                                             <  2:962509/97‑1 (MQ=255)                                                                                                                         aaCAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGttt  >  2:854969/1‑139 (MQ=255)                                                                                                                                   |                                                                                                                                 TGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTT  >  NC_000913/1497178‑1497436

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