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%	230.3 / NA	65	L284L (TTG→CTG)	glpQ	periplasmic glycerophosphodiester phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/32);  total (33/32)

GAGTGTATTCAGGCAGTTTATCTGACCGCACGGTATAAGGATGCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAACATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAACCCAGCTTCCATCCGGCTGTTTCTGCTGCGTTTCATTCCAG  >  NC_000913/2350039‑2350292                                                                                                                            |                                                                                                                                   gAGTGTATTCAGGCAGTTTATCTGACCGCACGGTATAAGGATGCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCgg                                                                                                                     >  1:858609/1‑139 (MQ=255)      tATTCAGGCAGTTTATCTGACCGCACGGTATAAGGATGCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCaa                                                                                                                >  1:1545814/1‑139 (MQ=255)              cAGTTTATCTGACCGCACGGTATAAGGATGCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCa                                                                                                                                  >  1:1012778/1‑113 (MQ=255)              cAGTTTATCTGACCGCACGGTATAAGGATGCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCa                                                                                                                                  <  2:1012778/113‑1 (MQ=255)                        gACCGCACGGTATAAGGATGCACTACCAGTTTATTCTGCTGAGCATCTTGCACCAAGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCg                                                                                              <  2:858609/139‑1 (MQ=255)                                           gCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAcc                                                                           <  1:983884/139‑1 (MQ=255)                                            cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGCTATTACCCGGCTGCGATGTCTCCTCAATCAGCa                                                                                                                                  >  1:975900/1‑83 (MQ=255)                                            cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCa                                                                                                                                  <  2:975900/83‑1 (MQ=255)                                            cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCa                                                                                                                                  <  1:371264/83‑1 (MQ=255)                                            cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATc                                                                                                                        <  1:768988/93‑1 (MQ=255)                                            cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATc                                                                                                                        >  2:768988/1‑93 (MQ=255)                                            cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAAt                                                                                                               >  2:1241512/1‑102 (MQ=255)                                            cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAAt                                                                                                               <  1:1241512/102‑1 (MQ=255)                                            cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATAc                                                                                                             >  2:599505/1‑104 (MQ=255)                                            cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATAc                                                                                                             <  1:599505/104‑1 (MQ=255)                                            cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGcc                                                                                            >  1:656583/1‑121 (MQ=255)                                            cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGcc                                                                                            <  2:656583/121‑1 (MQ=255)                                            cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAccc                                                                          >  2:609297/1‑139 (MQ=255)                                            cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAccc                                                                          >  2:168892/1‑139 (MQ=255)                                            cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAccc                                                                          >  2:391740/1‑139 (MQ=255)                                            cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAc                                                                            >  1:442811/1‑137 (MQ=255)                                            cACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAc                                                                            <  2:442811/137‑1 (MQ=255)                                                 ccAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCtt                                                                     >  1:304111/1‑139 (MQ=255)                                                      ttATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACa                                                                <  2:1545814/139‑1 (MQ=255)                                                         ttCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATcc                                                             <  2:304111/139‑1 (MQ=255)                                                              ctgAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATc                                                                                                                        >  2:272415/1‑75 (MQ=255)                                                              ctgAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATc                                                                                                                        <  1:272415/75‑1 (MQ=255)                                                              ctgAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCg                                                        <  1:168892/139‑1 (MQ=255)                                                              ctgAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCg                                                        >  1:236655/1‑139 (MQ=255)                                                               tgAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCg                                                                                              >  1:1126068/1‑100 (MQ=255)                                                               tgAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCg                                                                                              <  2:1126068/100‑1 (MQ=255)                                                                 aGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATc                                                                                                                        <  1:1148511/72‑1 (MQ=255)                                                                 aGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATc                                                                                                                        >  2:1148511/1‑72 (MQ=255)                                                                 aGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAg                                                     <  1:609297/139‑1 (MQ=255)                                                                   cATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGtt                                                   >  1:306875/1‑139 (MQ=255)                                                                           aCCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCAc                                                                                          >  2:525429/1‑92 (MQ=255)                                                                           aCCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCAc                                                                                          <  1:525429/92‑1 (MQ=255)                                                                           aCCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTaa                                           >  2:72834/1‑139 (MQ=255)                                                                             cATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGa                                                                                                                    >  1:228024/1‑64 (MQ=255)                                                                             cATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGa                                                                                                                    <  2:228024/64‑1 (MQ=255)                                                                               tGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAc                                                                            <  2:478973/102‑1 (MQ=255)                                                                               tGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAc                                                                            >  1:478973/1‑102 (MQ=255)                                                                               tGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAACCCa                                       >  2:636231/1‑139 (MQ=255)                                                                                gCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAAcc                                         >  1:1115627/1‑136 (MQ=255)                                                                                gCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAAcc                                         <  2:1115627/136‑1 (MQ=255)                                                                                   aGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAACCCAGCtt                                   <  1:636231/139‑1 (MQ=255)                                                                                       aGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAACCCAGCTTCCAt                               <  2:236655/139‑1 (MQ=255)                                                                                           tGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAACCCAGCTTCCATCCgg                           >  1:121201/1‑139 (MQ=255)                                                                                           tGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAACCCAGCTTCCATCCg                            >  1:400779/1‑138 (MQ=255)                                                                                           tGATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAACCCAGCTTCCATCCg                            <  2:400779/138‑1 (MQ=255)                                                                                            gATATTACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAACCCAGCTTCCATCCGGc                          >  2:24674/1‑139 (MQ=255)                                                                                                 tACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCAt                                                                                                          <  1:362116/54‑1 (MQ=255)                                                                                                 tACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCAt                                                                                                          >  2:362116/1‑54 (MQ=255)                                                                                                 tACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATcc                                                             >  2:1192498/1‑99 (MQ=255)                                                                                                 tACCCGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATcc                                                             <  1:1192498/99‑1 (MQ=255)                                                                                                   cccGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAcc                                                                           >  1:1582824/1‑83 (MQ=255)                                                                                                   cccGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAcc                                                                           <  2:1582824/83‑1 (MQ=255)                                                                                                   cccGGCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAACCCAGCTTCCATCCGGCTGTTTct                   >  2:782162/1‑139 (MQ=255)                                                                                                       gCTGCGATGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAACCCAGCTTCCATCCGGCTGTTTctgctg               <  2:121201/139‑1 (MQ=255)                                                                                                             aTGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAcc                                                                           <  2:278701/73‑1 (MQ=255)                                                                                                             aTGTCTCCTCAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAcc                                                                           >  1:278701/1‑73 (MQ=255)                                                                                                                  tcctcAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCg                                                        >  1:930617/1‑87 (MQ=255)                                                                                                                  tcctcAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCg                                                        <  2:930617/87‑1 (MQ=255)                                                                                                                    ctcAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAACCCAGCTTCCATCCGGCTGTTTCTGCTGCGTTTCATTCCAg  <  1:391740/139‑1 (MQ=255)                                                                                                                    ctcAATCAGCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAACCCAGCTTCCATCCGGCTGTTTCTGCTGCGTTTCATTCCAg  <  2:1271586/139‑1 (MQ=255)                                                                                                                         tcagCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAc                                                                            <  1:1495142/60‑1 (MQ=255)                                                                                                                         tcagCATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCAc                                                                            >  2:1495142/1‑60 (MQ=255)                                                                                                                            |                                                                                                                                   GAGTGTATTCAGGCAGTTTATCTGACCGCACGGTATAAGGATGCACTACCAGTTTATTCTGCTGAGCATCTTGCACCATGCCAGTGAGTTTGATATTACCCGGCTGCGATGTCTCCTCAATCAACATATGGTAATCCGGACCAATACCATCTGCATATTCCGCCACCTGTTTCATGGCACCCGGCTTAAACATCCAGTCGTAGTTGTAATTAACCCAGCTTCCATCCGGCTGTTTCTGCTGCGTTTCATTCCAG  >  NC_000913/2350039‑2350292

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