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%	228.0 / NA	70	*51* (TAG→TAA)	hokE	toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/38);  total (32/38)

CCCTTGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGC  >  NC_000913/607852‑608102                                                                                                                                         |                                                                                                                   cccTTGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAAcc                                                                                                                  <  2:10951/139‑1 (MQ=255)    ttGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGtt                                                                                                               <  1:616822/139‑1 (MQ=255)    ttGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGtt                                                                                                               <  2:16122/139‑1 (MQ=255)     tGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTtg                                                                                                              >  2:293247/1‑139 (MQ=255)         gCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCg                                                                                                          >  2:853937/1‑139 (MQ=255)         gCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCg                                                                                                          <  1:289910/139‑1 (MQ=255)           aGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCggg                                                                                                        <  1:293247/139‑1 (MQ=255)                   tGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCggg                                                                                                        >  2:984200/1‑131 (MQ=255)                   tGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCggg                                                                                                        <  1:984200/131‑1 (MQ=255)                       gtgtTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCAt                                                                                            >  2:314655/1‑139 (MQ=255)                        tgtTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATa                                                                                           >  2:1060384/1‑139 (MQ=255)                        tgtTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATa                                                                                           >  2:1025738/1‑139 (MQ=255)                                  tGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTaa                                                                                   >  1:200346/1‑137 (MQ=255)                                  tGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTaa                                                                                   <  2:200346/137‑1 (MQ=255)                                      gggATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAg                                                                             <  1:66639/139‑1 (MQ=255)                                        gATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                           <  1:314655/139‑1 (MQ=255)                                           ttACGCTTCTGGTCGGAGTCTGACTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtg                                                                        <  1:853937/139‑1 (MQ=255)                                                      gTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                                             <  1:1025738/139‑1 (MQ=255)                                                      gTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                                             >  1:258610/1‑139 (MQ=255)                                                       tCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGcc                                                            <  1:337709/139‑1 (MQ=255)                                                           agaCTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCggg                                                                                                        <  1:58278/91‑1 (MQ=255)                                                           agaCTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCggg                                                                                                        >  2:58278/1‑91 (MQ=255)                                                           agaCTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCt                                                                                     <  1:603538/110‑1 (MQ=255)                                                           agaCTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCt                                                                                     >  2:603538/1‑110 (MQ=255)                                                           agaCTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTaa                                                                              >  1:130290/1‑117 (MQ=255)                                                           agaCTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTaa                                                                              <  2:130290/117‑1 (MQ=255)                                                           aaaCTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                                                        >  2:1059438/3‑139 (MQ=255)                                                                cGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTATGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                                             <  1:530648/129‑1 (MQ=255)                                                                cGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTATGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                                             >  2:530648/1‑129 (MQ=255)                                                                cGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACCTAAGGGTTGTGGTATTACGCCTGATATGAt                                                   >  2:775506/1‑139 (MQ=255)                                                                     tgtgAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATa                                                                                           >  2:914175/1‑94 (MQ=255)                                                                     tgtgAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATa                                                                                           <  1:914175/94‑1 (MQ=255)                                                                          gTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGcc                                         <  1:775506/139‑1 (MQ=255)                                                                               cGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                                              <  2:949989/113‑1 (MQ=255)                                                                               cGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                                              >  1:949989/1‑113 (MQ=255)                                                                                    aaGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGt                                                                                                    >  2:530184/1‑70 (MQ=255)                                                                                    aaGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGt                                                                                                    <  1:530184/70‑1 (MQ=255)                                                                                     aGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                              <  2:219544/139‑1 (MQ=255)                                                                                     aGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                              <  2:258610/139‑1 (MQ=255)                                                                                        aaCGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc                           <  1:1059438/139‑1 (MQ=255)                                                                                            tAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGa                       >  2:690845/1‑139 (MQ=255)                                                                                               tatTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAgcg                                                                                       >  1:67535/1‑72 (MQ=255)                                                                                               tatTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAgcg                                                                                       <  2:67535/72‑1 (MQ=255)                                                                                                    aGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCg                                                                                                          <  2:1012436/48‑1 (MQ=255)                                                                                                    aGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCg                                                                                                          >  1:1012436/1‑48 (MQ=255)                                                                                                    aGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGt               >  2:524446/1‑139 (MQ=255)                                                                                                       ttAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAg            <  1:1060384/139‑1 (MQ=255)                                                                                                         aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGt                                                                                                                >  1:1037504/1‑37 (MQ=255)                                                                                                         aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGt                                                                                                                <  2:1037504/37‑1 (MQ=255)                                                                                                         aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCg                                                                                                          <  2:710783/43‑1 (MQ=255)                                                                                                         aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCg                                                                                                          >  1:710783/1‑43 (MQ=255)                                                                                                         aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCt                                                                                     <  2:377825/64‑1 (MQ=255)                                                                                                         aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCt                                                                                     >  1:377825/1‑64 (MQ=255)                                                                                                         aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTaa                                                                              <  1:446799/71‑1 (MQ=255)                                                                                                         aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTaa                                                                              >  2:446799/1‑71 (MQ=255)                                                                                                         aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtg                                                                        <  2:867483/77‑1 (MQ=255)                                                                                                         aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtg                                                                        >  1:867483/1‑77 (MQ=255)                                                                                                              tGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCt     >  1:865331/1‑139 (MQ=255)                                                                                                                ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCggg                                                                                                        >  1:913503/1‑38 (MQ=255)                                                                                                                ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCggg                                                                                                        <  2:913503/38‑1 (MQ=255)                                                                                                                ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                                        <  1:945025/102‑1 (MQ=255)                                                                                                                ttCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                                        >  2:945025/1‑102 (MQ=255)                                                                                                                 tctcGCTTACGAACCGAAGCAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc  <  2:727014/139‑1 (MQ=255)                                                                                                                 tctcGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc  >  1:146820/1‑139 (MQ=255)                                                                                                                 tctcGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc  <  1:524446/139‑1 (MQ=255)                                                                                                                 tctcGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc  <  2:146820/139‑1 (MQ=255)                                                                                                                          cGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                              >  2:1070405/1‑102 (MQ=255)                                                                                                                          cGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                              <  1:1070405/102‑1 (MQ=255)                                                                                                                                 aagaagTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc  <  1:645277/123‑1 (MQ=255)                                                                                                                                 aagaagTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc  >  2:645277/1‑123 (MQ=255)                                                                                                                                         |                                                                                                                   CCCTTGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGC  >  NC_000913/607852‑608102

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: ‑