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%	134.5 / NA	42	*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 (15/27);  total (15/27)

TCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAGCCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACT  >  NC_000913/4352747‑4353011                                                                                                                                      |                                                                                                                                    tCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagg                                                                                                                                <  2:219093/139‑1 (MQ=255)  cGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               <  2:958288/139‑1 (MQ=255)  cGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               <  1:352411/139‑1 (MQ=255)  cGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               <  1:1405497/139‑1 (MQ=255)       tAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGAGATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGAt                                                                                                                          <  2:820967/139‑1 (MQ=255)                   tAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATa                                                                                                              <  2:202446/139‑1 (MQ=255)                     gCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATACTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATaaa                                                                                                            <  1:760316/139‑1 (MQ=255)                           aGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATgg                                                                                                      >  1:397178/1‑139 (MQ=255)                             tGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    <  2:486685/139‑1 (MQ=255)                             tGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    <  1:920988/139‑1 (MQ=255)                                                  gCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    >  1:230388/1‑118 (MQ=255)                                                  gCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    <  2:230388/118‑1 (MQ=255)                                                            ttAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGa                                                                     <  2:397178/139‑1 (MQ=255)                                                                    aCCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               <  2:722070/73‑1 (MQ=255)                                                                    aCCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               >  1:722070/1‑73 (MQ=255)                                                                    aCCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATaag                                                                                                                       >  1:1061248/1‑81 (MQ=255)                                                                    aCCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATaag                                                                                                                       <  2:1061248/81‑1 (MQ=255)                                                                          tatTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTa                                                                                                 <  1:198783/97‑1 (MQ=255)                                                                          tatTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTa                                                                                                 >  2:198783/1‑97 (MQ=255)                                                                               aTTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    >  2:541848/1‑89 (MQ=255)                                                                               aTTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    <  1:541848/89‑1 (MQ=255)                                                                                   tGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGa                                                                     >  1:236611/1‑116 (MQ=255)                                                                                   tGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGa                                                                     <  2:236611/116‑1 (MQ=255)                                                                                     aCCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAAttttt                                                                            <  1:1069201/107‑1 (MQ=255)                                                                                     aCCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAAttttt                                                                            >  2:1069201/1‑107 (MQ=255)                                                                                         cAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATaagaa                                                                                                                     >  1:757806/1‑62 (MQ=255)                                                                                         cAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATaagaa                                                                                                                     <  2:757806/62‑1 (MQ=255)                                                                                                      gTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               >  2:51206/1‑39 (MQ=255)                                                                                                      gTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               <  1:51206/39‑1 (MQ=255)                                                                                                      gTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    <  2:1381257/66‑1 (MQ=255)                                                                                                      gTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    >  1:1381257/1‑66 (MQ=255)                                                                                                      gTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTcac                           >  2:412726/1‑139 (MQ=255)                                                                                                       ttGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    <  2:805272/65‑1 (MQ=255)                                                                                                       ttGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    >  1:805272/1‑65 (MQ=255)                                                                                                                  ttCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    >  1:1037176/1‑54 (MQ=255)                                                                                                                  ttCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    <  2:1037176/54‑1 (MQ=255)                                                                                                                     cGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACata            <  1:412726/139‑1 (MQ=255)                                                                                                                     cGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACata            <  2:1256837/139‑1 (MQ=255)                                                                                                                         aGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTGTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTGACATGAGAAAACACATATTCg        >  2:833389/1‑139 (MQ=255)                                                                                                                               aCAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACt  <  1:599644/139‑1 (MQ=255)                                                                                                                               aCAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACt  <  1:72328/139‑1 (MQ=255)                                                                                                                               aCAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACt  <  1:833389/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                    TCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAGCCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACT  >  NC_000913/4352747‑4353011

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