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

TGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAA  >  NC_000913/607856‑608114                                                                                                                                     |                                                                                                                               tGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTtg                                                                                                                          >  1:766341/1‑139 (MQ=255)           aTAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGt                                                                                                                <  1:9607/139‑1 (MQ=255)                   gtgtTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCAt                                                                                                        >  2:242712/1‑139 (MQ=255)                    tgtTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATa                                                                                                       >  2:465282/1‑139 (MQ=255)                        tAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAgcg                                                                                                   <  1:984993/139‑1 (MQ=255)                        tAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAgcg                                                                                                   <  1:443906/139‑1 (MQ=255)                          aCGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCt                                                                                                 <  2:217052/139‑1 (MQ=255)                                    gATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                                       <  1:120275/139‑1 (MQ=255)                                    gATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                                       <  2:766341/139‑1 (MQ=255)                                        tACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtgt                                                                                   >  1:1135139/1‑139 (MQ=255)                                        tACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAAAAATGAGTTTAAGGCTGGTCTCGATTACGAACGAAAGAAGTAAACGTTGTGCGGGGAGTAATCGAATAAGAGCAAACTTAAaggaagt                                                                                   >  1:422380/1‑135 (MQ=255)                                        tACGCTTCTGGTCGGAGACTCGCTGTGTAAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtgt                                                                                   >  1:385182/1‑139 (MQ=255)                                          cGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTgg                                                                                 >  1:316465/1‑139 (MQ=255)                                             ttCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                                                              >  2:372991/1‑139 (MQ=255)                                                  gTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTaa                                                                                               >  2:620956/1‑117 (MQ=255)                                                  gTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTaa                                                                                               <  1:620956/117‑1 (MQ=255)                                                       agaCTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                                                                    <  1:455734/139‑1 (MQ=255)                                                       agaCTCGCTGTGTAAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                                                                    <  2:385182/139‑1 (MQ=255)                                                         aCTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                                       <  2:454401/118‑1 (MQ=255)                                                         aCTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                                       >  1:454401/1‑118 (MQ=255)                                                                  gtgAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                                                         <  1:465282/139‑1 (MQ=255)                                                                  gtgAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                                                         <  1:319210/139‑1 (MQ=255)                                                                     aGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGc                                                      >  2:105607/1‑139 (MQ=255)                                                                       tttACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                                                    <  1:514938/139‑1 (MQ=255)                                                                             gTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAgcg                                                                                                   <  1:1196477/86‑1 (MQ=255)                                                                             gTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAgcg                                                                                                   >  2:1196477/1‑86 (MQ=255)                                                                                aaGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAATCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtg                                                                                    >  2:796907/1‑98 (MQ=255)                                                                                aaGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCCTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTTAGGGTtg                                                                                    <  1:796907/98‑1 (MQ=255)                                                                                          atatTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATa                                 >  1:860/1‑139 (MQ=255)                                                                                          atatTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATa                                 >  1:549249/1‑139 (MQ=255)                                                                                              tGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                                                    <  1:283129/116‑1 (MQ=255)                                                                                              tGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                                                    >  2:283129/1‑116 (MQ=255)                                                                                              tGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTg                             >  2:107652/1‑139 (MQ=255)                                                                                              tGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTg                             <  2:957682/139‑1 (MQ=255)                                                                                                  tttAAGGCTGTTCTCGCTTCCGAACCGAAGAAGTAACAGTTGTGCGGGGAGTAATCCCATAAGCGCTAACCTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACAAGAACTGGTCa                         >  1:799699/1‑139 (MQ=21)                                                                                                     aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGTGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCa                      >  1:863044/1‑139 (MQ=255)                                                                                                       ggCTGTTCTCGCTTACGAACCGAAGAAGTAACCGtt                                                                                                                           >  1:116125/1‑36 (MQ=255)                                                                                                       ggCTGTTCTCGCTTACGAACCGAAGAAGTAACCGtt                                                                                                                           <  2:116125/36‑1 (MQ=255)                                                                                                       ggCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAAt                    >  2:593965/1‑139 (MQ=255)                                                                                                        gCTGTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAAtt                   <  2:729736/139‑1 (MQ=255)                                                                                                           gTTCTCGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTg                <  2:316465/139‑1 (MQ=255)                                                                                                               tcGCTTACGAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGccc            <  1:372991/139‑1 (MQ=255)                                                                                                                       gAACCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGt    <  2:549249/139‑1 (MQ=255)                                                                                                                         aCCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTaa  <  1:105607/139‑1 (MQ=255)                                                                                                                         aCCGAAGAAGTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTaa  <  1:123640/139‑1 (MQ=255)                                                                                                                             aagaagTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc              >  1:478460/1‑123 (MQ=255)                                                                                                                             aagaagTAACCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc              <  2:478460/123‑1 (MQ=255)                                                                                                                                     |                                                                                                                               TGCGGCAGTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAA  >  NC_000913/607856‑608114

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