Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,008,075	A→G	100%	D189G (GAC→GGC)	ygeX →	2,3‑diaminopropionate ammonia lyase, PLP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,008,075	0	A	G	100.0%	272.7 / ‑6.7	80	D189G (GAC→GGC)	ygeX	2,3‑diaminopropionate ammonia lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (46/33);  total (47/33)

CTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGG  >  NC_000913/3007978‑3008170                                                                                                  |                                                                                                cTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACACGCAACAGCACAGCTAAGAAATGGTACAGGgca                                                                                               >  1:1040077/1‑100 (MQ=255)      tGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCa                                                                                          <  2:1369933/100‑1 (MQ=255)         cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                       <  1:488770/100‑1 (MQ=255)         cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATg                                                                                        >  1:945073/1‑99 (MQ=255)         cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTCACCATGCAACACGCGAAGCAGAAAGGCAGGGAAGTGGTACAGGGCACGGCATg                                                                                        >  1:917790/1‑99 (MQ=255)          gAGTGCATCGTCACGGATATGAACTATGGCGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATggg                                                                                      <  2:97674/100‑1 (MQ=255)              gCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAgg                                                                                  <  2:1111498/100‑1 (MQ=255)              gCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAgg                                                                                  <  1:2005874/100‑1 (MQ=255)               cATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGt                                                                                 <  2:1368180/100‑1 (MQ=255)               cATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGt                                                                                 >  2:1270873/1‑100 (MQ=255)                  cGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTac                                                                              >  1:1658955/1‑100 (MQ=255)                  cGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTac                                                                              >  2:378573/1‑100 (MQ=255)                   gTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTaca                                                                             >  1:1416940/1‑100 (MQ=255)                   gTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTaca                                                                             >  2:1783702/1‑100 (MQ=255)                   gTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTaca                                                                             >  2:2015991/1‑100 (MQ=255)                   gTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTaca                                                                             >  2:2657331/1‑100 (MQ=255)                     cACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACAcc                                                                           >  1:1765744/1‑100 (MQ=255)                       cGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCaa                                                                         >  2:2036774/1‑100 (MQ=255)                       cGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCaa                                                                         <  1:588580/100‑1 (MQ=255)                       cGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGAAGCACGGCTGGGAACTGGTACAGGGCACGGCATGGGAAGGTTACACaaa                                                                         >  2:764454/1‑100 (MQ=255)                         gATATGAACTATGACGATACCGTTCGCCTGACCATGCCACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGAAAGGTTACACCaaaa                                                                       >  1:424341/1‑100 (MQ=255)                         gATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCaaaa                                                                       <  2:585937/100‑1 (MQ=255)                           tatGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATc                                                                     <  2:2228163/100‑1 (MQ=255)                           tatGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATc                                                                     <  1:1630746/100‑1 (MQ=255)                             tGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATccc                                                                   <  1:1253939/100‑1 (MQ=255)                             tGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATccc                                                                   >  2:899544/1‑100 (MQ=255)                             tGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATccc                                                                   >  1:809931/1‑100 (MQ=255)                              gAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAAGACACGGCATGGGAAGGTTACACCAAAATCCCa                                                                  >  1:1385038/1‑100 (MQ=255)                               ccgatctaCGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCaa                                                                 <  2:1082213/93‑1 (MQ=255)                                aCTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAAc                                                                >  1:699993/1‑100 (MQ=255)                                aCTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAAc                                                                <  2:1062979/100‑1 (MQ=255)                                     gACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGa                                                           >  1:2121791/1‑100 (MQ=255)                                      aCGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCaaagatcgg                                                          >  1:1082223/1‑93 (MQ=255)                                      aCGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGAt                                                          >  2:994608/1‑100 (MQ=255)                                       cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                         >  1:36029/1‑100 (MQ=255)                                       cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCAAGGCATGGGAAGGATACACCAAAATACCAACCTGGa                                                           >  1:336069/1‑98 (MQ=255)                                         aTACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCAt                                                       <  2:2459080/100‑1 (MQ=255)                                              gTTCGCCTGACGATGCAACACGCGCAGCAGCCGGGCGGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACTAAAATCCCAACCTGGATCATGCAAg                                                  <  2:2251083/100‑1 (MQ=255)                                              gTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAg                                                  <  1:46132/100‑1 (MQ=255)                                                tCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGc                                                <  1:928515/100‑1 (MQ=255)                                                 cGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCt                                               >  1:2100115/1‑100 (MQ=255)                                                 cGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCt                                               >  2:1052675/1‑100 (MQ=255)                                                 cGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCt                                               >  2:1392121/1‑100 (MQ=255)                                                  gCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTa                                              >  1:1254096/1‑100 (MQ=255)                                                  gCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTa                                              >  2:378712/1‑100 (MQ=255)                                                  gCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTa                                              >  1:483033/1‑100 (MQ=255)                                                   ccTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGTAAAAGTAACAACAAAATCCCAACCTAGATCATGCAAAGCTAc                                             >  1:1581509/1‑100 (MQ=255)                                                      gACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCa                                          <  1:507238/100‑1 (MQ=255)                                                        ccATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAc                                        >  2:1782118/1‑100 (MQ=255)                                                            gcagcaCTCGTAGCAGCAAGGCTGGGAAGGGGTACAGGGCACAGCATGGGAAGGTTACACCAAAATCCCAACATGGATCATACAAGGCTACGCAACCCTg                                    >  2:417850/1‑100 (MQ=255)                                                            gCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTg                                    >  1:569803/1‑100 (MQ=255)                                                            gCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTg                                    <  2:2083270/100‑1 (MQ=255)                                                               acacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCa                                 <  1:191289/100‑1 (MQ=255)                                                               acacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCa                                 <  1:968262/100‑1 (MQ=255)                                                                cacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAg                                >  1:419058/1‑100 (MQ=255)                                                                caagcgtAGCAGTACGGCTGGGAAGTGGTATAGAGTACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAg                                <  1:168148/97‑1 (MQ=255)                                                                 acGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGa                               >  2:2290398/1‑100 (MQ=255)                                                                  cgcgCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGAt                              >  2:93217/1‑100 (MQ=255)                                                                   gcgcAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATg                             <  2:2626110/100‑1 (MQ=255)                                                                       agcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGc                         <  1:361479/100‑1 (MQ=255)                                                                        gcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGcc                        >  2:1981976/1‑100 (MQ=255)                                                                         cagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCg                       >  1:553237/1‑100 (MQ=255)                                                                           gcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGAATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTc                     >  1:678979/1‑100 (MQ=255)                                                                              cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAg                  <  2:1836767/100‑1 (MQ=255)                                                                              cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAg                  >  2:2387888/1‑100 (MQ=255)                                                                                gCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGGTGAAGCCGTCGAGCa                <  2:1015796/100‑1 (MQ=255)                                                                                gCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCa                <  1:43827/100‑1 (MQ=255)                                                                                   gggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAAt             >  2:1202441/1‑100 (MQ=255)                                                                                    ggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACACAACCCTGGCAGATGAAGCCGTCGAGCAAATg            >  1:699998/1‑100 (MQ=255)                                                                                      aaGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCg          >  2:390962/1‑100 (MQ=255)                                                                                      aaGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCg          <  2:665164/100‑1 (MQ=255)                                                                                        gTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTg        <  1:378575/100‑1 (MQ=255)                                                                                        gTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTg        <  1:1297630/100‑1 (MQ=255)                                                                                          ggTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGTTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaa      <  1:698469/100‑1 (MQ=255)                                                                                          ggTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaa      >  2:1206550/1‑100 (MQ=255)                                                                                          ggTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaa      >  2:2639458/1‑100 (MQ=255)                                                                                           gTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaaa     >  2:920634/1‑100 (MQ=255)                                                                                            tACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAAt    <  2:424338/100‑1 (MQ=255)                                                                                             aCAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaaa     >  2:1514149/1‑98 (MQ=255)                                                                                             aCAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATg   <  1:2657363/100‑1 (MQ=255)                                                                                              cAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAAACTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATgg  >  2:1420705/1‑100 (MQ=255)                                                                                                  |                                                                                                CTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGG  >  NC_000913/3007978‑3008170

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑