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

TATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCA  >  NC_000913/1497176‑1497427                                                                                                                                     |                                                                                                                        tatGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCaa                                                                                                                   >  1:668581/1‑139 (MQ=255)   tgggaatgggaTTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAg                                                                                                                 >  1:872436/1‑139 (MQ=255)      gaatgggaTTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGa                                                                                                              >  1:355703/1‑139 (MQ=255)                                      tCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAg                                                                              >  2:283919/1‑139 (MQ=255)                                          tAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGcccc                                                                          >  2:642787/1‑139 (MQ=255)                                          tAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGcccc                                                                          >  1:485145/1‑139 (MQ=255)                                            gATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCg                                                                        >  1:804603/1‑139 (MQ=255)                                               ccAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAAc                                                                     >  2:736537/1‑139 (MQ=255)                                               ccAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAAc                                                                     >  1:522931/1‑139 (MQ=255)                                               ccAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAAc                                                                     >  1:789674/1‑139 (MQ=255)                                                      cacaTCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGgcgcg                                                              <  2:804603/139‑1 (MQ=255)                                                         aTCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCg                                                           <  1:642787/139‑1 (MQ=255)                                                           cTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTg                                                         <  1:736537/139‑1 (MQ=255)                                                                 cgTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATgt                                                     >  1:723228/1‑137 (MQ=255)                                                                 cgTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATgt                                                     <  2:723228/137‑1 (MQ=255)                                                                   tGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACCATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATgtagta                                                 >  2:562911/1‑139 (MQ=255)                                                                   tGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATgtagta                                                 >  2:452832/1‑139 (MQ=255)                                                                                cGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGgcg                                    <  1:452832/139‑1 (MQ=255)                                                                                           gttCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTc                         >  2:454520/1‑139 (MQ=255)                                                                                                 aTACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGcc                   <  2:522931/139‑1 (MQ=255)                                                                                                 aTACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGcc                   <  1:90826/139‑1 (MQ=255)                                                                                                 aTACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGcc                   <  2:458836/139‑1 (MQ=255)                                                                                                      aCGATGCAGGTGTTGATTCAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTtga              <  1:572747/139‑1 (MQ=255)                                                                                                      aCGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCg                  >  2:916751/1‑135 (MQ=255)                                                                                                      aCGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCg                  <  1:916751/135‑1 (MQ=255)                                                                                                      aCGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTtga              <  2:477424/139‑1 (MQ=255)                                                                                                      aCGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTtga              <  1:463708/139‑1 (MQ=255)                                                                                                      aCGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTtga              <  2:789674/139‑1 (MQ=255)                                                                                                          tGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTg                                                         <  1:900993/92‑1 (MQ=255)                                                                                                          tGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTg                                                         >  2:900993/1‑92 (MQ=255)                                                                                                           gCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATa         <  1:454520/139‑1 (MQ=255)                                                                                                             aGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACa       <  1:283919/139‑1 (MQ=255)                                                                                                                  ttGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCa  >  1:501074/1‑139 (MQ=255)                                                                                                                    gATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATa         >  2:464112/1‑130 (MQ=255)                                                                                                                    gATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATa         <  1:464112/130‑1 (MQ=255)                                                                                                                         aaaaCAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGgcg                                    <  2:467824/98‑1 (MQ=255)                                                                                                                         aaaaCAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGgcg                                    >  1:467824/1‑98 (MQ=255)                                                                                                                                     |                                                                                                                        TATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCA  >  NC_000913/1497176‑1497427

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