Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,353,081	G→A	*58* (TAG→TAA)	ghoT →	toxin of GhoTS toxin‑antitoxin pair, membrane‑lytic protein, stimulator of persister cell formation

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,353,081	0	G	A	100.0%	120.5 / NA	38	*58* (TAG→TAA)	ghoT	toxin of GhoTS toxin‑antitoxin pair, membrane‑lytic protein, stimulator of persister cell formation
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/18);  total (20/18)

ACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCA  >  NC_000913/4352954‑4353216                                                                                                                                |                                                                                                                                        aCATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGtt                                                                                                                              <  1:892708/139‑1 (MQ=255)     atCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGt                                                                                                                          >  1:97580/1‑139 (MQ=255)         tttGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGAc                                                                                                                      <  1:356513/139‑1 (MQ=255)          ttGTCATTATCTGGTTTATCTCACATGAGAAAACACATAATCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGAcc                                                                                                                     >  1:61654/1‑139 (MQ=255)                       gTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGt                                                                                                        >  2:113941/1‑139 (MQ=255)                                 caTGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGc                                                                                              >  1:774856/1‑139 (MQ=255)                                                     tttACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGAt                                                                          >  1:496569/1‑139 (MQ=255)                                                            aGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAagag                                                                                                                                      >  1:964659/1‑72 (MQ=255)                                                            aGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAagag                                                                                                                                      <  2:964659/72‑1 (MQ=255)                                                              tGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAagag                                                                                                                                      >  1:1337209/1‑70 (MQ=255)                                                              tGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAagag                                                                                                                                      <  2:1337209/70‑1 (MQ=255)                                                              tGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTgg                                                                                                                        >  1:885899/1‑84 (MQ=255)                                                              tGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTgg                                                                                                                        <  2:885899/84‑1 (MQ=255)                                                              tGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAAc                                                                 >  1:855273/1‑139 (MQ=255)                                                               gCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACg                                                                <  2:774856/139‑1 (MQ=255)                                                                  ttCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGtct                                                                                          >  2:1371164/1‑110 (MQ=255)                                                                  ttCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGtct                                                                                          <  1:1371164/110‑1 (MQ=255)                                                                   tCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGt                                                                                            <  2:589774/107‑1 (MQ=255)                                                                   tCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGt                                                                                            >  1:589774/1‑107 (MQ=255)                                                                        gTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAagag                                                                                                                                      <  2:106421/60‑1 (MQ=255)                                                                        gTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAagag                                                                                                                                      >  1:106421/1‑60 (MQ=255)                                                                        gTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTg                                                                                                                                   >  1:949605/1‑63 (MQ=255)                                                                        gTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTg                                                                                                                                   <  2:949605/63‑1 (MQ=255)                                                                        gTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTg                                                                                                       >  1:318627/1‑91 (MQ=255)                                                                        gTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTg                                                                                                       <  2:318627/91‑1 (MQ=255)                                                                        gTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTg                                                                       >  1:589178/1‑123 (MQ=255)                                                                        gTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGCACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGaa                                                       >  1:896606/1‑139 (MQ=255)                                                                        gTCGGAATAACCTAGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTg                                                                       <  2:589178/123‑1 (MQ=255)                                                                               tAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCg                                                <  2:496569/139‑1 (MQ=255)                                                                               tAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCg                                                <  2:97580/139‑1 (MQ=255)                                                                                  ccTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTgg                                                                                                                        <  2:143633/64‑1 (MQ=255)                                                                                  ccTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTgg                                                                                                                        >  1:143633/1‑64 (MQ=255)                                                                                       ccAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTg                                                                                                       >  1:652699/1‑76 (MQ=255)                                                                                       ccAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTg                                                                                                       <  2:652699/76‑1 (MQ=255)                                                                                       ccAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACTAAGCAGTCGCCTGACTg                                        <  2:855273/139‑1 (MQ=255)                                                                                       ccAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTg                                        >  1:56872/1‑139 (MQ=255)                                                                                           tGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTc                                                 >  2:1138876/1‑126 (MQ=255)                                                                                           tGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTc                                                 <  1:1138876/126‑1 (MQ=255)                                                                                                       tgGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGt                                                                                            >  1:94257/1‑71 (MQ=255)                                                                                                       tgGCATTACTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGt                                                                                            <  2:94257/71‑1 (MQ=255)                                                                                                                         tctTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACg                                                                                                  <  2:594107/47‑1 (MQ=255)                                                                                                                         tctTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACg                                                                                                  >  1:594107/1‑47 (MQ=255)                                                                                                                         tctTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTggg      >  1:634039/1‑139 (MQ=255)                                                                                                                             ttAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCa  <  2:634039/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                        ACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCA  >  NC_000913/4352954‑4353216

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