Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,350,162	A→G	L284L (TTG→CTG)	glpQ ←	periplasmic glycerophosphodiester phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,350,162	0	A	G	100.0%	223.8 / NA	62	L284L (TTG→CTG)	glpQ	periplasmic glycerophosphodiester phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/32);  total (30/32)

GGAGTGTATTCAGGCAGTTTATCTGACCGCACGGTATAAGGATGCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAACATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAACCCAGCTTCCATCCGGCTGTTTCTGCTGCGTTTCATTCCAG  >  NC_000913/2350038‑2350292                                                                                                                             |                                                                                                                                   ggAGTGTATTCAGGCAGTTTATCTGACCGCACGGTATAAGGATGCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCg                                                                                                                      <  2:1089926/139‑1 (MQ=255)       tATTCAGGCAGTTTATCTGACCGCACGGTATAAGGATGCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCa                                                                                                                                  <  2:210011/121‑1 (MQ=255)       tATTCAGGCAGTTTATCTGACCGCACGGTATAAGGATGCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCaa                                                                                                                >  1:1469304/1‑139 (MQ=255)           cAGGCAGTTTATCTGACCGCACGGTATAAGGATGCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATAcc                                                                                                            <  2:1469304/139‑1 (MQ=255)             ggCAGTTTATCTGACCGCACGGTATAAGGATGCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCAt                                                                                                          >  1:788882/1‑139 (MQ=255)                 gTTTATCTGACCGCACGGTATAAGGATGCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGc                                                                                                      >  2:650447/1‑139 (MQ=255)                             gCACGGTATAAGGATGCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCAc                                                                                          <  1:650447/139‑1 (MQ=255)                             gCACGGTATAAGGATGCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCAc                                                                                          <  1:631140/139‑1 (MQ=255)                              cACGGTATAAGGATGCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCAcc                                                                                         >  1:290331/1‑139 (MQ=255)                                           tGCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAc                                                                            <  2:290331/139‑1 (MQ=255)                                            gCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAcc                                                                           <  2:788882/139‑1 (MQ=255)                                            gCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAcc                                                                           <  1:352219/139‑1 (MQ=255)                                             ctcTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCa                                                                                  <  2:1036314/129‑1 (MQ=255)                                             cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATAcc                                                                                                            >  1:44095/1‑105 (MQ=255)                                             cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATAcc                                                                                                            <  2:44095/105‑1 (MQ=255)                                             cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCa                                                                                  >  1:1036314/1‑131 (MQ=255)                                             cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCAt                                                                                 <  1:243948/132‑1 (MQ=255)                                             cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCAt                                                                                 <  2:1466198/132‑1 (MQ=255)                                             cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCAt                                                                                 >  2:243948/1‑132 (MQ=255)                                             cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCAt                                                                                 >  1:1466198/1‑132 (MQ=255)                                             cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAccc                                                                          >  1:250301/1‑139 (MQ=255)                                             cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAccc                                                                          >  2:1152561/1‑139 (MQ=255)                                             cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAccc                                                                          >  1:933537/1‑139 (MQ=255)                                             cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAccc                                                                          >  1:1348620/1‑139 (MQ=255)                                             cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAccc                                                                          >  2:165193/1‑139 (MQ=255)                                             cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAccc                                                                          >  1:674857/1‑139 (MQ=255)                                             cACTACCAGTTTATTCTGCTAAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATAc                                                                                                             <  1:93368/104‑1 (MQ=255)                                             cACTACCAGTTTATTCTGCTAAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATAc                                                                                                             >  2:93368/1‑104 (MQ=255)                                               cTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCgg                                                                        <  2:250301/139‑1 (MQ=255)                                                tACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGc                                                                       <  2:1348620/139‑1 (MQ=255)                                                   cAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTa                                                                    <  2:933537/139‑1 (MQ=255)                                                             tgctgAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGt                                                          <  2:674857/139‑1 (MQ=255)                                                             tgctgAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGt                                                          <  1:1152561/139‑1 (MQ=255)                                                               ctgAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATc                                                                                                                        <  1:393583/75‑1 (MQ=255)                                                               ctgAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATc                                                                                                                        >  2:393583/1‑75 (MQ=255)                                                               ctgAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCg                                                        >  2:915913/1‑139 (MQ=255)                                                                tgAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAAt                                                                                                               >  2:806605/1‑83 (MQ=255)                                                                tgAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAAt                                                                                                               <  1:806605/83‑1 (MQ=255)                                                                  aGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAg                                                     <  2:818491/139‑1 (MQ=255)                                                                  aGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAg                                                     <  2:760233/139‑1 (MQ=255)                                                                         tGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATgg                                                                                                                             <  1:930464/60‑1 (MQ=255)                                                                         tGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATgg                                                                                                                             >  2:930464/1‑60 (MQ=255)                                                                                tGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCa                                                                                                                                  >  1:279132/1‑48 (MQ=255)                                                                                tGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGc                                                                                             >  2:709256/1‑85 (MQ=255)                                                                                tGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGc                                                                                             <  1:709256/85‑1 (MQ=255)                                                                                tGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAACCCa                                       >  2:605910/1‑139 (MQ=255)                                                                                          tttGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAACCCAGCTTCCATcc                             >  2:506498/1‑139 (MQ=255)                                                                                           ttGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAACCCAGCTTCCATCCg                            >  1:1267457/1‑139 (MQ=255)                                                                                            tGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGtt                                                   <  1:1397340/115‑1 (MQ=255)                                                                                            tGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGtt                                                   >  2:1397340/1‑115 (MQ=255)                                                                                                   aCCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGcc                                                                                            <  1:186324/67‑1 (MQ=255)                                                                                                   aCCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGcc                                                                                            >  2:186324/1‑67 (MQ=255)                                                                                                    cccGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTg                                                                                                       >  2:1498793/1‑55 (MQ=255)                                                                                                    cccGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTg                                                                                                       <  1:1498793/55‑1 (MQ=255)                                                                                                        gCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCg                                                        >  2:938126/1‑98 (MQ=255)                                                                                                        gCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCg                                                        <  1:938126/98‑1 (MQ=255)                                                                                                         cTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATc                                                              >  2:740139/1‑91 (MQ=255)                                                                                                         cTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATc                                                              <  1:740139/91‑1 (MQ=255)                                                                                                                 tctcCTCAATCAGCATATGGTAATCCGGACCAATACCATATGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAACCCAGCTTCCATCCGGCTGTTTGTGCTGCGTTTCAtt      <  1:915913/139‑1 (MQ=255)                                                                                                                   tcctcAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCAc                                                                                          <  1:82117/53‑1 (MQ=255)                                                                                                                   tcctcAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCAc                                                                                          >  2:82117/1‑53 (MQ=255)                                                                                                                     ctcAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAACCCAGCTTCCATCCGGCTGTTTCTGCTGCGTTTCATTCCAg  <  1:165193/139‑1 (MQ=255)                                                                                                                      tcaatcagCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAcc                                                                           <  2:274994/65‑1 (MQ=255)                                                                                                                      tcaatcagCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAcc                                                                           >  1:274994/1‑65 (MQ=255)                                                                                                                             |                                                                                                                                   GGAGTGTATTCAGGCAGTTTATCTGACCGCACGGTATAAGGATGCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAACATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAACCCAGCTTCCATCCGGCTGTTTCTGCTGCGTTTCATTCCAG  >  NC_000913/2350038‑2350292

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