Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,352,791	A→G	I70V (ATC→GTC)	ghoS →	antitoxin of the GhoTS toxin‑antitoxin system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,352,791	0	A	G	100.0%	185.1 / NA	51	I70V (ATC→GTC)	ghoS	antitoxin of the GhoTS toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/26);  total (25/26)

TTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAGCCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATT  >  NC_000913/4352656‑4352918                                                                                                                                        |                                                                                                                                ttaagattaagaTCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAg                                                                                                                             >  1:102534/1‑140 (MQ=255)      attaagatCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCa                                                                                                                        >  2:272042/1‑140 (MQ=255)      attaagatCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCa                                                                                                                        <  1:198653/140‑1 (MQ=255)      attaagatCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCa                                                                                                                        <  1:272042/140‑1 (MQ=255)           gatCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGc                                                                                                                   <  2:308296/140‑1 (MQ=255)                     gTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTAc                                                                                                         >  2:158516/1‑140 (MQ=255)                      tatgatatgaATTTCTCCTCTATCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCt                                                                                                                  >  1:164711/1‑130 (MQ=255)                      tatgatatgaATTTCTCCTCTATCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCt                                                                                                                  <  2:164711/130‑1 (MQ=255)                        tgatatgaATTTCTCCTCTATCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCgg                                                                                                      >  1:40391/1‑140 (MQ=255)                            atgaATTTCTCCTCTATCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGAc                                                                                                             <  2:205662/129‑1 (MQ=255)                            atgaATTTCTCCTCTATCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGAc                                                                                                             >  1:205662/1‑129 (MQ=255)                                    ctcctcTATCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTg                                                                                          >  2:274139/1‑140 (MQ=255)                                     tcctcTATCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGa                                                                                         >  1:139785/1‑140 (MQ=255)                                     tcctcTATCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGa                                                                                         >  1:85497/1‑140 (MQ=255)                                     tcctcTATCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGa                                                                                         >  1:335415/1‑140 (MQ=255)                                         ctATCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTc                                                                                     >  1:179362/1‑140 (MQ=255)                                         ctATCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTc                                                                                     <  2:166127/140‑1 (MQ=255)                                         ctATCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTc                                                                                     <  1:157730/140‑1 (MQ=255)                                           aTCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCaa                                                                                   <  2:335415/140‑1 (MQ=255)                                                 gCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCaa                                                                                   >  1:33102/1‑134 (MQ=255)                                                 gCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCaa                                                                                   <  2:33102/134‑1 (MQ=255)                                                    gatgaATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTc                                                                                     >  1:321640/1‑129 (MQ=255)                                                    gatgaATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTc                                                                                     <  2:321640/129‑1 (MQ=255)                                                    gatgaATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAAc                                                                                  >  1:242417/1‑132 (MQ=255)                                                    gatgaATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAAc                                                                                  <  2:242417/132‑1 (MQ=255)                                                      tgaATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAg                                                                        >  2:151415/1‑140 (MQ=255)                                                             gggATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTtgatga                                                                 <  2:220713/140‑1 (MQ=255)                                                               gATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACt                                                               >  2:91034/1‑140 (MQ=255)                                                               gATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACt                                                               >  2:69090/1‑140 (MQ=255)                                                               gATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACt                                                               >  2:337225/1‑140 (MQ=255)                                                                    cACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTtgatga                                                                 <  2:179787/133‑1 (MQ=255)                                                                    cACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTtgatga                                                                 >  1:179787/1‑133 (MQ=255)                                                                    cACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTtg                                                                     >  2:18639/1‑129 (MQ=255)                                                                    cACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTtg                                                                     <  1:18639/129‑1 (MQ=255)                                                                           aTTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGt                                                   >  2:336835/1‑140 (MQ=255)                                                                             tgaatgaaAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGttt                                                 <  1:91034/140‑1 (MQ=255)                                                                             tgaatgaaAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGttt                                                 <  1:69090/140‑1 (MQ=255)                                                                               aatgaaAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTa                                               >  2:291441/1‑140 (MQ=255)                                                                                    aaaCTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGt                                          <  2:131794/140‑1 (MQ=255)                                                                                          cTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAGCCag                                    <  2:139785/140‑1 (MQ=255)                                                                                          cTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAGCCag                                    <  2:85497/140‑1 (MQ=255)                                                                                           ttCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAGCCagg                                   <  1:274139/140‑1 (MQ=255)                                                                                           ttCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAGCCagg                                   >  2:170425/1‑140 (MQ=255)                                                                                             cGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAGCCaggaa                                 <  1:353858/140‑1 (MQ=255)                                                                                                  tAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTAc                                              <  1:41766/122‑1 (MQ=255)                                                                                                  tAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTAc                                              >  2:41766/1‑122 (MQ=255)                                                                                                    aCGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAGCCAGGAAGGATaag                          <  2:179362/140‑1 (MQ=255)                                                                                                                       gTGAAATCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAGCCAGGAAGGATAAGAAGAACATAAAAAATGGCa       <  2:151496/140‑1 (MQ=255)                                                                                                                          aaaTCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAGCCAGGAAGGATAAGAAGAACATAAAAAATGGCACTa    <  1:170425/140‑1 (MQ=255)                                                                                                                          aaaTCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAGCCAGGAAGGATAAGAAGAACATAAAAAATGGCACTa    <  2:23015/140‑1 (MQ=255)                                                                                                                            aTCGAAGATTTAGTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAGCCAGGAAGGATAAGAAGAACATAAAAAATGGCACTAtt  >  1:61530/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                TTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATGGGATACCACGACAATTGAATGAAAACTCCTTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAGCCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATT  >  NC_000913/4352656‑4352918

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