Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,352,880	G→A	*99* (TAG→TAA)	ghoS →	antitoxin of GhoTS toxin‑antitoxin pair, endonuclease for ghoT mRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,352,880	0	G	A	100.0%	151.7 / NA	48	*99* (TAG→TAA)	ghoS	antitoxin of GhoTS toxin‑antitoxin pair, endonuclease for ghoT mRNA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/30);  total (18/30)

TCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAGCCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACT  >  NC_000913/4352747‑4353011                                                                                                                                      |                                                                                                                                    tCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGGCTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagg                                                                                                                                <  1:782797/139‑1 (MQ=255) tCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTAGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagg                                                                                                                                <  2:850125/139‑1 (MQ=255)  cGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               <  2:354600/139‑1 (MQ=255)  cGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               <  2:458365/139‑1 (MQ=255)  cGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               <  1:443561/139‑1 (MQ=255)  cGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCag                                                                                                                                 <  2:263807/137‑1 (MQ=255)  cGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCag                                                                                                                                 >  1:263807/1‑137 (MQ=255)  cGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGACTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               <  2:888092/139‑1 (MQ=255)      aTAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCaggaagga                                                                                                                           >  2:1068870/1‑139 (MQ=255)       tAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGAt                                                                                                                          <  1:729545/139‑1 (MQ=255)       tAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGAt                                                                                                                          <  1:1068870/139‑1 (MQ=255)                 gTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACa                                                                                                                >  2:130571/1‑139 (MQ=255)                  ttAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               <  2:522420/123‑1 (MQ=255)                  ttAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               >  1:522420/1‑123 (MQ=255)                          aaGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCa                                                                                                                                  >  2:945603/1‑112 (MQ=255)                          aaGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAAGAAGGATAAGAAGAACATAAAAAATg                                                                                                       >  1:774408/1‑139 (MQ=255)                          aaGTGAAATCGAAGAGTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCa                                                                                                                                  <  1:945603/112‑1 (MQ=255)                             tGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    <  1:460810/139‑1 (MQ=255)                             tGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    <  1:1216671/139‑1 (MQ=255)                                        aTTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAGCATTATGACAGTTGATGACTATTGCCGTCAGTTTTACAATTAACCAGGAATGATAAGAAGAACATAAAAAATGGCACTATTCTCTaa                                                                                         <  2:774408/139‑1 (MQ=255)                                                            ttAGTCTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGa                                                                     <  1:408976/139‑1 (MQ=255)                                                            ttAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGa                                                                     <  1:197259/139‑1 (MQ=255)                                                            ttAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGa                                                                     <  1:130571/139‑1 (MQ=255)                                                                  ttACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGgtg                                                               <  1:29964/139‑1 (MQ=255)                                                                  ttACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGgtg                                                               <  1:469747/139‑1 (MQ=255)                                                                   tACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCag                                                                                                                                 <  2:983481/72‑1 (MQ=255)                                                                   tACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCag                                                                                                                                 >  1:983481/1‑72 (MQ=255)                                                                    aCCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagg                                                                                                                                <  2:28454/72‑1 (MQ=255)                                                                    aCCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagg                                                                                                                                >  1:28454/1‑72 (MQ=255)                                                                    aCCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATaaga                                                                                                                      >  2:336478/1‑82 (MQ=255)                                                                    aCCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATaaga                                                                                                                      <  1:336478/82‑1 (MQ=255)                                                                    aCCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTAtt                                                                                               <  1:743735/105‑1 (MQ=255)                                                                    aCCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGgtg                                                               >  1:90442/1‑137 (MQ=255)                                                                    aCCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGgtg                                                               <  2:90442/137‑1 (MQ=255)                                                                    aCCGGATATTGATTTTGACCTCAACATTATGACAGTTGATAACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTAtt                                                                                               >  2:743735/1‑105 (MQ=255)                                                                              gATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCt                                                   <  2:617531/139‑1 (MQ=255)                                                                                  ttGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               <  2:1190448/59‑1 (MQ=255)                                                                                  ttGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               >  1:1190448/1‑59 (MQ=255)                                                                                   tGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               <  1:499791/58‑1 (MQ=255)                                                                                   tGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               >  2:499791/1‑58 (MQ=255)                                                                                          aaCATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGa                                                                     <  2:1005844/109‑1 (MQ=255)                                                                                          aaCATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGa                                                                     >  1:1005844/1‑109 (MQ=255)                                                                                                      gTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACa                                                         <  2:1002743/109‑1 (MQ=255)                                                                                                      gTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACa                                                         >  1:1002743/1‑109 (MQ=255)                                                                                                         gatgaCTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATg                        >  2:379889/1‑139 (MQ=255)                                                                                                         gatgaCTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATg                        >  2:897763/1‑139 (MQ=255)                                                                                                           tgaCTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATgag                      >  1:275007/1‑139 (MQ=255)                                                                                                                               aCAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTGATCTCACATGAGAAAACACATATTCGTTTACt  <  1:725040/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                    TCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAGCCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACT  >  NC_000913/4352747‑4353011

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