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%	212.9 / NA	65	*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 (23/42);  total (23/42)

TTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAGCCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACT  >  NC_000913/4352746‑4353011                                                                                                                                       |                                                                                                                                    ttCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCag                                                                                                                                 <  1:537531/139‑1 (MQ=255) ttCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCag                                                                                                                                 <  1:1133561/139‑1 (MQ=255) ttCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCag                                                                                                                                 <  2:647589/139‑1 (MQ=255)  tCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagg                                                                                                                                <  1:640551/139‑1 (MQ=255)  tCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagg                                                                                                                                <  2:1347219/139‑1 (MQ=255)   cGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               <  2:618034/139‑1 (MQ=255)   cGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTCTACAAGTAACCagga                                                                                                                               <  1:769474/139‑1 (MQ=255)        tAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGAt                                                                                                                          >  1:1469504/1‑139 (MQ=255)          aCGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATaa                                                                                                                        >  2:590766/1‑139 (MQ=255)             aCATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATaagaa                                                                                                                     <  2:1158067/139‑1 (MQ=255)             aCATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAACGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATaagaa                                                                                                                     >  1:1337542/1‑139 (MQ=255)              cATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATaagaag                                                                                                                    <  2:1416624/139‑1 (MQ=255)                    tAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATa                                                                                                              >  1:688201/1‑139 (MQ=255)                     aGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATaa                                                                                                             >  2:1192793/1‑139 (MQ=255)                         gCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATaaaaaa                                                                                                         <  2:207965/139‑1 (MQ=255)                            aGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               <  2:258727/114‑1 (MQ=255)                            aGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               >  1:258727/1‑114 (MQ=255)                            aGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAATAAGAACATAAAAAATgg                                                                                                      <  1:840949/139‑1 (MQ=255)                              tGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    <  1:1319769/139‑1 (MQ=255)                              tGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    <  2:688201/139‑1 (MQ=255)                              tGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    <  1:250185/139‑1 (MQ=255)                                              aTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATaaaaaa                                                                                                         <  2:1154350/118‑1 (MQ=255)                                              aTCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATaaaaaa                                                                                                         >  1:1154350/1‑118 (MQ=255)                                                             ttAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGAt                                                                                                                          >  2:1433460/1‑86 (MQ=255)                                                             ttAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGAt                                                                                                                          <  1:1433460/86‑1 (MQ=255)                                                              tAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGAt                                                                    <  2:421606/139‑1 (MQ=255)                                                                   ttACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGgtg                                                               <  2:1439890/139‑1 (MQ=255)                                                                   ttACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGgtg                                                               <  1:1262617/139‑1 (MQ=255)                                                                    tACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    >  1:339565/1‑101 (MQ=255)                                                                    tACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    <  2:339565/101‑1 (MQ=255)                                                                     aCCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACa                                                                                                                <  1:691024/88‑1 (MQ=255)                                                                     aCCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACa                                                                                                                >  2:691024/1‑88 (MQ=255)                                                                     aCCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGa                                                                     >  2:959335/1‑131 (MQ=255)                                                                     aCCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGa                                                                     <  1:959335/131‑1 (MQ=255)                                                                              tGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATcc                                                    <  2:1201198/139‑1 (MQ=255)                                                                                     gaggtCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGAtt                                                                   <  2:146388/113‑1 (MQ=255)                                                                                     gaggtCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGAtt                                                                   >  1:146388/5‑117 (MQ=255)                                                                                           aaCATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    <  2:15618/78‑1 (MQ=255)                                                                                           aaCATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    >  1:15618/1‑78 (MQ=255)                                                                                               ttATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGtt                                   <  1:1461070/139‑1 (MQ=255)                                                                                                       gTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               >  1:199869/1‑39 (MQ=255)                                                                                                       gTTGATGACTATTTCCGTCAGTTTTACAAGTAACCagga                                                                                                                               <  2:199869/39‑1 (MQ=255)                                                                                                       gTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTACAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTcac                           >  1:733677/1‑139 (MQ=255)                                                                                                       gTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGAt                                                                    >  1:590088/1‑98 (MQ=255)                                                                                                       gTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGAt                                                                    <  2:590088/98‑1 (MQ=255)                                                                                                       gTTGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTcac                           <  2:1448208/139‑1 (MQ=255)                                                                                                        ttGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTaa                                                                                         >  2:1337219/1‑76 (MQ=255)                                                                                                        ttGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTaa                                                                                         <  1:1337219/76‑1 (MQ=255)                                                                                                        ttGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGa                                                                     <  1:415554/96‑1 (MQ=255)                                                                                                        ttGATGACTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGa                                                                     >  2:415554/1‑96 (MQ=255)                                                                                                          gatgaCTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    <  1:648603/63‑1 (MQ=255)                                                                                                          gatgaCTATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCa                                                                                                    >  2:648603/1‑63 (MQ=255)                                                                                                                tATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGgt                                                                >  2:170600/1‑93 (MQ=255)                                                                                                                tATTTCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGgt                                                                <  1:170600/93‑1 (MQ=255)                                                                                                                    tCCGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAAcaca              <  1:1192793/139‑1 (MQ=255)                                                                                                                      cGTCAGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACata            <  2:1330168/139‑1 (MQ=255)                                                                                                                          aGTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCg        >  2:959223/1‑139 (MQ=255)                                                                                                                           gTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGt       <  1:453161/139‑1 (MQ=255)                                                                                                                           gTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGt       <  2:1429011/139‑1 (MQ=255)                                                                                                                           gTTTTACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGt       <  1:959223/139‑1 (MQ=255)                                                                                                                               tACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTcaca                          >  2:1260710/1‑116 (MQ=255)                                                                                                                               tACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTcaca                          <  1:1260710/116‑1 (MQ=255)                                                                                                                               tACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACata            <  2:1224425/130‑1 (MQ=255)                                                                                                                               tACAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACata            >  1:1224425/1‑130 (MQ=255)                                                                                                                                aCAAGTAACCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACt  <  1:590766/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                    TTCGCAATAACGACATCGTTAGCCGCAAGTGAAATCGAAGATTTAATCAGGCTCAAATGCTTAGACTTACCGGATATTGATTTTGACCTCAACATTATGACAGTTGATGACTATTTCCGTCAGTTTTACAAGTAGCCAGGAAGGATAAGAAGAACATAAAAAATGGCACTATTCTCTAAAATATTAATTTTTTATGTGATTGGTGTGAACATATCCTTTGTCATTATCTGGTTTATCTCACATGAGAAAACACATATTCGTTTACT  >  NC_000913/4352746‑4353011

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