Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	607,988	G→A	*51* (TAG→TAA)	hokE →	toxic polypeptide, small

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	607,988	0	G	A	100.0%	208.8 / NA	64	*51* (TAG→TAA)	hokE	toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/30);  total (34/30)

TGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAA  >  NC_000913/607856‑608121                                                                                                                                     |                                                                                                                                      tGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTtg                                                                                                                                 >  2:390030/1‑139 (MQ=255)       aGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCggg                                                                                                                           <  1:275698/139‑1 (MQ=255)           aTAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGt                                                                                                                       <  2:306919/139‑1 (MQ=255)            tAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTa                                                                                                                      <  2:713006/139‑1 (MQ=255)                   gtgtTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTGCGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCAt                                                                                                               >  1:344006/1‑139 (MQ=255)                   gtgtTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCAt                                                                                                               >  1:730548/1‑139 (MQ=255)                   gtgtTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCAt                                                                                                               >  2:395709/1‑139 (MQ=255)                    tgtTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATa                                                                                                              >  1:234286/1‑139 (MQ=255)                          aCGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCt                                                                                                        <  1:390355/139‑1 (MQ=255)                           cGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTa                                                                                                       >  2:102942/1‑139 (MQ=255)                                cTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTa                                                                                                  <  1:102942/139‑1 (MQ=255)                                cTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTa                                                                                                  >  2:762810/1‑139 (MQ=255)                                cTGGGATTTACGCTTCTGGTCGGAGACTCACTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTa                                                                                                  >  1:514645/1‑139 (MQ=255)                                    gATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                                              <  1:757373/139‑1 (MQ=255)                                    gATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                                              <  1:117653/139‑1 (MQ=255)                                    gATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACCTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                                              <  2:514645/139‑1 (MQ=255)                                        tACGCTTCTGGTCGGTGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAg                                                                                                >  2:515060/1‑133 (MQ=255)                                        tACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAg                                                                                                <  1:515060/133‑1 (MQ=255)                                        tACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtgt                                                                                          >  1:611033/1‑139 (MQ=255)                                               cTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                                                                   >  1:703700/1‑139 (MQ=255)                                                  gTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                                                                <  2:654037/139‑1 (MQ=255)                                                  gTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                                                                <  1:16846/139‑1 (MQ=255)                                                   tCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGcc                                                                               <  1:395709/139‑1 (MQ=255)                                                       agaCTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                                                                           <  2:730548/139‑1 (MQ=255)                                                         aCTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGatat                                                                         >  1:6924/1‑139 (MQ=255)                                                                 tgtgAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                                              >  2:497919/1‑110 (MQ=255)                                                                 tgtgAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                                              <  1:497919/110‑1 (MQ=255)                                                                       tttACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                                                           <  2:234286/139‑1 (MQ=255)                                                                         tACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt                                                         >  2:159767/1‑139 (MQ=255)                                                                                 aGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                                                 >  1:185009/1‑139 (MQ=255)                                                                                        tAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGa                                          >  1:692163/1‑139 (MQ=255)                                                                                        tAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGa                                          >  2:508552/1‑139 (MQ=255)                                                                                                aGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                                              >  1:504414/1‑79 (MQ=255)                                                                                                aGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                                              <  2:504414/79‑1 (MQ=255)                                                                                                aGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGt                                  >  2:619080/1‑139 (MQ=255)                                                                                                aGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGt                                  >  1:597291/1‑139 (MQ=255)                                                                                                 gTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg                                                        >  2:150484/1‑116 (MQ=255)                                                                                                 gTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg                                                        <  1:150484/116‑1 (MQ=255)                                                                                                 gTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTc                                 >  1:274295/1‑139 (MQ=255)                                                                                                  tttAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTg                                    >  1:91870/1‑135 (MQ=255)                                                                                                  tttAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTg                                    <  2:91870/135‑1 (MQ=255)                                                                                                  tttAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCa                                >  1:402452/1‑139 (MQ=255)                                                                                                     aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCa                             >  2:640389/1‑139 (MQ=255)                                                                                                     aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCa                             >  1:548948/1‑139 (MQ=255)                                                                                                        gCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAAtt                          <  2:692163/139‑1 (MQ=255)                                                                                                        gCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAAtt                          <  2:344006/139‑1 (MQ=255)                                                                                                          tGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCt                        >  1:709564/1‑139 (MQ=255)                                                                                                           gTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                                                                  >  2:94948/1‑96 (MQ=255)                                                                                                           gTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                                                                  >  2:420782/1‑96 (MQ=255)                                                                                                           gTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                                                                  <  1:94948/96‑1 (MQ=255)                                                                                                           gTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                                                                  <  1:420782/96‑1 (MQ=255)                                                                                                           gTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTg                       <  2:597291/139‑1 (MQ=255)                                                                                                            ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtgt                                                                                          >  2:182327/1‑71 (MQ=255)                                                                                                            ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtgt                                                                                          <  1:182327/71‑1 (MQ=255)                                                                                                            ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTgg                      >  2:53533/1‑139 (MQ=255)                                                                                                             tctcGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc                     <  2:548948/139‑1 (MQ=255)                                                                                                             tctcGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc                     <  1:719642/139‑1 (MQ=255)                                                                                                                  cTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGACCata                <  1:390030/139‑1 (MQ=255)                                                                                                                       gAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGt           <  1:159767/139‑1 (MQ=255)                                                                                                                           cGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGc       >  2:762917/1‑139 (MQ=255)                                                                                                                             aagaagTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                                                 >  2:117101/1‑95 (MQ=255)                                                                                                                             aagaagTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                                                 <  1:117101/95‑1 (MQ=255)                                                                                                                             aagaagTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGccc     >  2:347568/1‑139 (MQ=255)                                                                                                                                aagTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaa  <  1:762810/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                      TGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAA  >  NC_000913/607856‑608121

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