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%	166.9 / NA	47	E151E (GAA→GAG)	opgD	OPG biosynthetic periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/27);  total (20/27)

GGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCA  >  NC_000913/1497186‑1497427                                                                                                                           |                                                                                                                        ggaTTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCg                                                                                                         >  1:544377/1‑139 (MQ=255)   aTTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGc                                                                                                       >  2:13549/1‑139 (MQ=255)     tCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCtt                                                                                                     >  2:44855/1‑139 (MQ=255)      ccGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCttt                                                                                                    >  2:506968/1‑139 (MQ=255)       cGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTg                                                                                                   >  2:805210/1‑139 (MQ=255)           gccgcGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCgg                                                                                               >  2:13225/1‑139 (MQ=255)            ccgcGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGt                                                                                              >  2:602145/1‑139 (MQ=255)                  tCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTcgc                                                                                        <  2:544377/139‑1 (MQ=255)                                        gCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGTTTTTCGCGTGTTTAAAGCCCCCGAACTgg                                                                  >  2:692594/1‑139 (MQ=255)                                        gCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTgg                                                                  >  1:126522/1‑139 (MQ=255)                                        gCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTgg                                                                  <  1:73151/139‑1 (MQ=255)                                          aaCACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGgcg                                                                <  2:812380/139‑1 (MQ=255)                                          aaCACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGgcg                                                                <  1:33923/139‑1 (MQ=255)                                          aaCACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGgcg                                                                <  2:506011/139‑1 (MQ=255)                                          aaCACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGgcg                                                                <  2:233242/139‑1 (MQ=255)                                            cacaTCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGgcgcg                                                              <  1:13225/139‑1 (MQ=255)                                              caTCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGcc                                                            >  1:831059/1‑139 (MQ=255)                                              caTCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGcc                                                            <  2:954155/139‑1 (MQ=255)                                              caTCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGcc                                                            <  2:554703/139‑1 (MQ=255)                                               aTCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCg                                                           <  2:126522/139‑1 (MQ=255)                                                   ggCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTa                                                                                 <  2:984260/113‑1 (MQ=255)                                                   ggCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTa                                                                                 >  1:984260/1‑113 (MQ=255)                                                    gcgcGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATg                                                      <  1:44855/139‑1 (MQ=255)                                                                      cGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTg                                                         <  1:454833/118‑1 (MQ=255)                                                                      cGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTg                                                         >  2:454833/1‑118 (MQ=255)                                                                      cGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGgcg                                    <  1:805210/139‑1 (MQ=255)                                                                                    cAACTACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCcgc                      <  1:692594/139‑1 (MQ=255)                                                                                       aTACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGcc                   <  1:13549/139‑1 (MQ=255)                                                                                            aCGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGa                                                        <  2:377630/97‑1 (MQ=255)                                                                                            aCGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGa                                                        >  1:377630/1‑97 (MQ=255)                                                                                            aCGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGcc                   >  1:584198/1‑134 (MQ=255)                                                                                            aCGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGcc                   <  2:584198/134‑1 (MQ=255)                                                                                            aCGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTtga              <  2:762218/139‑1 (MQ=255)                                                                                                tGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTgatgat          <  2:831059/139‑1 (MQ=255)                                                                                                 gCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTg                                                                   >  2:429097/1‑81 (MQ=255)                                                                                                 gCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTg                                                                   <  1:429097/81‑1 (MQ=255)                                                                                                  cAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTtgatga           <  1:70292/136‑1 (MQ=255)                                                                                                  cAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTtgatga           >  2:70292/1‑136 (MQ=255)                                                                                                        ttGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCa  <  1:602145/139‑1 (MQ=255)                                                                                                          gATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGgcg                                                                >  2:537419/1‑75 (MQ=255)                                                                                                          gATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGgcg                                                                <  1:537419/75‑1 (MQ=255)                                                                                                                aaaCAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGgcg                                                                <  2:572929/69‑1 (MQ=255)                                                                                                                aaaCAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGgcg                                                                >  1:572929/1‑69 (MQ=255)                                                                                                                 aaCAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATg                                                      <  2:80084/78‑1 (MQ=255)                                                                                                                 aaCAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATg                                                      >  1:80084/1‑78 (MQ=255)                                                                                                                      ttAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGcc                   <  2:853377/108‑1 (MQ=255)                                                                                                                      ttAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGcc                   >  1:853377/1‑108 (MQ=255)                                                                                                                           |                                                                                                                        GGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCA  >  NC_000913/1497186‑1497427

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