Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,897,439	G→A	*222* (TAG→TAA)	yieH →	phosphoenolpyruvate and 6‑phosphogluconate phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,897,439	0	G	A	100.0%	242.5 / NA	74	*222* (TAG→TAA)	yieH	phosphoenolpyruvate and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/37);  total (37/37)

AGGTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTG  >  NC_000913/3897307‑3897572                                                                                                                                     |                                                                                                                                      aGGTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAAttcttc                                                                                                                                 >  1:995941/1‑139 (MQ=255)    tATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTcaca                                                                                                                              <  2:166662/139‑1 (MQ=255)     aTGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCacac                                                                                                                             <  2:164159/139‑1 (MQ=255)         aaGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTccc                                                                                                                         <  2:995941/139‑1 (MQ=255)           gtgtTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCtt                                                                                                                       >  2:1210475/1‑139 (MQ=255)             gtTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCa                                                                                                                     >  2:21075/1‑139 (MQ=255)              ttCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCAc                                                                                                                    >  1:1069510/1‑139 (MQ=255)              ttCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCAc                                                                                                                    >  1:342147/1‑139 (MQ=255)               tCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACt                                                                                                                   >  2:1041820/1‑139 (MQ=255)                   cTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTAcc                                                                                                               >  1:1152433/1‑139 (MQ=255)                      cTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTcac                                                                                                                               <  2:23759/120‑1 (MQ=255)                      cTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTcac                                                                                                                               >  1:23759/1‑120 (MQ=255)                        gcgcCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCt                                                                                                          >  2:1357612/1‑139 (MQ=255)                               cccGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTg                                                                                                   <  2:1069510/139‑1 (MQ=255)                               cccGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTg                                                                                                   <  1:427739/139‑1 (MQ=255)                               cccGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTg                                                                                                   <  2:342147/139‑1 (MQ=255)                                   cacaATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGgcgc                                                                                               >  2:1274381/1‑139 (MQ=255)                                         aaGCCGATGGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATCGGCGCTCAAAg                                                                                         >  2:954896/1‑139 (MQ=255)                                          aGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAgg                                                                                        <  1:1210475/139‑1 (MQ=255)                                                         cGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGtt                                                                         <  2:1206155/139‑1 (MQ=255)                                                                 accaccTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGt                                                                 <  2:1315774/139‑1 (MQ=255)                                                                  ccaccTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTc                                                                >  1:376529/1‑139 (MQ=255)                                                                         tACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTccc                                                                                                                         >  2:784843/1‑75 (MQ=255)                                                                         tACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTccc                                                                                                                         <  1:784843/75‑1 (MQ=255)                                                                         tACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTAcc                                                                                                               >  2:931128/1‑85 (MQ=255)                                                                         tACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTAcc                                                                                                               <  1:931128/85‑1 (MQ=255)                                                                           cccATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTTGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACg                                                       <  1:1041820/139‑1 (MQ=255)                                                                           cccATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACg                                                       <  2:1152433/139‑1 (MQ=255)                                                                           cccATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACg                                                       <  2:855246/139‑1 (MQ=255)                                                                             cATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGCTATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTc                                                     <  1:954896/139‑1 (MQ=255)                                                                               tCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTaaa                                                                                <  1:1105763/110‑1 (MQ=255)                                                                               tCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTaaa                                                                                >  2:1105763/1‑110 (MQ=255)                                                                                cTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCt                                                                                              >  1:1029889/1‑95 (MQ=255)                                                                                cTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCt                                                                                              <  2:1029889/95‑1 (MQ=255)                                                                                cTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCggg                                                  >  2:510807/1‑139 (MQ=255)                                                                                cTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCggg                                                  >  2:1298554/1‑139 (MQ=255)                                                                                   tCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTaa                                               >  2:920091/1‑139 (MQ=255)                                                                                           aCCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTcac                                                                                                                               <  1:243468/51‑1 (MQ=255)                                                                                           aCCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTcac                                                                                                                               >  2:243468/1‑51 (MQ=255)                                                                                            ccTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACg                                      <  1:1274381/139‑1 (MQ=255)                                                                                                     tgGAAAGCGCGTGGTTGGGATATTACGGCATAAttc                                                                                                                                    <  1:243253/36‑1 (MQ=255)                                                                                                     tgGAAAGCGCGTGGTTGGGATATTACGGCATAAttc                                                                                                                                    >  2:243253/1‑36 (MQ=38)                                                                                                     tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTccc                                                                                                                         <  2:1298894/47‑1 (MQ=255)                                                                                                     tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTccc                                                                                                                         >  1:1298894/1‑47 (MQ=255)                                                                                                     tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGt                                                             >  2:433682/1‑107 (MQ=255)                                                                                                     tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGt                                                             <  1:433682/107‑1 (MQ=255)                                                                                                         aaGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGc                         <  1:1298554/139‑1 (MQ=255)                                                                                                            cgcgTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGgta                                                                                      <  2:1326190/75‑1 (MQ=255)                                                                                                            cgcgTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGgta                                                                                      >  1:1326190/1‑75 (MQ=255)                                                                                                               gTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAgg                                                                              >  1:993449/1‑80 (MQ=255)                                                                                                               gTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAg                                                                               <  2:350025/79‑1 (MQ=255)                                                                                                               gTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAg                                                                               >  1:350025/1‑79 (MQ=255)                                                                                                               gTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGt                                                      <  1:594848/104‑1 (MQ=255)                                                                                                               gTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGt                                                      >  2:594848/1‑104 (MQ=255)                                                                                                               gTGGTTGGGATATTACGGCATAATTCTGCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAgt                                                                              <  2:993449/80‑2 (MQ=255)                                                                                                                   ttGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGa               <  1:1357612/139‑1 (MQ=255)                                                                                                                   ttGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGa               >  1:1364017/1‑139 (MQ=255)                                                                                                                      ggATATTACGGCATAATTCTTCACACTCCCTTCACTTAcc                                                                                                               <  2:1404856/40‑1 (MQ=255)                                                                                                                      ggATATTACGGCATAATTCTTCACACTCCCTTCACTTAcc                                                                                                               >  1:1404856/1‑40 (MQ=255)                                                                                                                        atatTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGtt          <  2:1364017/139‑1 (MQ=255)                                                                                                                           ttACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTaaa                                                                                >  2:100829/1‑66 (MQ=255)                                                                                                                           ttACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTaaa                                                                                <  1:100829/66‑1 (MQ=255)                                                                                                                            tACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCt                                                                                              >  2:277819/1‑51 (MQ=255)                                                                                                                            tACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCt                                                                                              <  1:277819/51‑1 (MQ=255)                                                                                                                            tACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAgg                                                                                        >  2:604661/1‑57 (MQ=255)                                                                                                                            tACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAgg                                                                                        <  1:604661/57‑1 (MQ=255)                                                                                                                             aCGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTg     >  2:92102/1‑139 (MQ=255)                                                                                                                                gCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGgtg  >  1:231815/1‑139 (MQ=255)                                                                                                                                 cATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGgta                                                                                      >  1:213728/1‑54 (MQ=255)                                                                                                                                 cATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGgta                                                                                      <  2:213728/54‑1 (MQ=255)                                                                                                                                 cATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACg                                                       <  2:100455/85‑1 (MQ=255)                                                                                                                                 cATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACg                                                       >  1:100455/1‑85 (MQ=255)                                                                                                                                 cATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGa               >  1:1317561/1‑125 (MQ=255)                                                                                                                                 cATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGa               <  2:1317561/125‑1 (MQ=255)                                                                                                                                     |                                                                                                                                      AGGTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTG  >  NC_000913/3897307‑3897572

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