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%	97.0 / NA	31	*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 (21/10);  total (21/10)

TCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGC  >  NC_000913/4352959‑4353213                                                                                                                            |                                                                                                                                     tCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACCTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGt                                                                                                                       >  2:699792/1‑139 (MQ=255) tCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTg                                                                                                                      >  1:483722/1‑139 (MQ=255) tCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTACCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTg                                                                                                                      >  1:580743/1‑139 (MQ=255)            tATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGt                                                                                                           >  1:401744/1‑139 (MQ=255)             aTCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGtt                                                                                                          >  1:38731/1‑139 (MQ=255)                  gTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGt                                                                                                     >  2:776779/1‑139 (MQ=255)                  gTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGt                                                                                                     >  1:261108/1‑139 (MQ=255)                      ctctcACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGaaa                                                                                                 >  1:499762/2‑139 (MQ=255)                                agaAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGtct                                                                                       >  1:35752/1‑139 (MQ=255)                                                     ttAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCg                                                                  <  1:96514/139‑1 (MQ=255)                                                       aGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGa                                                                <  2:401744/139‑1 (MQ=255)                                                         tGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGc                                                                                                               >  1:513843/1‑90 (MQ=255)                                                         tGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGc                                                                                                               <  2:513843/90‑1 (MQ=255)                                                         tGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAAc                                                              >  1:439959/1‑139 (MQ=255)                                                         tGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAAc                                                              >  2:144325/1‑139 (MQ=255)                                                              tCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGa                                                         >  1:130813/1‑139 (MQ=255)                                                                   gTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGtct                                                                                       >  2:13391/1‑104 (MQ=255)                                                                   gTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGtct                                                                                       <  1:13391/104‑1 (MQ=255)                                                                     cGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGc                                                  <  1:776779/139‑1 (MQ=255)                                                                             ccTGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCt                                          <  2:261108/139‑1 (MQ=255)                                                                               tGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGa                                        <  1:144325/139‑1 (MQ=255)                                                                                  ccAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTg                                     >  2:20394/1‑139 (MQ=255)                                                                                  ccAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTg                                     >  2:215529/1‑139 (MQ=255)                                                                                  ccAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTg                                     >  2:9068/1‑139 (MQ=255)                                                                                  ccAATGAGTCTACCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTg                                     <  2:580743/139‑1 (MQ=255)                                                                                           cTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCa                            <  1:688539/139‑1 (MQ=255)                                                                                                   ggCATTA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTg                    >  2:190551/1‑139 (MQ=255)                                                                                                        tA‑CTTTTTTCTCTCTTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAatt               >  1:780913/1‑139 (MQ=255)                                                                                                                      ctTTTAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGgc  >  1:590023/1‑139 (MQ=255)                                                                                                                        tttAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGt                                                                                         <  1:607313/50‑1 (MQ=255)                                                                                                                        tttAAGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGt                                                                                         >  2:607313/1‑50 (MQ=255)                                                                                                                            |                                                                                                                                     TCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACTTAGTGCATTCCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTA‑CTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGC  >  NC_000913/4352959‑4353213

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